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StlllUN  I.  INTRODUCTION 


1.  INTRODUCTION  AND  DEFINITION  OF  TEKKO 

An  understanding  of  the  role  of  negative  iors  in  low  pressure  gas 
discharges  is  important  in  several  fields  such  as  electric  discharge 
lasers,  high  power  switches,  and  plasma  deposition  of  thin  films. 
Discharges  such  as  these  which  contain  negative  ion  densities  in 
addition  to  the  usual  electron  and  positive  ion  densities  are  termed 
electronegative  discharges.  They  are  so  named  since  the  attaching  gas 
yielding  negative  ions  is  referred  to  as  an  electronegative  gas. 
Discharges  containing  only  electrons  and  positive  ions  are  termer 
electropositive  discharges. 

Both  types  of  discharges  are  characterized  by  the  distinctive 
luminous  and  dark  regions  shown  in  Figure  1.  The  actual  position  of 
the  various  regions  depends  upon  the  gas  or  gas  mixture,  pressure,  and 
current  or  voltage  across  the  discharge.  Starting  at  the  cathode  there 
is  sometimes  visible  a  very  narrow  dark  space  termed  Aston's  dark  space 
followed  by  a  feeble  layer  of  light  called  the  cathode  glow.  However, 
the  term  cathode  glow  often  includes  the  tota1  region  from  the  cathode 
t.o  the  beginning  of  the  negative  glow.  The  region  between  the  cathode 
glow  ar.a  negative  clow  is  ret  totally  dark,  but  only  appears  so  to  the 
eye  in  comparison  to  its  luminous  neighboring  reqiors  and  is  termed  the 
cathode  dark  space.  The  negative  ylow  is  probably  the  most  luminous 
region  of  the  discharge.  A  sharp  luminous  boundary  exists  with  the 
cathode  nonlumincus  space,  but  the  bounoary  on  the  anode  side  is 
diffuse  and  gradually  grows  dimmer  as  it  merges  into  the  relatively 
dark  laraday  dark  space.  This  is  normally  followeo  by  a  long  luminous 


1 


Figure  1  Discharge  Parameters  as  a  Function  of  Discharge  Length 


positive  column  region  extending  up  to  the  anode  glow.  The  anoct  glow 
appears  adjacent  to  the  anode  surface.  Depending  on  the  distance 
between  the  electrodes,  only  Aston's  dark  space,  the  cathode  glow,  and 
the  cathode  dark  space  are  necessary  lor  a  discharge  to  exist. 

It  has  been  observed  that  the  presence  of  negative  ions  can 
spatially  alter  the  above  regions  ana  in  some  geces  such  as  iodine, 
negative  ions  can  become  the  dominant  negative  charge  carrier  through 
most  of  the  discharge  regions.  Formation  of  negative  ions  in  the 
positive  column  of  the  discharge  often  results  in  striations  which  can 
cause  fluctuations  in  the  light  output  of  the  positive  column.  The 
most  undesirable  effect  of  negative  ions  is  in  causing  severe  radial 
construction  of  the  positive  column  limiting  the  efficient  usage  cf  the 
available  container  volume.  Up  to  now,  however,  the  exact  role  nega¬ 
tive  ions  play  in  the  cathode  fall  region  of  a  discharge  was  unknown. 
The  literature  frequently  cites  experimental  evidence  indicating  that 
the  cathode  fall  region  contracts  longitudinally  with  the  addition  of 
even  a  trace  of  an  attaching  gas.  This  was  interpreted  to  mean  that 
the  electronegative  gas  can  have  a  significant  impact  in  the  region. 

The  original  question  which  initiated  this  study  was  if  the  formation 
of  negative  ions  in  the  cathode  fall  led  to  the  longitudinal  contrac¬ 
tion  of  the  cathode  fall  region. 

before  proceeding  further,  definitions  of  equilibrium  and  non- 
equilibrium  are  necessary  to  alleviate  any  confusion.  ‘  Fouil  ibriuni'  in 
high  pressure  discharges  is  termed  local  thermodynamic  equil.brium 
iLTE)  and  means  that  all  states  of  the  charged  ano  neutral  species  can 
be  described  by  one  temperature .  In  contrast,  at  iow  pressures  the 
term  'eaui  1  ibriuni' ,  more  properly  called  'field  equilibrium',  means 


that  the  electron  energy  distribution  function  is  characterized  by  the 
local  value  of  E/N  (E  is  the  electric  field,  N  is  the  gas  number 
density)  and  the  species  in  the  volume  element  ot  interest.  The 
electron  distribution  function  need  not  be  Maxwellian,  as  long  as 
it  can  be  defined  in  terms  of  E/N.  Nonequilibrium  at  low  pressures 
occurs  when  the  above  criterion  for  field  equilibrium  does  not  apply. 
For  example,  electrons  are  not  in  equilibrium  in  large  spatial  gradi¬ 
ents  where  the  mean  velocity  of  the  electrons  would  not  be  equal  to  the 
mean  velocity  of  electrons  in  an  equivalent,  but  spatially  uniform 
field. 

Basic  analytical  analysis  (References  ?0,  21,  60,  79,  86)  now 
indicates  that  the  electrons  are  not  in  eouilibrium  with  the  electric 
field  in  the  cathode  fall  region.  Since  the  electrons  are  being  accel¬ 
erated  across  the  cathode  fall  region,  they  are  gaining  energy  from  the 
field  faster  than  they  lose  energy  to  elastic  scattering  or  inelastic 
processes  in  the  gas.  The  electrons  leave  the  cathode  at  low  energies 
of  few  electron  volts  where  the  field  is  actually  the  highest.  The 
electrons  approach  the  negative  glow  with  a  multi-peaked  distribution 
function  as  the  field  is  approaching  a  minimum.  Since  the  electrons 
are  not  in  equilibrium  with  the  field,  their  multiplication  cannot  be 
described  by  the  usual  Townsend  ionization  coefficient.  This  coeffi¬ 
cient  is  measured  in  uniform  fields  where  the  energy  gained  by  the 
electrons  from  the  field  is  balanced  by  their  energy  losses  to  various 
elastic  and  inelastic  processes  with  the  gas  molecules.  By  analogy, 
the  equilibrium  attachment  coefficient  also  is  not  expected  to  cor¬ 
rectly  describe  the  formation  of  negative  ions  acuss  the  catb  de  fall 
region.  The  attachment  coetficient  is  expected  to  be  large  adjacent  to 


the  cathode  and  then  to  decrease  as  the  electrons  are  accelerated  to 
higher  energies.  The  formation  cf  negative  ions  in  this  nenequilibrium 
sheath  region  has  not  been  studied  either  analytically  or  experimental¬ 
ly.  Concentrations  of  negative  ions  could  alter  the  space  charge  in 
the  sheath  region  and  thus  affect  the  voltage  drop  across  the  cathode 
tall  region.  This  phenomena  of  noriequi  1  ibrium  in  the  cathode  fall 
region  requires  a  'nonequilibrium  theory',  such  as  calculating  the 
electron  energy  distribution  function  as  ?  function  of  distance  in 
order  to  properly  model  the  electron  number  density  and  ionization  and 
attachment  coefficients.  However,  'equilibrium  theory'  has  often  been 
applied  to  Gtscribe  the  electric  field,  voltage  drop,  and  electron 
current  because  its  results  have  agreed  fortuitously  with  the  macro¬ 
scopic  parameters  ihat.  have  been  able  to  be  measured  across  the  cathode 
fall  region. 

£.  OBJECTIVES  OF  THIS  STUDY 

This  analytical  study  investigu ‘ the  role  small  concentrations 
of  an  attaching  gas  play  in  the  cathode  fall  region  of  a  glow  dis¬ 
charge.  To  do  this,  c  new  equilibrium  model  based  on  current  continu¬ 
ity  equations  and  Poisson's  equations  was  developed.  Additionally,  a 
nonequilibrium  model  based  on  the  Boltzmann  transport  equation  and 
Frisson's  equation  was  modified  to  be  self-consistent,  include  negative 
ions  and  be  more  flexible  The  electron  energy  distribution  function 
was  calculated  as  a  function  of  distance  through  the  cathode  fa  1 1 
region  and  the  effects  of  negative  ion  formation  were  examined.  The 
first  and  most  important  objective  was  to  examine  'he  effects  n  the 
(athode  fa1!  region  resuiting  from  the  addition  cf  small  amounts  of  an 
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attaching  gas  such  as  HC1  to  a  rare  gas  discharge.  The  second  objec¬ 
tive  was  to  model  the  nonequi 1 ibrium  behavior  of  the  three  charged 
particle  densities  through  the  cathode  fall  region.  The  third  objec¬ 
tive  was  to  compare  and  contrast  the  equilibrium  and  nonequilibrium 
approaches  and  determine  their  regions  of  applicability. 

3.  OVERVIEW  OF  CONTENTS  AND  RESULTS 

Section  I  describes  the  purpose  and  objective  of  this  study, 
defines  the  regions  of  a  discharge,  and  outlines  the  contents  of  the 
remaining  chapters. 

In  Section  II,  the  sufficient  condition  for  the  existence  of  a 
discharge  is  described  in  terms  of  electron  multiplication  in  the 
cathode  fall  region.  A  historical  review  is  presented  on  how  trace 
amounts  of  negative  ions  affect  an  electropositive  discharge.  A 
qualitative  theory  is  described  which  simplistically  describes  the 
result  of  negative  ion  formation  in  glow  discharges. 

Section  III  reviews  the  previous  equilibrium  theories  of  the 
cathode  fall  region,  primarily  for  electropositive  gases.  These 
equilibrium  theories  are  based  on  continuity  type  equations  coupled 
with  Poisson's  equation.  Boundary  conditions  are  also  discussed.  A 
new  equilibrium  model  of  the  cathode  fall  region  allowing  for  negative 
ion  formation  is  described  and  numerical  results  are  presented. 

In  comparison  to  the  equilibrium  approaches  described  in  Section 
III,  Section  IV  reviews  the  previous  ncnequilihrium  theories  of  the 
cathode  fall  region  based  on  calculating  the  electron  distribution  as  a 
function  of  position  through  the  cathode  fall  region.  A  numerical 
technique  for  calculating  the  electron  distribution  as  a  function  of 


position  in  the  cathode  fall  is  described.  This  study  extends  this 
technique  to  include  attachment,  the  presence  of  negative  ions  in  the 
relevant  equations,  gas  mixtures  with  a  variable  number  of  collision 
processes,  and  the  calculation  of  the  Townsend  ionization  and 
attachment  rates  from  the  distribution  function.  A  convergent  tech¬ 
nique  was  developed  which  leads  to  a  self  consistent  solution  between 
Poisson's  equation  and  the  electron  distribution  function.  The  boun¬ 
dary  conditions,  electron  and  ion  kinetics,  and  the  self  consistency  of 
the  solution  are  discussed.  Characteristic  cathode  fall  parameters  are 
calculated  as  a  function  of  distance  and  percentage  ot  attaching 
species. 

Section  V  compares  the  equilibrium  theory  with  the  nonequilibrium 
theory  and  discusses  the  discrepancies  between  the  two  theories  in  the 
cathode  fall  region.  The  ancGe  fall  region  is  also  discussed. 

Finally,  Section  VI  summarizes  the  new  results  of  this  study  and 
presents  some  considerations  for  future  study. 
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QUALITATIVE  THEORY 


1.  DESCRIPTION  OF  THE  CATHODE  FALL  REGION 

The  regions  of  a  discharge  were  ioentificc  in  Figure  1.  Of  these, 
the  cathode  fall  region  is  the  most  important  segment  of  the  discharge 
in  that  the  largest  potential  drop  per  unit  distance  usually  occurs 
there.  However,  it  is  probably  the  region  least  understood  theoret¬ 
ically.  It  is  the  region  of  the  discharge  closest  to  the  cathode  and 
consists  of  the  three  subregions  identified  in  the  previous  chapter  as 
the  primary  dark  space  or  Aston's  dark  space,  the  cathode  glow,  and  the 
cathode  dark  space.  The  properties  of  the  cathode  fall  region  are 
almost  independent  of  the  rest  of  the  discharge,  since  a  discharge  can 
exist  without  a  positive  column,  Faraday  dark  space,  or  negative  glow. 

A  discharge  cannot  ex^st,  however,  without  a  cathode  dark  space  (Refer¬ 
ence  6).  Specifically,  the  length  of  the  cathode  dark  space,  the 
voltage  drop  across  this  distance,  and  the  discharge  current  do  not 
depend  on  any  other  external  physical  parameters  such  as  the  length  of 
the  discharge  (providing  the  anode  is  not  so  small  as  to  inhibit  the 
axial  current).  The  cathode  dark  space  is  the  longest  of  the  three 
cathode  fall  regions. 

The  essential  processes  in  the  cathode  fall  region  are  as  lollows. 
i fit  purpose  of  the  cathode  is  to  provide  electrons  to  maintain  the 
discharge.  Several  processes  exist  by  which  the  catfiode  can  emit 
electrons.  Direct  processes  include  thermionic  err  ssion  and  t  eld 
emission.  For  applied  voltages  up  to  several  thousand  volts  these 
processes  can  he  ignored.  Indirect  emission  of  electrons,  usually 


called  secondary  emission,  results  from  bombardment  of  positive  ions, 
photons,  excited  and  neutral  molecules  on  the  cathode.  Due  to  the  low 
collision  rate  of  excited  ano  neutral  molecules  with  the  cathode  in 
most  discharges,  the  physical  processes  contributing  most  to  the 
emission  of  electrons  from  the  cathode  are  ion  and  photon  bombardment. 
Ihe  electrons  produced  by  secondary  emission  i » •: ■  accelerated  across  the 
cathode  dark  space,  whereupon  they  begin  to  excite  and  ionize  the  gas 
creating  new  electrons  and  positive  ions.  The  positive  ions  in  turn 
are  accelerated  back  towards  the  cathode  where  they  produce  new  elec¬ 
trons.  For  every  (M~l )  electrons  produced  through  ionization  in  the 
gas,  y  electrons  are  released  at  the  cathode.  The  maintenance  condi¬ 
tion  for  a  self  sustained  discharge  becomes 

Y (M-l )  =  1  or  M=1  +  |  (1) 

The  multiplication  M  of  electrons  depends  upon  the  nature  and  pressure 
of  the  gas.  It  is  commonly  expressed  (in  ar  electropositive  gas)  as 

dc 

M  =  exp  J  a(x)dx  (p) 

o 

where  a(x)  is  Townsend's  ionization  coefficient.  The  number  ot  elec¬ 
trons  released  at  the  cathode  per  incident  particle  (y)  is  influenced 
by  the  size  and  erergy  of  the  ion  or  photon  as  well  as  the  cathode 
material  and  the  state  ot  its  surface.  Thus,  the  cathode  fall  region 
is  where  processes  essentia1  to  the  maintenance  of  the  rest  of  the 
discharge  occur.  The  next  section  describes  the  sealing  relai  trships 


for  clow  discharges. 


2 .  SIMILARITY  RULES 


Although  the  various  parameters  describing  a  olow  discharge  are 
all  interrelated  arid  often  depend  on  each  other  through  complex  rela¬ 
tionships,  the  parameters  describing  two  different  discharges  can  be 
related  to  each  other  using  similarity  or  scaling  rules.  Two  dis¬ 
charges  are  said  to  be  similar  if  they  are  i cl « •  i . ;  ical  iri  one  or  two 
parameters  ana  different  in  others.  For  example,  two  discharges  in  the 
same  gas,  at  the  same  gas  temperature,  with  the  same  electrode  materi¬ 
al,  arid  with  the  same  potertial  difference  between  the  electrodes  are 
similar  even  if  all  the  corresponding  linear  dimensions  differ  by  a 
constant  factor.  This  includes  the  vessel  dimensions,  the  electrode 
dimensions,  arc  the  properties  of  the  gas  such  as  the  mean  free  path. 
Francis  (Reference  6)  has  comprehensively  derived  a  set  of  similarity 
rules  which  result  in  several  invariant  quantities.  Jf  the  discharge 
voltage,  current,  and  gas  temperature  are  the  same  for  both  discharges, 
he  has  shown  that  the  following  Quantities  are  also  irvariant: 


E/p  - 

electric  f ield/pressure 

pd  - 

pressure  *  distance 

J/P?  - 

current  oe nsity/pressure 

These  parameters  remain  invariant  for  respective  points  in  the  dis¬ 
charges  as  long  as  the  following  processes  dominate  the  kinetics. 

a)  Single  stage  processes  cuch  as  electron  impact  ionization, 
attachment,  or  detachment. 

b)  Renning  ionization 

r)  Charge  transfer 

d )  Drift  and  diffusion  from  a  volume 
if  two  stage  processes,  photoioriization .  or  recombination  are  the 


dominant  processes  ir  the  discharge,  then  the  similarity  rules  will  net 
be  followed  ana  the  sea'll  no  parameters  above  nay  no  longer  be  appro¬ 
priate.  The  above  invariant  quantities  apply  both  to  the  whole 
discharge  as  well  as  to  specific  regions  of  the  discharge.  Tor  exam¬ 
ple,  given  two  discharges  with  the  same  discharge  voltage,  current,  and 
gas  temperature,  the  invariant  quantities  for  the  cathode  fail  region 
would  be: 

Vc  -  cathode  fall  voltage  drop 

F/p  -  ratio  of  electric  field  to  pressure  at  respective 
positions  in  the  cathode  fall 
pd^  -  pressure  *  cathode  fall  length 
J/p 2  -  current  densi ty/(pressure ) * 
basically,  given  V,  I,  arid  T,  then  V^.,  F/p,  pd^,  and  J/p2  are 
constants  and  relate  how  the  electric  field,  cathode  fall  length  and 
current  density  scale  with  pressure.  The  next  section  describes  the 
experimental  effects  observed  and  attributed  to  negative  ions  in  low 
pressure  glow  discharges  including  the  cathode  fall  region. 

3.  EFFECTS  ATTRIBUTED  TO  NEGATIVE  IONS  IN  CAS  DISCHARGES 

Several  authors  have  reviewed  the  effects  of  adding  small  amounts 
of  negative  ions  to  a  glow  discharge  (References  5,  6,  8),  but  all 
iicjrtt  with  the  initial  description  advanced  by  Emeleus  and  Sayers  in 
1?.'X  (Reference  42).  Emeleus,  et  al.  observed  the  effects  on  the  dis¬ 
charge  when  a  trace  of  an  attaching  gas  such  as  chlorine  was  added  to  a 
neon  discharge  containing  a  little  helium.  li’Ji  the  adding  of  a  trace 
of  chlorine,  the  fluorescence  from  f.ht  pure  neon  ’scharge  i>  ..strated 
in  Figure  2a  changed  to  that  illustrated  in  figure  2b.  The  following 


a.  All  the  negative  zones  (cathode  dark  space,  4;  cathode  glow, 
3;  negative  glow,  5;  and  Faraday  dark  space,  6)  contracted  towards  the 
cathode  such  that  Aston's  dark  space,  or  the  primary  dark  space,  2, 
became  so  small  as  to  be  unobservable. 

b.  Curved  striations  appeared  in  the  positive  column,  8,  often 
accompanied  by  an  increase  in  the  potential  gradient. 

c.  A  new  dark  sheath,  7,  appeared  between  the  Faraday  dark 
space,  6,  which  was  moderately  luminous  and  the  positive  column,  8. 

d.  A  new  dark  sheath,  9,  also  appeared  between  the  positive 
column,  8,  and  anode  glow,  10. 

1? 


Spencer-Smith  (Reference  81)  and  more  recently,  L-avis  and  ['in; 
(Refercrxe  34)  observed  similar  hut  more  pronounced  effects  in  pure  io- 
oine  discharges,  and  in  helium  contairing  a  trace  of  iodine.  The  posi¬ 
tive  column  was  observea  to  constrict  radially  to  either  a  stationary 
form  or  a  mobile  ribbon.  This  type  of  constricu-c:  positive  column  has 
often  been  found  in  other  discharges  contain  ketrunegative  gases. 

fol lowing  the  discussion  by  Emeleus  and  Sayers  (Refererce  39),  the 
origin  of  the  characteristic  effects  of  elect  renegative  gas:es  was 
attributed  to  the  decrease  in  the  average  mobility  of  the  negative 
charge  carriers  as  electrons  are  converted  to  heavier  negative  ions. 
Three  almost  independent  effects  ensue: 

a.  Creation  or  addition  ui  negative  iors  corresponds  to  an 

increase  in  pressure  in  an  electropositive  discharge.  The  reduction  of 
the  drift  velocity  cf  the  electrons  by  their  attachment  to  gas  mole¬ 
cules  is  equivalent  to  the  more  continuous  retardation  caused  by  more 
frequent  collisions  from  an  increase  in  gas  pressure.  The  contraction 
or  the  negative  zones,  the  occurrence  of  str  .  in  the  positive 

cclumn,  and  the  constriction  of  the  positive  column  are  all  charac¬ 
teristics  of  an  ek-ctroposi  ti  ve  discharge  at  higher  pressures,  for 
example  100  torr  compared  to  1  frrr. 

b.  Since  negative  ins  are  more  massive  than  electrons,  they 
diffuse  nr  re  slowly  from  regions  r  +  space  charge1.  Conseon.  ntly , 
negative  iors  are  thoucl  t  to  set  u|  regions  of  resultant  negative  spare 
charge  and  to  make  existing  reokns  thinner.  This.  U:adr  k  ;  tie  forma 
tior.  cf  double  space  charge  layer  ^  cf  positive  iris  next  to  negative 
i'-ns  aria  the  conversion  cf  a  urururu  positive  column  to  a  sc  -a ted 
column.  The  high  field  gradient  between  these  layers  enhaiies  the 


formation  of  stationary  stnations  and  thus  the  enhanced  boundaries 
between  dark  anu  luminous  regions.  The  axial  field  resulting  from 
these  stationary  double  (+-;  space  charcr  regions  accelerates  electrons 
to  sufficient  energy  to  procuce  maximum  excitation  cf  the  gas  atoms  and 
molecules  leading  to  increased  light  emission,  lhe  dark  sheaths,  /, 
between  the  positive  column,  8,  and  the  Fare,  •  dark  space,  6,  and  the 
anode  glow,  9,  are  believed  to  he  the  result  e  +  space  charge  due  to 
negative  ions  created  in  previous  regions  of  the  tube  (References  5,  6, 
42). 

c.  Since  some  of  the  elections  become  attached  to  gas  molecules, 
the  radial  Gistribution  of  negative  charge  carriers  will  charge.  Due 
to  the  lower  negative  ion  mobility,  the  negative  ions  will  be  less 
likely  to  diffuse  away  from  the  slightly  positive  axis  it  the  positive 
column.  Since  detachment  of  electrons  from  negative  ions  can  often  be 
an  inefficient  process  in  many  gases,  large  concentrations  of  negative 
ions  can  build  up  near  the  axis.  This  causes  positive  ions  also  to 
drift  to  the  axis  to  balance  the  total  charge  resulting  in  a  con¬ 
stricted  Gischarge. 

The  above  description  is  based  mainly  on  qualitative  observations 
■..Inch  have  been  reviewed  in  Reference  5.  Several  quantitative  measure¬ 
ments  have  been  made  most  of  which  have  been  accomplished  since  1938, 
confirming  the  above  hypotheses.  These  will  be  reviewed  briefly. 

Spencer- Smith  (Inference  81)  (.rd  Woolsey  (Kt tereu*  <  Vfc )  used  probe 
measurements  and  magnetic  methods  for  measuring  charge  to  nx-ss  ratios 
to  investigate  the  negative  glow  and  positive  column  n-qions  of  low 
pressure  iodine  disc  he  rye.  Probe  measurements  in  .he  presence  of  large 
guai.tit.es  (1  negative  ions  are  vetg  difficult  to  interpret,  but 
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Spencer- Snn th  found  the  positive  and  negative  ion  concentration  m  the 
negative  glow  to  be  about  equal,  and  about  100  t. i roes  larger  than  the 
electron  concentration.  In  the  positive  column,  the  ratio  of  positive 
ions  to  electrons  remained  about  the  same,  but  the  absolute  concentra¬ 
tions  were  lower.  Spencer-Smith  also  found  the  electron  temperature  in 
the  negative  glow  to  be  normal  (~1  to  1.5  eV  -or  F/N  —  6?0Td.  The 
electron  temperature  in  the  positive  column  was  much  greater,  ~13  eV. 
This  car.  be  accounted  for  by  some  of  the  elections  being  accelerated 
through  several  potential  drops  as  a  result  of  striations.  However, 
the  conventional  concept  of  electron  temperature  under  such  ncnequi- 
libriur.  ci rcumstances  probably  is  not  valid.  Woolsey's  results  were 
similar  to  Spencer-Smith's  in  that  he  found  the  ratio  of  negative  ions 
to  electrons  to  be  2C0  in  both  the  negative  glow  and  the  positive 
column.  Woolsey  also  determined  the  radial  distribution  of  negative 
charge  carriers  to  be  much  flatter  than  the  normal  zero  order  Bessel 
function  found  in  the  negative  clow.  The  number  of  negative  ions  was 
found  to  be  much  larger  and  more  concentrated  around  the  axis  in  the 
positive  column  than  the  negative  glow,  in  agreement  with  Emeleus  and 
Sayer's  speculation.  In  both  experiments  the  positive  column 
fluorescence  about  the  center  of  the  discharge  tube  was  constricted. 

In  a  similar  experiment  using  probes,  Zimmermann  (Reference  98)  found 
the  positive  ion  concent  rations  more  than  fifty  times  greatei  than  the 
electron  concentration  in  the  constricted  positive  co’umn  of  a  hot 
cathode  discharge  in  hydrogen  chloride  vapor. 

Lurit  and  Gregg  (Referer.ee  63),  and  Thompson  (Reverence  84),  have 
investigated  oxygen  discharges.  Acain  using  probe,  and  r.a.ss  dialysis, 
Lunt  and  f.regq  found  the  concentration  of  0  in  the  positive  column 
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with  an  E/N  c:250  Td  to  be  about  equal  to  the  concentration  of 
electrons.  In  the  Faraday  dark  space  and  negative  glow,  however,  hard¬ 
ly  any  negative  ions  were  detected  (<  1 y ) .  Thompson  used  a  radio 
frequency  mass  spectrometer  probe  to  measure  the  particle  concentra¬ 
tions.  He  found  the  concentration  cf  0  was  about  equal  to  the 
concentration  of  0o+,  thus  confirming  Lunt  ai  'reg's  observation  that 
H~  was  the  most  abundant  negative  ion.  Thompson  also  measured  negative 
particle  energy  distributions  in  various  regions  of  the  discharge  and 
compared  them  to  the  electron  energy  distributions  in  the  same  regions 
in  a  pure  nitrogen  discharge.  The  energy  distributions  in  oxygen 
contained  an  apparent  large  low  energy  peak  due  to  negative  ions,  as 
expected.  He  determined  the  radial  space  potential  in  the  positive 
column  to  be  much  flatter  for  oxygen  than  tor  nitrogen  and  the  axial 
space  potential  changed  much  more  abruptly  through  a  striatior  than  in 
nitrogen.  Although  both  discharges  contained  standing  striations,  the 
striations  in  the  oxygen  discharge  were  much  closer  together.  The 
radial  distribution  of  the  electron  concentre  n  ,n  the  oxygen  dis¬ 
charge  was  almost  constant  whereas  the  radial  variation  of  positive  and 
negative  ions  both  decreased  almost  to  zero  supporting  the  idea  that 
the  positive  and  negative  ions  would  concentrate  near  the  axis. 

More  recently  Davis  and  king  (Reference  4)  observed  the  striations 
formed  in  a  He-I^  glow  discharge  laser.  They  calculated  the  ratio  of 
negative  ions  to  electrons  to  be  on  the  order  of  20-110  for  a  current 
density  of  100  0mp/m2  and  on  the  order  oi  2-11  for  a  current  densify  of 

4 

K  amp/m2.  1  he  negative  ions  with  their  lower  mobil'!  +  '  ;adily  led  to 
the  formation  of  both  stable  and  unstable  stationary  ana  moving 
striations  as  well  as  the  constriction  of  the  positive  column  when 


small  quartities  (<.  i  torr)  of  iodine  was  admitted  to  pure  hoi  i  ur; 
cischaryes  (7-10  torr).  The  forward  wave  (anode  to  cathode)  st  nations 
occurred  at  a  lUkH?  trequency  lar<.  results  ir.  a  partial  modulation  of 
the  laser  intensity.  The  radius  of  otic  constrictec.  positive  column 
decreased  with  an  increase  in  iodine  pressure.  Even  when  the  column 
remained  stable,  it  was  not  always  axially  synetric,  but  often  adopted 
a  twisted  conf iguration  in  the  discharge  tube  caused  by  fluctuations  in 
the  local  wall  temperature  and  iodine  density.  Thus  constricted 
discharges  are  undesirable  for  laser  applications,  since  they  prevent 
efficient  usage  of  the  available  resonator  volume.  The  instabilities 
which  result  in  these  striations  and  constriction  of  the  positive 
column  has  been  investigated  extensively  by  Nighan  and  Wiegand 
(Reference  70).  This  topic  will  not  be  investigated  in  this  study. 

The  occurrence  of  negative  ions  in  the  negative  glow  ana  positive 
column  is,  however,  relevant  to  the  cathode  fall  region  and  the 
boundary  conditions  relating  these  regions. 


SECTION  III.  ELECTRON  KINETICS  IN  THE  CATHODE 


FALL  REGION:  EQUILIBRIUM  ANALYSIS 

1.  REVIEW  OF  EQUILIBRIUM  THEORIES 

There  are  several  theories  for  the  cathode  fall  region  in  an 
electropositive  gas  which  relate  the  cathode  fall  voltage  (V(),  the 
discharge  current  (J),  and  the  cathode  fall  distance  (d  ).  These 
theories  will  be  referenced  later  as  they  are  discussed.  In  compari¬ 
son,  there  is  a  distinct  lack  of  any  kind  of  analytical  analysis  for 
electronegative  discharges.  These  theories  for  the  cathode  fall  region 
can  be  divided  into  equilibrium  analyses  based  on  Poisson's  equation 
and  hydrodynamic  current  continuity  equations,  and  nonequilibrium  anal¬ 
yses  based  on  calculating  an  electron  distribution  function.  The  equi¬ 
librium  analyses  will  be  analyzed  and  expanded  upon  in  this  section. 

The  nonequi 1 ibrium  analyses  will  be  considered  in  the  next  section. 

The  equilibrium  analyses  are  of  two  types;  those  which  regard  the 
cathode  dark  space  as  an  isolated  independent  unit  considered  separate 
from  the  negative  glow  and  those  which  include  the  negative  glow  so 
that  together  their  processes  sustain  the  discharge.  These  theories  of 
glow  discharges  are  thoroughly  summarized,  but  not  necessarily  criti¬ 
cally  reviewed  in  several  references  (References  6,  7,  89).  The  cath- 
ude  fall  current-voltage  relationship  derived  by  von  Engel  and 
Steenbeck  (Reference  90)  for  a  self-sustaining  glow  discharge  has  pro¬ 
vided  a  benchmark  for  all  following  papers  ard  will  be  briefly  reviewed 
first.  The  emphasis  will  then  be  on  results  presented  since  1v/5. 

The  von  Fngel  and  Steeribeck  theory  assumed  the  electric  field  in 


the  cathode  tall  tc  be  a  linearly  decreasing  function,  uf  the  font 


E=E  (1  -  (3) 

o  d 

c 

where  E  is  the  electric  field  at  the  cathode  ano  d  ib  the  cathode 
c  c 

tall  distance.  The  icr  flow  to  the  cathode  is  assumed  to  be  mobility 
limited  so  that 

\  =  ii+E  (4) 

where  v+  is  the  positive  ion  velocity  ar>c  u  their  mobility. 

Von  Engel  and  Steenbeck  used  boundary  conditions  that  set  the  positive 
ion  current  to  zero  at  the  cathode  fal 1 -negative  glow  boundary  and 
assumed  each  positive  ior,  incident  on  the  cathode  produced  y  electrons, 
yielding 

J+(d)=0  (5) 

and 

J  (0)  =  y J+( 0)  (6) 

Equation  (5)  is  not  entirely  correct  in  that  the  positive  ion  current 
does  decrease  at  the  cathode  fall-negative  glow  bourdary  but  dees  not 
necessarily  go  to  zero.  The  question  of  whether  the  number  of  positive 
ions  entering  the  cathode  fal1  from  the  negative  glow  is  significant  or 
not  has  developed  into  an  area  of  controversy  and  u  described  later. 


J  9 


Using  the  empirical  lorm  of  Townsend's  ionization  coefficient, 


p  A  exp  (- 


(/) 


they  derived  equations  relating  the  total  current  J,  the  cathode  fall 
voltage  V^,  and  the  cathoae  fall  width  dc>  The  cathode  fall  thickness 
is  derived  by  applying  the  maintenance  condition.  Equation  1,  to  the 
cathoae  fall  region  resulting  in 

i^c  ] 

T  dc  :  /  u(x)dx  =  ?.n(  1  +  -)  ^ 


Definition  of  the  width  of  the  cathode  fail  is  not  exact  but  it  is 
usually  defined  as  the  distance  from  the  cathoae  to  the  point  in  the 
discharge  where  the  electric  field  is  extrapolated  to  zero.  The 
voltage  is  determined  from  the  electric  field. 

d  Ed 

Vc  =  /  E(x)dx  -|--C  (9) 

o 


if  E(x)  is  given  by  Equation  3.  Since  the  ion  current  to  the  cathode 
is  generally  mobility  limited,  they  aerivea  a  corrent-vol tage  relation¬ 
ship  for  the  cathode  fall  from  the  high  pressure  space-charge-1 imited 
current  voltage  relationship  yielding 


JV  =  (1 


V+p 


(ID) 


Equations  8,  9,  and  1C  uniquely  describe  the  cathoae  lull  region 
and  ohry  similarity  or  scaling  rules  described  ir.  Eection  I. 

a 


Solving  these  equations  numerical ly  for  a  linearly  decreasing  electric 
field,  provides  good  agreement,  with  experiment  over  a  wide  range  of 
conditions  in  electropositive  gases.  Helium  is  an  exception  and  one 
reason  may  be  that  the  Townsend  ionization  coefficient  cannot  be  fitted 
with  an  exponential  fcrri  like  Equation  7.  However,  the  above 
assumptions  have  recently  become  controversial  areas  (Relerence  7). 

Druyvesteyn  and  Perming  (Reference  38)  pointed  out  that  if  a  sig¬ 
nificant  number  of  posit, ve  iers  entered  the  cathode  fall  region  from 
the  negative  glow,  the  maintenance  condition  for  the  discharge  was 
altered  to  the  following: 

1  +  1 

"dc  =  fnr+i  (id 

where  is  the  ratio  of  the  positive  ion  current  to  the  electron 
current  at  the  cathode  fall -negative  clow  interface.  Thus  the  value  of 
Tdc  could  be  reduced  if  the  positive  ion  current  was  comparable  to  the 
electron  current  at  the  cathode  f al 1 -negative  glow  boundary. 

Little  and  Von  Engel  (Reference  59)  refute  the  idea  that  a  signif¬ 
icant  number  of  ions  enter  the  cathode  fall  region  from  the  negative 
glow.  A  linear  electric  field  requires  a  constant  net  space  charge 
density  in  the  cathode  fall,  whereas  if  a  significant  number  o* 
positive  ions  originated  in  the  negative  glow,  then  the  positive  ion 
density  would  be  expected  to  decrease  towards  the  cathode.  As  Ingold 
(Reference  7)  pointed  out,  these  are  not  conclusive  arguments,  since 
such  field  measurements  are  very  difficult  and  any  s’owly  varying  field 
would  appear  approximately  linear  to  the  experimentalist.  Reviews  by 
Francis  (Reference  6)  and  Ueston  (Reference  93)  survey  many  attempts  to 


measure  the  electric  field  in  the  cathode  fell  region.  The  weight  of 
the  evidence  suggests  an  electric  field  decreasing  linearly  through  the 
cathode  dark  space  becoming  almost  zero  in  the  negative  glow.  Several 
recent  theoretical  papers  have  tried  to  show  the  eltviric  field  to  be 
linear  in  the  cathode  tall.  These  papers  will  be  reviewed  later,  in 
chronological  order. 

The  other  area  of  controversy  is  the  concept  of  Townsend's  ioni¬ 
zation  coefficient  (  0  and  the  assumption  that  the  electrons  are  in 
equilibrium  with  the  field.  t  is  the  number  of  ionizing  collisions 
per  centimeter  of  path  length  and  is  usually  determined  from  dr’ft  tube 
experiments.  These  experiments  assume  the  electrons  are  in  equilibrium 
with  a  uniform  electric  field.  This  condition  probably  is  not  not  in 
the  cathode  fall  where  the  electric  field  changes  rapidly.  Evidence  of 
this  appears  in  the  nonequilibrium  analyses  of  the  cathode  fa11  by  Lony 
(Reference  60)  and  Tran  Ngoc  An  (Reference  06)  which  are  discussed  in 
Section  IV.  Thus  the  theoretical  models  based  on  the  concept  that  the 
electron  distribution  function  is  ir  eoui 1 ibrium  with  the  field  may  not 
be  valid.  As  von  tngel  noted  about  his  own  theory,  the  equilibrium 
theories  give  "physically  the  correct  picture,  but.  the  numerical 
agreement  with  observation:  is  largely  fortuitous"  (Reference  89:224). 
This  observation  is  also  discussed  later  ir'  Section  V  where  the 
ect"i1ibrium  and  neneaui  1  ibrium  results  are  compared. 

Returning  to  the  review  of  the  literature,  since  1976,  11  'yariov 
fetererne  IT)  analytically  investigated  the  effect  of  the  positive  ion 
i  entriluit  ion  from  the  creative  gHw.  he  based  his  approach  or.  experi- 
nitcil  evidence  of  Gunther  shul  tze  (Reference  °6)  indicating  that  as  the 
,  rode  approaches  the  plasma  side  ei  the  cathode  sheath,  the  potential 


some  of 


drop  across  the  cathode  la’ll  increases.  This  indicates  that 
the  positive  ions  required  to  carry  -Jie  current  the  cathode  tall 
region  originate  in  the  negative  glow.  Ill 'yanov  solved  the  following 
set  of  equations  for  both  high  and  low  ior  mobility  cases,  defined 
bel  ow 

dE  _  en+ 


dx  10 

(12) 

dJ 

dT  =  (J“J+,u(E) 

(13) 

cx(E)  =  Ap  e  E 

(14) 

=  en+u+E  (low  field  case) 

(15) 

en+k+E^2  (high  field  case) 

(16) 

The  boundary  conditions  which  neglect  any  potential  drop  in  the  nega¬ 
tive  glow  are  given  by 

E(d)  =  0  07) 

J,(d)  =  (18) 

where  was  a  variable  parameter  representing  in  his  theory  the  ratio 
of  positive  ion  current  density  to  electron  curror  density  at  the 
cathode  1  a  1 1 -nega t i ve  ylow  boundary.  For  both  mobility  cases,  his 


final  transcendental  equations  contained  integrals  winch  had  to  be 
solved  numerically.  His  results  indicated  that  if  the  fraction  of  ion 
current  coming  from  the  negative  glow  (A)  decreases,  then  the  electric 
field  at  the  cathoae  increases.  For  a  given  secondary  emission  coeffi¬ 
cient  (a)  ,  he  determined  that  the  normal  cathode  fail  potential  (ly), 
normal  current  density  (J^),  and  normal  cathode  fall  length  (x^,) 
changed  quite  slowly  over  a  wide  range  of  <5  values.  He  found  the 
greatest  discrepancy  between  his  theory  and  von  Engle's  theory  at  small 
where  a  singularity  existed  in  his  equations  for  \y  and  J^.  His 
equations  for  V..,  J^,  and  depended  only  on  gas,  cathode  parameters, 
and  'v>  ,  once  the  transcendental  equation  for  the  electric  field  was 
solved.  From  the  fact  that  \y,  d^,  and  x^  were  a  function  of  5  ,  he 
concluded  that  a  complete  theory  for  the  cathode  fall  normal  glow  dis¬ 
charge  must  take  into  account  ion  production  in  the  negative  glow. 

Neuringer  (Reference  69)  made  an  analytical  investigation  of  the 
cathode  fall  region  in  high-voltage  low-current  discharges.  He 
determined  that  to  first  order  the  electric  field  decreased  linearly  in 
the  cathode  fa!1  region,  the  voltage  drop  depended  on  the  two-thirds 
power  of  J/p2  and  the  cathode  fall  thickness  was  independent  of  the 
operating  current  and  varied  inversely  with  pressure.  His  model 
combinec  the  electron  current  continuity  equation. 


Poisson's  equation, 


dF  -!  rjF(x)  Je(x> 

dx  "  Fo  [rf+(E(x))  "  v'TCRIT 

witli  the  requirement  for  total  current  conservation. 


(? 


v*J  =  0  or  ,1  =  J  +  J  =  constant  (2 

+  e 

and  J+  are  the  electron  anci  positive  ion  current  oensities  and  vg 
and  v+  are  the  electron  and  positive  ion  drift  velocities.  He  used 
Townsend's  ionization  coefficient  (Equation  7),  Ward’s  (Reference 
92)  modified  formula  for  the  positive  ion  velocity  for  ions  in  a  hi  oh 
field. 


v 


+ 


<+(  /P) 


1/2 


(2; 


and  the  fact  that  v+/v  «  1  throughout  the  cathode  fall  to  obtain  a 
differential  eouation  separable  in  J  and  E.  This  equation 


(?■ 


can  be  integrated  attuning  the  electric  field  goes  tc  vie  and  the 
electron  current  density  is  equal  to  the  total  current  density  at  the 
cathode  fa1 1 -negati ve  glow  boundary. 


he  was  able  to  derive  an  expression  for  the  electric  tield  as  a  func- 


drop  terns.  Additional  analysis  in  the  next  section  includes  negative 
ions  and  reveals  that  the  lirst  order  approximations  do  not  have  to  be 
made.  Nevertheless,  his  results  add  credence  to  von  Fred's  assumption 
that  the  electric  field  decreases  linearly  in  the  cathode  fall  region. 

Davies  and  Evans  (Deference  ?■'.')  extended  Neurinoer’s  analysis  and 
showed  it  could  be-  applied  to  the  whole  discharge  by  changing  the 
boundary  conditions  and  assuming  the  ratio  of  positive  icr  to  electron 
mobility  to  be  a  constant. 

d 

r  =  vd  ( 29) 

e 

This  assuiTiption  was  found  to  breakdown  only  at  electric  fit-ids  for 
which  1  -  J  «  rJ^,  which  normally  was  ret  the  case.  The  boundary 
condition  at  the  cathode  was  the  same  as  Neurincer's  i.e.  Equation 
? 7,  hut  the  other  boundary  condition  was  changed  to  represent  the 
electric  field  at  the  anode  or  to  specify  a  particular  applied  voltage. 


E(d)  --  E 

a 


(30) 


d 

V  =  f  E(x)dx  ( 

o 

based  upon  combining  Frisson's  Equation  2L  and  the  electron  current 
continuity  Equation  19  into  a  transcendental  type  or  equal inr  . 
including  the  above  changes,  they  developed  a  new,  fast  numerical 
technique  for  predicting  the  variation  of  E  and  j  as  a  function  of 
distance  from  cathode  to  anode.  Although  Equatior  i9  is  true  in 
the  cathode  fall  region  where  ionization  is  the  dominant  procers  and 


27 


the  electrons  are  scattered  predominantly  in  the  forwarG  cl  i  recti  or:,  it 
is  not  correct  in  the  positive  column  re pi  on  wheai  the  election'  are 
scattered  isotropically  and  the  gain  ol  electrons  is  balanced  by  losses 
to  radial  diffusion,  recomb  inatiori ,  or  attachment,  ir,  essence,  when  r 
' s  hie  1  d  constant,  this  is  tantamount  to  includirc,  diffusion  n  the 
positive  column.  However,  diffusion  was  not  'included  as  a  loss  process 
in  the  electron  current  continuity  equation.  Ir,  addition  to  this 
aiscrepancy,  the  anode  f a  1 1  voltage  appears  to  be  much  too  small. 
Estimating  from  their  Figure  1,  the  anode  fall  in  helium  appears  to  be 
less  than  b  volts  which  is  much  less  thar  the  nominal  value  of  26  volts 
(Reference  6:14f).  Later  in  this  chapter  the  same  problem  was  encoun¬ 
tered  when  this  method  was  useu.  Further  discussion  of  the  anode-  fall 
is  reserved  for  Section  V. 

Recently,  Mitchell.  Kline,  ano  Denes  (Reference  67)  calculated  the 
electric  field,  particle  fluxes,  ana  particle  densities  as  a  function 
of  position  in  the  cathode  and  anoae  region'  rf  a  high  pressure  glow 
discharge  containing  an  attaching  gas  mixture.  They  used  a  one  di men¬ 
s’  cnal  steady  state  continuity  equation  model  including  Poisson's 
eourtion.  The  following  boundary  conditions  were  used: 

Cat !■' de  Anode 

•  n  r.  n  0 

e+  e  +  - 

J4  =  0 

where  J  /J  and  E  were  adjustable  initial  paramett.s.  They  rf,  cited 


Positive  Column 

dJ  dn  dE 

dv  dx  dx  (J2) 


that  the  cathode  sheath  characteri  sties  are  insensitive  to  vari.it  ions 


in  the  concentrations  of  the  attachino  species,  whether  diffusion  or 
electrons  were  included  or  not.  However,  they  found  the  anode  sheath 
wa s  affected  by  the  inclusion  of  electron  diffusion  as  ve1!  as  the 
attaching  concentration.  As  will  be  shewn  in  the  next  section,  the 
boundary  condition  for  the  positive  column  in  Equation  (32)  may  not 
lead  to  a  unique  solution.  Thus  there  remains  some  question  on  t he 
convergence  arid  lienee  self-consistency  of  their  solutions. 

2.  NEW  EQUILIBRIUM  ANALYSIS  OF  THE  CATHODE  REGION  INCLUDING  NEGATIVE 
IONS 

This  section  expands  Neuringer's  analysis  of  the  cathode  fall 
region  to  include  negative  ions.  Since  the  ratio  of  positive  ions  to 
electron  mobility  is  held  constant,  the  transverse  diffusion  ct  elec¬ 
trons  was  added  to  the  electron  continuity  equation.  A  one-dimensional 
steady  state  distribution  is  assumed  so  all  quantities  are  functions  of 
the  distance  x  from  the  cathode  surface  and  not  of  the  transverse 
geometry.  When  attachment  and  transverse  dil fusion  processes  are 
included,  the  equations  governing  the  steady-state  discharge  are: 

Pi  i  sson  '  s  equation 


the  electron  current  continuity  equation 


the  negative  ion  current  mntinuity  equation 

i.  .1 

dJ  ,  _  r  -  * 

dx  J r  PV.V, 

i.o  the  total  current  conservation  emotion 

V -J  =  0  or  J  -  +  J  *  J_ 


^  ah  J 


136) 


where  L  is  the  electric  field  strength,  J+,  J  ,  J  and  J  are  the 
positive  ion,  electron,  negative  i or  and  total  current  densities 
respectively,  v+,  v  ,  and  v  are  the  positive  icr,  electron,  and 
negative  ion  drift  velocities,  r  is  the  permittivity  of  free  space,  a 
and  n  are  Townsend's  first  ionization  and  attachment  coefficients,  6 

e 

is  the  electron  free  diffusion  coefficient  (cu^/sec),  a  is  the  electron 
mean  free  path,  k  is  the  ion-ion  recombination  rate  (crn3/sec ) ,  and  x 
is  the  distance  from  the  cathode  to  the  anode.  These  equations  differ 
from  Davies  and  Evan's  equations  in  that  negative  ions  have  been  aooeci 
to  Poisson's  equation,  attachment  and  electron  diffusion  have  been 
added  iu  the  electron  continuity  equation,  the  negative  ion  current 
continuity  equation  is  now  required,  and  the  negative  ion  current 
density  has  been  added  to  the  total  current  conservation  equation. 

Normalizing  the  current  density  to  unit  implies: 

J  J 

j  =  _1  i  =  " 

e  J  T  (.37) 

Then  dividing  Equation  33  by  34  and  usiny  37  to  eliminate 


the  positive  ion  current  density  yields 


_ I  _ 

v  1/  D  \ 

. .0 


Equations  3'1 ,  35,  and  3E  u»e  tne  set  wl  ich  will  he  analyzed  instead  of 

the  original  system  oi  Equations  33,  34,  and  35.  Because  '*>  n,  v_.  , 

v  are  strong  functions  of  the  electric  field,  this  system  of  couplec 

equations  is  highly  nonlinear.  The  system  becomes  solvable  if  the 

to  1  lowing  assumptions  are  made: 

v  v 

^  =  r  and  -±  -  1  (39) 

e 

Ihe  variation  in  r  with  the  electric  field  is  very  small  and  becomes 
significant  only  when  J+  »  rJg.  Thus  these  are  good  assumptions  over 
the  complete  range  of  E.  Making  use  of  these  assumptions ,  separating 
variables  and  using  the  boundary  conditions  E/p  -  E  /p  and  j  =  j  at 
the  cathode  results  in  the  following  integral  equation 


t  /  P 

J  (“' 


%-  )v ,pd(E/p)  = 


-T  U 

x  v  / 


f  U-(Wr),ie-2j_]  j 


The  Townsend  first  ionization  coefficient  that  was  used,  was  the 


more  general  form  given  by 


-  =  A  e 
P 


(B(p/En 


where  A  and  B  are  empirically  determined  constants  characteri st.i c  of 
the  gas,  and  s  =  l,  for  molecular  gases  and  s=$  for  mu;, atomic  gases. 
Similarly,  the  attachment  coefficient  can  be  expressed  ac  an  exponen¬ 
tial  + 

,/P  =  C  e-(D(E/P>  > 


(4k) 


ur  as  d  series 


r'/P  =  ao  + 


^  (E/p) 


where  C,  L,  and  t  or  ai  are  errpi rical  !y  determined  constants  depending 
on  the  attaching  gas.  these  are  gocc  assumptions  provided  the  elec¬ 
trons  are  in  equilibrium  with  the  field.  The  ion  dritt  velocity  in  a 
high  f’eld  region  can  be  described  by  Ward's  modified  formula. 


v+  =  k+  (E/p) 


Substituting  for  a,  n,  anci  v+,  Equation  40  becomes: 
°/Tflp  e-B(p/E)s.Cp  e-D(E/p)^_  V  1,.  ,.,.0/2 
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/  p  . 


+(E/p)1/2_ 


k+( E/p)  '  pd( E/p) 


Jeo  r . 

,/  b 


-  -  (1+r )-2  dj 
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Changing  variables  by  letting  ip  =  B(p/E)s,  the  first  term  on  the  left  side 
of  the  equation  becomes  the  difference  of  the  two  incomplete  gamma  functions 


°f  Ap  e"B(p/e)Sk+(E/p)1/2pd(E/p)  =  Ap2k+B3/2s  f  e^dil 

E/p  ¥  ^  +1 


A  2k  BJ/  -  s 

-B-±- -  [r(-3/2s,B(p/E  )S)  -  r(-3/2s  ,  B  ( p/E  )  )] 

S  0 


:  a-1 
t  dt 


,  .  t  -t  a-1 

Where  r(a,x)_  let  dt  (Hererence  9:?6t)  Similarly 

chancing  variables  so  that  <t>  =  D(E/p)f  ,  the  second  ten,,  becomes: 


2 
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The  third  term  integrates  easily  to  yiele 

H  rp  D  rp 

_e_  d ( E/p)  =  [EQ/p 

r  /  X  ' 

E/P 


Integrating  the  right  hand  siGe  of  45,  yields 


eo 

■t  it 


L_  -  0+r)-2 


-Idj  =  p  n(4^  )+(l+r)(j  -j  )+2 
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Therefore  45  can  now  be  rewritten  as 


B  f  y  3  c  1  +  r  3  t1  AB“r 

~  [  (-  2?,  B(p/E)S)J  +  to  7t  |_’(2t,  D(E/p)r)J  '  >\p  E/p  s  ^ 

C _  ~  q  ^  j 

toM  1  lt,0(Eo/p)t)J  •  7^  (E0/p)  •  T7r[;n(i>  *  (Hr)(jeo' 


f  2s,B(p/Eo) 


.  (50) 
3eo  3- 

(l+r)( j  -j  )+2  f  j\  dj 
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The  gamma  functions  can  be  simplified  for  numerical  calculat  „  with 
the  help  of  the  following  recurrence  relation  (Refererce  15:457). 


r(a+l ,x)  =  ar(a,x)  +  xap  x 


For  s  =  £,r(o/2s,  x)  reduces  to 


-x  i  12 


'  ( -3  ,x )  =  -  g[E1 (x)  -  e‘  (  x  -  ■  2  +  -3)] 


where  E^(x)  /  ^  dt  ,  the  exponential  integral  (Reference  fJ:2?8). 

x 

For  numerical  calculations  t,(x)  oar  easily  It  ropresenteo  Fy  polynomi¬ 
al  app'oxTation  (Reference  9:251).  For  s  1,  "u/<-S,  x )  reduces  to 


4  .  2  -x 

^  >  erfc  ix  +  3  e 


(  '3/2  ’  1  /  ) 

'X  X  / 


where  erfc1/.) 


■’  r  _  t 

J  e  dt  is  the  complementary 


•y  error  ruu.!ion  for 


which  there  are  tables. 


Similarly  lor  t  -  i,  '  (+3/2t,  x)  reouces  tc 


-x 


—  erfc 


✓x 


(6«) 


These  simp"!  if  ications  for  the  qanria  functions  above  were  used  in  the 
numerical  calculations  to  ease  computation.  Error  functions,  exponen¬ 
tials,  and  square  roots  were  more  accessible  in  computer  libraries  than 
incomplete  gamma  functions. 

These  terms  still  form  a  trarcendental  equation  involving  j  and 
E.  Giver  i  at  the  cathode,  various  mathematical  methods  can  be  used 
to  find  E  which  is  a  root  to  the  equation.  This  field  can  then  be  used 
in  solving  the  continuity  equations,  3 e  ana  35,  very  near  the 
cathode.  Iteratively,  one  can  spatial 1>  step  from  the  cathode  to  the 
anode  by  alternatively  solving  for  E(je)  and  j  .  This  method  elimi¬ 
nates  the  integration  problems  conventional  differential  equation 
solvers  have  in  solving  the  original  set  of  equations,  which  included 
Poisson's  equation  explicitly.  E(j  )  is  a  smooth  furction  and  can  be 
determined  using  successive  numerical  approximations  to  solve  the 
transcer. dental  Equation  50.  See  Figure  3. 

An  expression  for  the  variation  of  the  electric  field  car  be 
obtained  by  first  taking  the  derivative  of  Equation  50.  The 
den v alive  ot  each  of  the  terms  in  50  is  given  below. 
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The  derivatives  involving  j  car:  be  eliminated  by  substituting  in 


Equation  34  in  normalized  form  yielding 


(E/p) 
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In  an  electropositive  gas,  j_  =  C  and  «  1  very  near  the  cathode  so 
the  last  two  terms  in  59  are  negligible.  Integrating  Equation 


59  under  these  conditions  yields 


(E/p)£-  =  (Ep/p) 


2  _  3 Jx 
‘  2k+PF, 


or  taking  the  first  two  terms  in  a  binomial  expansion  yields 


E/p  -  E  /p  -  — - - 

0  Me0(E0/p)1/2 

This  is  the  same  result  heuringer  obtained,  26,  using  many  more 
approximations  in  his  derivation,  lie  retained  cnly  the  first  term  in 
Equation  49  as  he  assumed  j  «  J  near  the  cathooe.  This  resulted 
in  c  r.  equation  similar  to  60  which  contained  terms  proportional  to 
( E/p ) ^  which  he  conveniently  dropped.  The  binomia.  expression  leading 
to  Equation  61  cr  iverges  as  long  as  (£  /p)  ^  - 


This  inequality  holds  in  the  cathode  fall  region  but  may  rot  hclG  in 


the  negative  glow  or  positive  column  region  where  E^/p  becomes  very 
small  and  no  longer  a  linearly  decreasing  function. 

Returning  to  Equation  59,  the  dominant  terms  for  an  electro¬ 
positive  gas  are 


j  (i-y  =  jj. _ 

k+p2eo  (E/p)1 ktp2co(E/p)1/? 


which  is  consistent  with  Poisson's  equation  and  the  fact  that  the  field 
in  the  cathode  fall  region  is  determined  primarily  by  the  positive  ion 
number  density. 

In  an  electronegative  gas,  the  situation  is  more  complex.  How¬ 
ever  very  near  the  cathode  even  in  an  electronegative  gas,  j  «  Je  and 
Je«  i.  Thus  the  field  there  is  determined  primarily  by  the  positive 
ion  current  in  both  electronegative  and  positive  gases.  As  one  pro- 
ceeos  across  the  cathode  fall,  the  electron  and  negative  ion  current 
dens-* t ips  increase  in  a  related  manner  given  by  35.  This  has  the 


affect  of  increasing  the  slcpe  of  the  electric  field  as  can  be  seen 
from  the  following  analysis. 

It  was  found  from  the  results  of  the  next  section  that  the  nega¬ 
tive  ion  current  density  through  the  cathcde  tall  region  could  be 
related  to  the  electron  current  density  by  an  exponential  function  of 
the  following  form, 


where  c  is  a  constant  determined  by  the  boundary  condition  that,  j  =  C 
the  cathcde  and  in  is  an  empirically  dctermineu  constant.  Results  ol 
the  nonequi  1  ibriun  analysis  i r.  the  next  Section  also  confirm  that 


37 


Eaeat.ion  63  is  c  good  appro iriation  throughout  the  cathode  tail 
region  txcept  for  a  very  smal1  region  a  few  electron  me  a  t  tree  paths  ' r 
length  in  front  of  the  cathode.  Using  this  approximation,  the  deriva¬ 
tive  of  the  integral  on  the  right  side  of  59  becomes 


(64) 


After  integrating  59  and  using  this  expression,  one  can  solve  for 
the  field  as  a  function  of  distance 


■S  J 

(E/p)  =  (E0/p)  ~  2T+pf.o  |x(l+2c)-(l+r)  ^  Jedx-2  J' 


Taking  the  first  two  terius  in  a  binomial  expansion  as  before  yields 


k  pc  (E  /  p) 

+  00 


1 -77  [x(l  +2c)-(1  +r )  /  Jed‘  '"I 


n'Je  d> 


In  the  electronegative  case,  the  field  is  shitted  higher  at  each 
corresponding  x  near  the  cathode  by  the  factor  (1  •>-  Lc).  This  shift  is 
small  but  can  be  observed  in  the  results  presented  in  the  next  section 
as  well  cis  the  results  of  the  neitequi  1  ibrium  analysis  i r.  Section  IV  for 
mixtures  of  Xe  and  HC1 .  Thus  this  analysis  predicts  that  the  attaching 
gas  slightly  lengthens  the  cathode  fall  arc  slightly  increases  the 
cathode  fall  voltage.  As  the  eiectron  n  1  rent  dens' grows  with 
respect  to  distance  the  integral  terms  become  larget  ..no  have  the 
effect  of  increasing  the  fie’d  in  the  positive  column  region  of  the 
discharge.  Thus  increasing  the  amount  of  an  attac+er  also  inc. eases 
the  field  m  the  positive  column  region  in  an  attempt  to  balance  the 
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qair  and  loss  of  electrons. 

Returning  to  Lquation  59  and  solving  for  the  minimum  slope  of 
the  electric  field  yields 


d ( E/p) 
dx 
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The  minimum  field  ir  an  electropositive  gas  occurs  when  j  =  1/1+r.  In 
the  electronegative  case,  the  minimum  field  occurs  when 
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e 


mj 

1 +2(c-e  e) 

(T+FJ 


l-2j 

1  +r 


(68) 


Both  of  these  expressions  are  equivalent  to  charge  neutrality. 


n+  =  ne  +  n 


(69) 


3.  DESCRIPTION  OF  PROG  RAN.  GLOW 

The  GLOW  program  was  oeveloped  during  this  study  to  provide  a 
fast,  yet  complete  description  uf  the  electric  field  in  a  gas  discharge 
from  cathode  to  anode,  similar  to  the  Davies  and  Evan's  model  { Refer¬ 
ence  32).  This  type  of  model  was  chosen  because  it  can  give  a  complete 
description  of  the  electric  field  as  a  function  cf  dislarce  Irom 
cathode  to  erode.  The  advantage  of  this  type  ot  technique  is  that  it 
minimizes  the  number  or  boundary  conditions.  This  code  differs  from 
most  other  approaches  including  Davies  arid  Evans  by  starting  a...  the 
cathcce  and  proceeding  towaros  the  anode.  Most  other  approaches  start 
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■it  the  anode  or  positive  col  urn  and  work  Lack  toward?  the  cathooc  coo 
to  stability  problems.  This  method  is  sufficiently  stable  to  start  at 
the  cathode  where  the  boundary  conditions  on  the  negative  tor,  current 
are  well  known.  The  theory  upon  which  the  model  is  based  has  already 
been  described  in  the  preceding  section.  This  section  will  describe 
the  boundary  conditions  and  the  iteration  tec!  i  (,i  e  s  in  detail. 

The  boundary  conditions  usee  at  the  cathode  and  anode  are  ones  on 
whirl  there  is  fair  agreement.  The  boundary  conditions  used  at  the 
cathode  are 


where  is  the  secondary  emission  coefficient  at  the  cathode.  The 
other  boundary  condition  requires  current  continuity  at  the  anede  such 
that 


Je  +  J_  =  J  or  J+  =  0 

Other  techniques  such  as  starting  in  the  positive  column  using  a 
differential  equation  solver  and  starting  toward  the  cathode  usua  Hy 
re  quirt  one  cr  more  boundary  conditions  in  addition  <  Equation  /() 
and  /l  'o  match  solutions  for  each  region  being  cescnbec. 

The  cede*  is  based  on  a  shooting  technique  which  varies  the  elec- 
mc  i  if  Id  at  the  Cathode  to  converge  upon  c  calculated  value  uf  the 
electric  field  at  the  anode.  The  value  of  the  eiertric  field  at  which 


the  electron  current  density  p1  us  the  negative  ion  current  density 
equuis  the  total  discharge  current  density  was  found  to  converge  very 
quickly  to  a  constant  value.  This  value  of  the  electric  field.,  at 
which  the  anooe  boundary  conditions  fold  was  the  field  that  v.os  the 
target  in  the  shooting  technique.  After  each  iteration  frnm  cathrac  to 
anode  the  target  value  of  the  electric  field  at  the  anode  was  updatco 
with  a  refined  value*  for  which  Equation  72  hole's.  If  F(d'  was  too 
high  such  that 

(73) 

J  +  j  ■>  J 
e 

then  E(l)  was  decreased.  If  E(o)  was  too  low.  nr  the  calculation  did 
not  reach  the  anode,  then  E(0)  was  increased.  Fiyurp  A  shews  the 
electric  field  as  the  program  converged.  For  the  second  ana  twelfth 
iterations  ir:  Figure  4,  no  solution  existed  for  the  transcendental 
equation  past  the  points  plotted.  This  implies  that  this  technique 
converges  not  only  to  a  solution  which  satisfies  fhe  boundary  equa¬ 
tions,  tut.  is  also  the  lowest  voltage  for  the  discharge  lor  the  initial 
conditions. 

Any  numerical  technique  for  finding  the  roots  of  an  equation  could 
be  used  in  solving  the  transcendental  equation,  such  as  the  method  of 
halving  the  interval.  Flowever,  o  modified  linear  interpolation  method 
(Reference  4b:lC'-ll)  was  used  since  u  s  supposed  tc  te  sligthly  fast¬ 
er  than  the  halving  the  interval  method.  The  only  rec  r  .  r  ernen  i  was  Itat 
the  first  twe  estimates  for  E(o)  should  bound  the  soiutior,  : . o .  one 
estimate  be  too  large  and  the  other-  ‘eo  small.  Rapid  convergence  to 
4-E  digit  accuracy  was  achieved  withir  TO  iteratic  .  from  ele<  nc 
fields  differing  as.  much  as  UJ-lf.v/crr.  at  the  cathode.  Numerical 
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calculations  were  terminated  when  the  target,  E(d),  was  within  1C;-  o4- 
the  calculated  value. 

Within  each  of  these  iterations  from  cathode  to  anode,  a  differen¬ 
tial  equation  solver,  and  this  modified  linear  interrelation  method  was 
used  to  solve  the  system  of  current  continuity  equations,  34  and 
40,  ard  find  the  roc;  to  the  transcendental  to nation,  50, 
respectively.  The  variable  order  Adams  predictor-corrector  method 
which  was  part  of  a  differential  equation  solver  package  called  DGEAR 
(Reference  36)  for  solving  stiff  differential  equations  was  used  to 
solve  tor  j  given  an  electric  field.  The  over?1’  iteration  scheme 
from  cathode  to  anode  went  as  follows; 

a.  Given  E(0),  DGEAR  was  called  to  solve  the  set  of  Equations 
34  and  35  for  jp(x)  a  short  distance  from  the  cathode. 

b.  Given  this  value  of  jt(x),  the  modified  linear  interpolation 
method  was  used  to  find  the  root  or  equilibrium  value  for  E(;<)  from  the 
transcendenta I  equation  50. 

c.  Given  this  value  of  E(x),  DGEAR  was  Celled  again  to  solve  for 
a  rew  j  (x)  at  the  next  small  increment  in  space.  The  spatial  step 
size  h  was  usually  picked  such  that  (E(0)  -  E(h))/E(0)  <  1G.  This 
cycle  of  calculating  j  (x)  and  E(x)  was  repeated  as  the-  program  stepped 
from  cathode  to  anode. 

The  program  can  be  run  without  the  help  of  the  differential 
equation  solver,  using  the  field  at  the  Iasi  step  to  calci lure  the 
field  from  Equation  50.  However,  this  method  was  very  sensitiie  to 
the  spatial  step  size  used.  The  addition  of  a  differential  equation 
solver  to  predict  the  electron  current  density  at  me  next  step  where 
the  transcendental  Equation  50  was  to  be  solved  next,  enabled  much 
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larger  steps  to  be  taker  with  a  savings  in  computer  time. 

Tlie  GLOW  program,  incorpcrati ng  both  a  differenti al  equation 
solver  and  the  modified  linear  interpolation  method  was  usually  run  cn 
a  Cyber-  750  computer  system.  Each  iteration  took  on  an  average  less 
than  4  sec  of  execution  time  for  250-300  spatial  steps.  It  required 
less  than  52k  ol  memory  to  load  and  run.  For  deleters  used  in  the 
numerical  calculations  are  described  in  Appendix  C. 

To  establish  a  baseline  and  verify  the  ccuracy  of  the  technique, 
Davies  ana  Evans'  boundary  conditions  and  ■>' n i t i a  1  values  were  used  and 
the  resultant  field  plotted  in  Figure  5.  They  set  the  electron 
current  density  equal  to  the  total  discharge  current  and  then  varied 
the  electric  field  at  the  anode  until  a  fixea  value  of  the  field  was 
obtained  at  the  cathode.  By  starting  at  the  anode,  Y  was  not  required 
as  a  boundary  condition.  The  GLOW  code  was  first  run  from  ar.cde  to 
cathode  to  determine  <  at  the  cathode.  The  code  was  then  run  from 
cathode  to  anode  using  this  value  of  <  =  .293.  In  this  instance,  1  was 
fixed  ana  the  field  at  the  cathoae  was  varied  until  the  electron 
current  density  equaled  the  total  discharge  current  density  at  the 
anode.  As  can  be  seen  from  Fiyure  5,  very  gooo  agreement  of  the 
GLOW  code  icathode  to  anode)  is  obtained  with  Davies  and  Evan's  result 
(anodt  to  cathode). 

4.  hLSUITS  OF  Nl'f, EPICAL  CALCULATIONS 

This  section  presents  the  results  of  calculations  usinc  the  GLOW 
code  ir  various  mixtures  of  rare  gases  with-  PCI.  All  data  presented 
for  Me  o r iC  Ar  was  obtained  for  a  gas  pressure  of  I  torr.  The  ;ata  tor 
Xe  was  calculated  for  a  gat  pressure  of  ICO  torr.  however,  it  was 

s4 


found  thut  for  al  i  three  rare-cas  mixtures.,  the  results  at  ditTtrert 
pressures  obey  ex.  the  following  pressure  scaling  relationships. 

V  "  constant 
J/p1’  =  constant 
pd  =  constant 

This,  finding  is  in  accordance  with  the  similarity  rules  or  scaling  laws 
ret  gas  discharges  discussed  in  Section  I.  Dates  was  taken  at  1.  10, 
and  iU(J  torr  for  two  gas  mixtures  and  was  found  to  deviate  less  than 
•  Olio  from  these  scaling  relationships. 

A  current-voltage  relationship,  similar  to  l.quation  10  can  be 
derived  for  the  high  i i e 1 d  case  where  the  positive  ion  velocity  is 
proportional  to  the  square  root  of  the  field.  Following  von  Fngel  and 
Steenbeck's  derivation,  but  using  Equation  44  for  the  ion  drift 
velocity  results  in  the  following  current-vo1 tage  relationship. 


j/p2  =  i  . 4344r  k  ( 1  +  Y)  -c  --0 
0  (pd)V< 


is  a  constant  for  a  given  Y  and  in  the  high  field  case.  This 
scaling  relationship  is  useful  ii  ccmparing  data  at  uiltercnt  voltage 
and  currents. 

Results  will  be  presented  in  the  secuence  Me-  matures,  Ar  nix- 
lures,  ana  then  Xu  mixtures.  Results  will  he  riistwsed  in  th.  >  section 
lor  the-  cathode  fall  e  rd  positive  column  regions,  but  not  tor  the  anode 
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fall  region.  At  present  the  CLOU  code  does  nr t  adequately  describe  the 
anoae  region.  The  anode  tali  region  is  discussed  ir.  Section  V. 

Helium  Mixtures 

This  section  presents  the  results  obtained  +  or  He/HCl  mixtures; 

given  y  -  1  torr,  d  --  3cm,  J  =  1.6F-5  amp/cm2,  and  Y  =  .2.  A  summary 

of  numerical  values  tor  the  following  figures  is  given  in  Tatle  1  for 

L  / p,  E.„ /p,  pd  ,  J  /p2,  and  V  .  E  is  the  electric  field  at  the 
o  minrrc  r  c  o 

cathode,  E^ ^  is  the  minimum  electric  field  which  in  these  calculations 
occurs  in  the  positive  column,  p  is  the  total  gas  pressure,  dc  is  the 
length  of  the  cathode  fall  region,  J  is  the  total  current  through  the 
discharge,  and  V  is  the  cathode  fall  voltage.  Consistent  with  the 
somewhat  arbitrary  but  common  scheme  for  calculating  the  cathode  tali 
distance,  the  linear  cathode  field  region  was  extrapolated  to  zero  to 
determine  the  cathode  fall  length  dr.  The  cathode  fall  voltage  was 
obtained  from  the  vo-tage  between  the  cathode  and  thic  point.  (Other 
definitions  can  he  used,  but  for  the  comparison  with  other  published 
models  this  delinition  was  applied.  Experimental  comparisons  would 
have  slightly  lower  values).  The  last  column  in  Table  1  is  the 
scaling  function  (SF)  given  in  Equation  75.  This  scaling  relation¬ 
ship  help!,  to  put  in  perspective  the  various  data  on  cathode  fall 
parameters.  The  value  of  SF(i,  K  /  for  100?  He  agrees  only  within 
0 2.5?  of  Ward's  theoretical  data  (Reference  S2).  This  discrepancy  is 
largely  due  to  thr  efferent  approximation  lor  the  positive  ion  drift 
vclr  c i ty  which  he  used : 

v (  .  +E/p( 1 -C  E/p)  E/p  ■  W 1 

r; 


(70) 


Cdta  was  ttl.en  froin  Reference  31. 
Data  was  calculated  1'n  this  study. 


V  =  k  {L/p)1/L':i-  -D—  )  E/p  -  w  O'D 

+  +  (E/p)1^  1 

Inhere  for  Hr  i<+  •••  8.E3  cn^-torr/vol t-sec , 

k4  -  4.1E4  cm  -ton ^/volt^-sec  ,  C  =  8.E-3  torr-cm/vol t , 

C  -  27.44  vol t‘/cmf -torr 2  and  1^  -  ‘dbv/cm-torr) . 

At  aGUv/cm-torr ,  Ward'?  approximation  yields  a  positive  ion  velocity 
that  is  only  1%  larger,  whereas  at  30  v/cm-ton  it  is  17?  larger. 

Since  the  slope  of  the  electric  field  is  inversely  proportional  to  the 
positive  ion  drift  velocity,  the  cathode  fall  length  and  voltage  are 
also  sensitive  functions  of  the  positive  ion  drift  velocity  ana  largely 
accounts  for  the  above  discrepancy.  Very  goto  agreement  is  reached 
with  Davies  ana  Evans  data  (Reference  32;  for  y  =  .23.  Care  shoula  be 
taker,  in  comparing  the  values  of  SFlY,kj)  since  this  'constant'  is 
proportional  to  (1  +  Y ) . 

As  can  be  seen  iri  Table  1.  for  He/HCl  mixtures  and  also  later  for 

Ar/HCl  and  Xe/EiCl  mixtures  as  well,  the  theoretical  values  for 

SF(y>E4)  do  not  compare  well  with  values  derived  from  experimental 

results.  This  is  aue  mainly  to  the  dttiriition  cf  V  and  d  and  the 

J  c  c 

great  difficulty  in  making  accurate  field  or  voltage  measurements  in 

the  cathode  fall  region.  Field  or  voltage  measurements  ’n  the  cathode 

rail  region  arc  exceptionally  difficult  due  to  the  high  fielos  and  high 

field  gradient.  V  art  c.  theoretically  are  define^  (  .  she  print  at 

which  the  electric  field  car,  be  extrapolated  to  zero.  In  order  to  compare  with 

theoretical  results,  accurate  measurement  of  dc  and  also  Vc 

regimes  trial  the  electric  field  in  the  cathode  f  1  be  mease  ed.  Only 

recent with  nonpet^orbinq  optical  measurement  techniques  h.is  arcirate 
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field  measurenients  in  the  cathode  fa1.1  reaion  been  fessible.  V  is 

c 

usually  measured  as  the  anode  is  moved  closer  to  the  cathode  while 
keeping  the  discharge  current  constant.  The  potential  difference 
across  the  electrodes  decreases  linearly  as  the  anocc  is  moved  through 
the  positive  column  region.  As  the  anode  is  moved  into  the  Faraday 
dark  space,  a  sudden  drop  in  potential  occurs  as  the  anode  fall 
disappears.  The  potential  difference  at  this  point  is  often 
experimental I3  defined  as  the  cathode  fall  voltage.  Further  movement 
ot  the  anode  through  the  negative  glow  causes  the  voltage  to  decrease 
very  little  due  to  the  very  low  field  in  that  region.  At  the  edge  oi 
the  cathode  dark  space,  the  voltage  begins  to  rise  rapidly  in  order  to 
maintain  a  constant  discharge  current.  The  cathode  fall  length  is 
usually  equated  to  this  determination  of  the  length  of  the  cathode  dark 
space.  The  experimental  definitions  of  V(  ano  d  are  thus  not 
consistent  with  the  theoretical  definitions,  however,  the  deflection 
of  an  electron  beam  shot  across  the  discharge  has  been  used  to 
determine  the  field  in  the  cathode  fall.  The  length  at  the  point  where 
the  field  is  extrapolated  to  zero  using  this  technique  has  agreed  well 
with  the  cathode  fall  length  determined  by  moving  the  anode.  This 
electron  beam  method  appears  to  vary  greatly  in  accuracy  ts  a  result  of  the 
design  of  the  apparatus.  Depending  upon  the  energy  spread  of  the 
electron  beam,  the  beam  could  also  affect  the  maintenance  condition  of 
the  discharge  and  thus  V  and  dc  it  the  beam  electrons  ionize  any  of 
the  gas  molecules.  Thus  throughout  the  remaining  results,  bF  (t,  k+) 
will  be  discussed  in  terms  of  other  theoretical  results  where  possible. 

F  i  ci;  re  6  shows  a  2%  decrease  in  the  electric  field  through 
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Figure  6  ELECTRIC  FIELD  IN  He/HCl  MIXTURES 


most  ot  the  cathode  fall  region  (0.-1 .26  cm)  as  the  amount  of  HC 1  is 
increased  from  0  to  5?.  This  is  contrary  to  the  results  that  might  be 
anticipated  from  Equation  66.  However,  HC 1  has  an  ionization 
potential  approximately  half  of  that,  for  he  (12.74eV  as  comparea  to 
24.69  eV),  the  HC 1  contributes  more  electrons  through  ionization  than 
the  discharge  loses  to  attachment  to  HC1  resulting  in  a  slightly 
reduced  cathode  fall  voltage  of  up  to  b%  for  the  95/5  He/HCl  mixture. 
There  is  also  a  trend  for  the  cathode  fall  width  ano  voltage  to  de¬ 
crease:  1%  and  1.4%  respectively  for  99/1  mixture  and  5%  and  6.3%  for 
the  95/5  mixture.  In  the  positive  column  region  (around  2.5  cm)  the 
electric  field  is  5%  higher  for  99/1  mixture,  and  25%  higher  for  a  95/5 
mixture  of  he/HCl  as  a  result  of  the  increase  in  the  number  of  attach¬ 
ing  molecules.  The  fields  in  the  positive  column  region  are  consistent 
with  the  value  of  the  field  required  to  have  the  electron  gain  and  loss 
processes  equal.  In  the  nearly  uniform  field  region  of  the  positive 
column  any  additional  loss  of  electrons  must  be  balanced  by  an  increase 
in  ionization  in  order  to  maintain  constant  current. 

Figure  7  shows  the  same  calculation  in  95/5  He/HCl  with  and 
without  the  contribution  of  ionization  from  HC1.  As  expected,  the 
electric  field  is  increased  in  the  cathode  fall  region  when  ionization 
of  HC 1  is  not  included.  This  additional  ionization  in  the  cathode  fall 
region  as  the  percentage  of  HC1  is  increased  results  in  a  small  rela¬ 
tive  increase  in  the  electron  current  density  in  Figure  8  in  the 
cathode  fall  region.  The  current  densities  plotted  in  Figures  8  and  9 
have  been  normalized  as  in  Equation  37.  The  electron  current  density 
decreases  in  the  positive  column  region  as  the  percentage  of  HC1  is 
increased  indicating  electrons  are  being  lost  to  the  formation  of  nega- 
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NORMALIZED  NEGATIVE  ION  CURRENT  DENSITY  IN  He/HCl  MIXTURES 


tive  ions.  Pz  expected  the  negative  ion  current  density  shewn  in 
Figure  9  becor.it s  a  larger  fraction  ol  the  discharge  current  otr.sity 
through  both  the  cathode  fall  region  and  the  positive  column  region, 
however  it  still  contributes  loss  than  1°  to  the  total  discharcr 
current. . 


Argon 

The  results  for  Ar/HCl  mixtures  will  be  presented  next  where 
p-ltorr,  d  =  1.5cm,  u  -  i.F-5  amp/cm?  and  <  ~  .04167.  A  summary  of 
numerical  values  for  the  following  figures  is  given  in  Table  2.  as 
was  done  for  He/HCl  mixtures. 

In  100%  Ar  the  value  of  SF(Y,  k  )  agrees  very  well  with  Ward's 
value,  within  1%.  Better  agreement  is  achieved  in  Ar  than  ir  (e 
because  the  slope  of  the  electric  field  is  steeper  in  Ar  resisting  in 
less  error  in  the  cathode  fall  length  and  voltage.  Approximately  the 
same  discrepancy  exists  in  the  positive  ler  drift  velocity  in  Ar  as  in 
lie.  At  400v/cn-turr ,  Ward's  approximation  in  Equation  76 
(where  for  Ar  u  =  1 .  FTcrr^-torr/ vol  t-sec  , 

k^  -  8.2513  cm5 -torr'/vol  t-sec ,  C  -  2.22E-3  t.orr-cm/vol  i, 

J  j  r 

0  =  86.52  vol t.‘/cnr -torr*  ,  and  W,  =  60v/cm-torr) 

i 

yields  a  positive  irr:  velocity  that  again  is  about  1W  larger,  whereas 
at  6fv/cm-tcrr  it  is  18*  larcer. 

In  contrast  to  He  mixtures  thert  is  a  negligible  f.T'ri  l  m  the 
electric  field  in  tiie  cathode  fall  region  (O.-.ffcpr  in  Fioi.ie  pj 
when  up  to  G1’  IC1  is  added  to  an  Ar  discharge.  li  this  case  since 
I  Cl  has  a  slightly  lower  icr.ization  potential  th..r  that  tor  f  r  M;;.7eV 
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able  l.  Summary  r.f  GLOW  Code  Results  for  Ar/HCl  Mixtures 


Data  was  taken  from  a  summary  presented  in  Reference 
Data  was  taken  from  Reference  85. 

Data  was  calculatec  in  this  study. 


ELECTRIC  FIELD  IN  Ar/HCl  MIXTURES 


as  compared  to  lb. 7c V ) ,  the  small  additional  ionization  aue  to  the 
small  percentage  cf  MCI  is  insufficient  tc  he  seen  as  any  effect  on  the 
electric  field  in  the  cathode  fall  region.  Figure  11  shows  the  same 
calculation  in  9b/ 5  Ar/HCl  with  and  without  the  contribution  of  ioni¬ 
zation  from  HC1 .  As  expected,  the  electric  field  is  increased,  but  not 
as  much  as  for  the  He/HCl  case  since  the  thresholds  for  ionization  art 
much  closer.  Similarly  no  change  is  seen  in  Figure  12  in  the  elec¬ 
tron  current  density  in  the  cathode  fall  region.  The  negative  ion 
current  density  in  Figure  13  never  becomes  a  significant  fraction  of 
the  total  discharge  current  density.  The  negative  ion  density 
increases  much  more  slowly  than  for  the  He  case,  since  the  fields  in 
the  cathode  fall  and  positive  column  are  higher  in  Ar  than  in  He. 

Ir.  the  positive  column  (around  1.2cm)  the  electric  field  is 
increased  approximately  3%  for  a  59/1  Ar/HCi  mixture  (not  plotted)  ar.d 
almost  12%  increase  for  a  95/5  mixture.  Again  this  is  the  result  of  an 
increase  in  the  number  of  attaching  molecules  in  the  positive  column 
region. 


Xenon 


The  results  for  Xe/HCl  mixtures  with  p  =  100  torr,  d  =  .01cm, 

J  r  . 1  amp/cm2,  and  Yr  .004  will  be  presented  next.  This  value  of  > 
was  chosen  so  results  coulo  be  compared  with  Ward's  re'ults  (reference 
05).  A  summary  of  numerical  values  for  the  following  figures  is  given 
in  ioble  3.  as  wac  r.orie  for  the  previous  mixtures. 

in  HUT  Xe,  there  is  (apparent  agreement  (0%,  between  the  •  cal  inc 
factors  calculated  from  Wards  data  ariG  the  data  from  the  GLOW  coot. 
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EFFECT  OF  IONIZATION  OF  HC1  IN  Ar/HCI  MIXTURES 


99%  fir 


NORMALIZED  ELECTRON  CURRENT  DENSITY  IN  Ar/HCl  MIXTURES 


Rr/HCI 


DISTRNCe  (cm) 

Figure  13  NORMALIZED  NEGATIVE  ION  CURRENT  DENSITY  IN  Ar/HCl  MIXTURES 


Approxi  rifely  the  Sdi.it  discrepancy  exists  in  the  positive  ion  drift 
velocity  as  in  the  previous  tart  gases.  At  690  v/cm-torr,  Ward's 
approximation  us1  ng  Enuation  76 

(where  in  Xe  a 4  -  4.F2  cm  torr/volt-cec, 

k+  =  4.E3  cm^-Torr2/vol t^-sec ,  C  -  ^.25F-J  torr-cm/vol t , 

L  -  225  vol t^/cm^-torr2 ,  anti  Wj  -  100  v/cm-tcrr) 

yields  a  positive  ton  velocity  that  is  about  lit  larger,  whereas  at  100 
v/cn-torr  ,  it  is  23"  larger. 

In  Xe  mixtures  there  is  a  very  slight  trenc  tor  the  electric  field 
to  increase  in  the  cathode  fall  region  (0.-.crnl  as  the  percentage  of 
HC1  is  increased.  This  is  observable  in  Moure  14.  Xe/HCl  mixtures 
are  unique  in  that  both  constituents  have  almost  the  same  ionization 
potential  (12.74eV  for  HC1  and  12.13eV  for  Xe).  Thus  only  effects  of 
attachment  are  observed  and  the  results  agree  with  the  analytical 
description  given  earlier.  The  effect  of  attachment  is,  to  raise  the 
field  slightly  throughout  the  cathode  fall  region. 

Lven  though  the  field  increases  as  HC1  is  c.oced  to  Xe,  the  elec¬ 
tron  current  density  decreases  slightly  throughout  most  of  the  cathode 
fall  region  in  Figure  15  due  to  the  pure  loss  of  electrons  to  at¬ 
tachment.  As  in  Ar  mixtures,  the  negative  ion  current  density  ir. 

Figure  16  grows  very  slowly  near  *-he  cathode  but  begins  to  grow 
rapioly  as  the  electric  field  levels  eft  into  the  positive  column.  It 
is  still  insignificant  in  comparison  to  the  rontribut ; or.  from  the 
electron  and  positive  ion  current  densities  to  the  tola'  discharge 
current  density. 

iri  the  positive  column  (.around  .00?  cm)  the  t  metric  fieV  in¬ 
creases  approx imately  ?°'  for  99/1  mixturp,  for  F5/5  mixture,  and  9°( 


65 


NORMALIZED  ELECTRON  CURRENT  DENSITY  IN  Xe/HCl  MlXTUkh 


for  90/10  mixture.  Aqair  this  is  the  result  of  an  increase  in  the  num¬ 
ber  of  attaching  species  in  the  positive  column  and  the  imposed  boun¬ 
dary  condition  of  constant  current. 

In  summary,  the  above  calculation?,  which  implicitly  assumed  that 
the  electrons  were  in  eauilibrium  with  the  local  field,  have  shown 
that : 

a.  The  tonnation  of  small  amounts  of  negative  ions  in  the 
cathode  fall  region  is  not  the  reason  that  the  cathode  a  length  con¬ 
tracts  in  the  axial  direction. 

b.  The  observed  contraction  of  the  cathoce  full  length  is  a  re¬ 
sult  of  the  increased  ionization  from  the  more  efficient  ionization  of 
the  attaching  gas  in  some  gas  mixtures.  The  contraction  will  only 
appear  in  gas  mixtures  in  which  the  ionization  rate  uf  the  attaching 
gas  is  larger  than  that  of  the  background  gas. 

c.  When  the  ionization  rates  and  thresholds  of  both  gases  in  the 
mixture  are  similar,  then  the  formation  of  negative  ions  tends  to  in¬ 
crease  the  voltage  drop  across  the  cathode  tall  region  and  to  increase 
the  field  at  the  cathode. 

d.  The  field  in  the  positive  column  reaion  increases  as  the 
percentage  of  attacher  is  increased.  TMs  is  required  in  order  tc 
balance  the  electron  gain  and  loss  processes  and  to  maintain  a  constant 
current. 
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SECTION  IV.  ELECTRON  KINETICS  IN  THE 


CATHODE  FALL  REGION:  NOMEOU IL I BR I UK  ANALYSIS 

1.  REVIEW  OF  NONEQUIL IBRIUK  ANALYSES 

The  analyses  in  Section  III  were  ell  based  on  Poisson's  equation 
and  current  continuity  or  flux  equations.  These  analyses  presuppose 
that  the  electron  energy  distribution  function  is  in  equilibrium  with 
the  local  electric  field.  As  mentioned  iri  Section  I,  there  is  analyti¬ 
cal  evidence  that,  the  electrons  are  not  in  ecu i  1  i bri uni  with  the  electric 
field  in  the  cathode  fall  region.  Since  the  electron  current  or  flux 
is  an  average  over  the  product  of  the  electron  number  density  times  the 
electron  velocity,  <r;v>,  it  is  not  possible  to  derive  n  or  v  from  this 
average  unless  the  nonequilibrium  electron  energy  distribution  function 
is  known.  There  are  two  methods  of  calculating  the  electron  energy  dis¬ 
tribution  function: 

1)  by  using  a  Monte  Carlo  simulation  of  the  electrons  as  they  drift 
from  cathode  to  aroce,  or 

L)  by  calculating  the  distribution  function  from  the  Boltzmann 
equation. 

This  chapter  will  critically  review  ano,  where  appropriate,  point 
out  detects  in  the  few  papers  that  have  investigated  the  nonequi 1 ifcrium 
electron  kinetics  usiny  either  of  these  two  methods.  The  new  results  of 
this  study,  which  usee  a  self-consistent  Boltzmann  method  and  included 
dissociative  attachment,  will  be  presented  at  the  end  of  the  chapter. 

First,  a  comparison  will  be  made  between  the  Monte  Carlo  method  and 
the  Boltzmann  method.  A1 thouyh  both  methods  should  yield  similar 
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rtsults,  th ert  are  impoi ■  i ant  computer  and  scientific  differences.  The 

Ronte  Carlo  method  is  much  more  costly  in  terms  of  computer  time.  In 

order  to  obtain  rood  statistics,  a  large  number  of  collisions  are 

5 

required,  typically  10  or  greater.  Ike  advantage  of  the  Korte  Carlo 
method  lies  in  its  ease  of  implementation  as  a  computer  coot,  however, 
the  Boltzmann  method  normally  enables  a  greater  insight  into  the  physics 
of  the  situation.  The  sensitivity  to  changes  in  inputs  is  more  easily 
investigated  primarily  due  to  the  shorter  run  time.  This  enables 
approximate  analytical  solutions  to  be  found  which  illustrate  the  Gepen- 
oer.ee  on  the  relevant  parameters.  In  uniform  electric:  fields,  the 
agreement  between  lionte  Carlo  results  and  Boltzmann  analyses  have  been 
excellent  (Reterences  62,  72,  77). 

The  theory  and  first  calculation  of  the  electron  distribution 
function  as  a  function  of  distance  through  the  cathode  fall  region  was 
reported  by  Allis  and  co-workers  (References  13,  14)  in  1975  and  1977. 
iheir  approach  involved  transforming  Boltzmann's  equation  into  energy 
spa.ee  where  the  new  variables  are  the  electron  potential  energy, 

=  ef  E(x)dx  and  the  electron  kinetic  energy  C  =  In  addition, 

they  treated  the  cathode  fall  region  as  a  separate  entity  allowing  for 
no  flux  of  positive  icr.s  from  the  negative  glow  and  forcing  the  elec¬ 
tric  field  to  zero  at  the  arede.  Their  calculations  were  based  on  an 
'ideal1  oar  defined  as  having  a  constant  momentum  transfer  cross 
section  and  constant  excitation  and  ionization  _ross  -  .’cts  with 
thresholds  of  oeV  and  16^V  respectively.  Ir  addition  their  calcula¬ 
tions  are  not  considered  so1! -cons  ister.t ,  meaning  the  calculations  of 
the  t  lectron  number  density  and  electric  field  we’rr  nut  iteram  .  on 
successively  to  converge  to  a  self-consistent  solution.  The  nonegui- 
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libriiin:  phenomena  winch  their  calculations  demons'll  ated  included: 

a.  A  d’p  in  the  electron  number  density  aojacent  to  the  cathode. 
This  results  from  the  electrons  Le i rig  accelerated  away  irom  the  cathode 
and  not  yet  having  sufficient  energy  fur  ionization. 

b.  A  peaked  -ionization  coett icier. t  which  is  snail  near  the  cath¬ 
ode  where  the  field  is  highest  ana  peaking  between  one  hall  and  two 
'birds  the  total  cathode  fall  length. 

ibis  technique  of  transforming  Boltzmann's  equation  into  energy  space 
has  become  the  basis  oi  several  investigations  aiscussed  later. 

Almost  simultaneously  with  Allis,  Tran  tiger,  et  al.  (Reference  £6) 
used  a  one  dimensional  Monte  Carlo  method  to  investigate  the  nonequi¬ 
librium  electron  kinetirs  of  the  cathode  fall.  They  began  by  assuming 
a  linearly  decreasing  electric  field.  Although  no  attempt  was  made  tc 
correct  the  field  through  Poisson’s  equation,  their  results  did  give  a 
good  physical  picture  of  the  nonequilibrium  process  occurring  in  a  high 
field  gradient.  They  found  the  electron  distribution  function  evolved 
from  a  sharply  peaked  distribution  near  the  cathode  to  a  multiply  peak¬ 
ed  distribution  entering  the  negative  glow.  This  was  the  first  theo¬ 
retical  prediction  of  the  three  groups  of  electrons  observed  experimen¬ 
tally  by  Pringle  and  Farvis  (Reference  73)  in  the  negative  glow  (high 
energy  electrons  which  have  undergone  few  collisions,  intermediate  cn- 
frgy  electrons  which  form  the  bulk  of  the  distribution,  and  low  energy 
electrons  which  were  primarily  a  result  ot  icrizatioi  r  ocesscs).  Tran 
hqr.c,  et  al.  concluded  th.it  the  use  01  Townsend's  ionisation  coeffi¬ 
cient  given  by  Equal u.r  41  was  inappropriate  ’n  the  cathode  tall  region 
‘’riie  the  electrons  never  reacheo  equilibrium  witl  the  field.  In  addi¬ 
tion  they  determiner  that  the  major  features  of  the  cathode  tail  are 
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actually  determined  by  a  ..mall  numbei  of  inelastic  col  1  i s i r. r. s  ana  that 
elastic  collisions  play  only  a  subsidiary  role. 

The  rrost  definitive  work  usiig  a  Holtzmann  analysis  was  accom- 
pl  i  she  a  ty  Long  (Reference  60).  It  was  based  upon  an  approach  out  1 'red 
by  Allis  (Reference  13)  of  a  direct  soittion  of  the  F'cltzmanr  equation 
without  the  conventional  expansion  methods.  Long  simplified  Allis'  re¬ 
sult  by  simplifying  the  angular  dependence  of  the  distribution  function 
in  the  cathode  fall  region.  I ; i s  method  of  solution  was  first  to  esti¬ 
mate  the  distribution  function  and,  tc  use  this  in  the  collision  terms 
to  integrate  the  distribution  function  alone  the  various  energy  charac¬ 
teristics.  The  distribution  function  was  iterated  urt.il  there  was  less 
than  .'e.%  change  in  the  distribution  function  at  each  energy  bin  and 
then  the  electron  number  densities  and  electron  n t rent  densities  were 
calculated  from  the  distribution  function  as  a  function  of  distance, 
lie  attempted  to  use  these  values  to  solve  Poisson's  equation  for  the 
electric  field  before  repeating  the-  process  to  converge  on  a  self-con¬ 
sistent  solution.  Long's  results  and  conclusions  for  mercury  and  argon 
agreed  very  well  with  those  of  Tran  figoc  for  hel’um.  However,  in  his 
integration  of  Prisson's  equation  for  the  electric  field,  he  failed  to 
include  any  change  in  the  ion  drift  velocity  due  to  changes  in  the 
electric  field  that  occurred  from  assuming  current  continuity.  This 
limited  the  convergence  oi  his  code  to  unique  cases  where  she:  anode  was 
placed  inside  the  caUmde  sheath.  T he  resolution  of  ti  o  problem  and 
further  improvements  of  Long's  computer  code  are  described  in  the  fol¬ 
lowing  section.  His  progran  required  it'i'h  of  mercery  to  load  and  rtn. 
Pun  tines  were  typically  5-10sec  per  f  ield  iterate  <•  or  a  CUC  LoOU 
compute  ", 
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Ml. i e  recently  Segur,  et,  al.  (ketererirp  79)  coupled  a  Boltzmarr  so¬ 
lution  fur  the  electrons  with  ar  .-nalyticul  approximation  to  the  Boltz¬ 
mann  equation  to*  positive  ions  in  the  cathode  fall  region  for  various 
hel  ium-rnercuty  mixtures.  They  assumed  the  electror  collisions  were 
isotropic  and  only  considered  charge  transfer  collisions  with  a  ccr- 
stant  cross  section  for  _he  positive  ions.  If.  did  not  have  to  use 
conservation  ct  current  to  calculate  t h.e  positive  ion  density  as  a 
function  of  position  as  Long  had  done.  Segur,  et  al.  initially  assumed 
a  linearly  decreasing  electric  field  vanishing  at  the  anode.  Within 
5-4  iterations  of  calculating  the  electron  and  ion  densities,  id  lowed 
by  solving  Poisson's  equation  for  the  field,  their  numerical  technique 
achieved  a  stable  convergence.  Thus  their  technique  was  self-consis¬ 
tent.  As  the  percentage  of  mercury  was  increased  from  OX  to  9.1%,  the 
peak  of  the  ionization  coefficient  increased  and  the  threshold  for 
ionization  shifted  towards  the  cathode.  This  is  due  to  the  fact  that 
mercury  atoms  have  a  larger  ionization  cross  section  as  well  as  a  lower 
ionization  threshold  than  helium  atoms.  They  found  the  ionization  in 
the  cathode  fall  region  took  place  in  two  stages:  mercury  is  ionized 
iirst  near  the  cathode  and  helium  is  ionized  slightly  further  away  from 
the  cathode.  Eventually  both  processes  exist  together  during  the 
remairiro  part  of  i  he  discharge.  Since  t tie  mobi  lity  of  the  mercury 
atoms  is  much  lower  that  hel  i  urn  atoms,  they  also  found  a  smaller  ionic 
current  neat  +  he  cathode  ar.G  a  larger  ionic  current  i.ea.  tin  ,  rede  when 
mercury  was  acute1. 

At  the  Some  time,  Boeuf,  et  al.  (Inference  L ' 1  investigated  the 
cathGCie  fall  in  a  helium  dischatye  using  a  se1  i  -cci  si  stent  Mor.e  Carlo 
teohr ioue .  They  also  investigated  some  features  in  an  oxygen  dis- 
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charge.  They  pc  rfonned  one  Kente  Carlo  calculation  for  positive  ions 
and  found  that  the  positive  i c r  'rift  velocity  reached  r nui  1  ibr-ium  with 
the  field  in  le  ss  than  HJ”C  or.:.  7 1  e>-eaf ter  ;l  ey  assured  that  the'  i t-r.s 

drifted  towards  the  cathode  with  their  equilibrium  drift  velocity,  in 
their  self-consistent  iterations  betwees  :lt  number  density  and  the 
field,  their  calculations  normally  reached  coi.-. .  rgence  after  5  itera¬ 
tions.  They  concluded  from  their  ionization  coefficient  and  electron 
drift  velocity  data  that  the  nonecui 1 ibrium  region  in  a  helium  dis¬ 
charge  extends  just  over  half  (54%)  of  the  cathode  fall  length  which 
agreed  with  the  previous  investigations. 

Z1  three  dimensional  Morte  Carlo  calculation  was  accomplished  very 
recently  by  Boeuf,  et.  al.  (Reference  201.  They  made  extensive  use  of 
the  rail  collision  technique  using  a  fictitious  total  cross  section  to 
maintain  a  pseudo-constant  collision  freroency  so  that  numerical  inte¬ 
grations  could  be  avoided  in  determining  the  time  interval  between  two 
collisions.  This  allowed  them  to  calculate  a.  three-dimensional  elec- 
trcr'  distribution  fundi  or  more  efficiently  takiro  angular  scattering 
into  account.  Radial  fields  were  not  included  and  no  attempt  was  made 
to  ’terate  successively  or  the  electron  number  density  and  electric 
field.  Thus  their  calculations  wer-e  not  self-consistent.  They  exam¬ 
ined  two  cases:  1  ’  a  normal  glcv.  discharge  in  helium  at  one  torr  with 
an  absorbing  anode  placed  in  the  regative  glow,,  and  2)  an  abnormal  dis¬ 
charge  in  helium  at  one  torr.  They  found  that  there  • s  little  error  if 
isotropic  scattering  is  assumed  arid  the  momentum  transfer  cross  section 
is  u sea  as  the  clastic  cross  section  irstead  of  assuming  angular  scat¬ 
tering  arc  using  the  elastic  differential  cress  sen;  ion.  fioeui  el  a7, 
shewed  that  the  original  assumption  of  forward  scattering  made  by  Tran 
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Ngoc  (Pe'erence  iv)  led  to  too  high  of  electron  energies  in  the  latter 
two  thirds  ot  the  cathode  fail.  As  a  resist  of  all  the  electrons  being 
directed  in  the  same  direction  ann  at  too  high  of  energy,  the  Townsend 
ionization  coefficient  for  the  forward  Mattering  case  was  as  much  at 
tlb'i  tee  low  throuch  most  ot  the  lot1  or  two  thirds  of  the  cathode  f a  1 1 
region.  They  colculateo  the  Townsend  ionizut  m  •  coefficient  by 
ii tegratirg  the  ionization  cross  section  ano  electron  distribution 
hunctioii.  Its  maximum  cecurred  somewhat  beyono  the  maxi  nun:  of  the  mean 
energy  distribution.  They  concluded  that  as  a  result  of  the  bimodal 
character  of  the  electron  distribution  function , the  mean  <  nergy  and 
ionization  coefficient  were  not  related  to  the  same  part  of  the 
di stribut ion  function.  Finally,  their  angular  distribution  functions 
show  that  the  high  energy  segment  of  the  distribution  function  is 
largely  forward  directed,  while  that.  ol  the  low  energy  part  is  much 
more  isotropic.  As  a  function  of  distance  through  the  cathode  fall 
region,  this  trend  became  much  more  evident  clo.t  to  the  negative  glow. 
This  concludes  the  critical  review  of  nonequ i i ibrii  i  ■  techniques. 

2.  DESCRIPTION  AND  MODIFICATIONS  TO  PRGGPAM  SHEATH 

The  ShLAiH  program  was  developed  by  W.  H.  Lorg  to  solve  the 
Boltzmann  equation  tor  electrons  in  ron-uniforr.i  fields  without  requir¬ 
ing  the  electron  distribution  function  to  be  in  equilibrium  with  the 
local  field.  The  collisioral  Boltzmann  equation  e'ese  ■'■ih-s  tit  evolu- 
.  ion  nf  a  swarr  of  electrons  in  both  space  and  time,  end  is  giver  h> 

!  '  •  ‘  -  ‘ L 
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where  f  -  f(r,v,t)  is  the  electron  distribution  function  which 

represents  a  complete  statisric.a 1  description  of  the  electron  in  a 

weakl>  ionized  gas.  The  variable  t  is  time,  v  is  the  electron  velocity 

vector  ,  r  is  the  electron  spat  ial  vector,  a  =  -e/m  (F  +  v  x  b)  is  the 

acceleration  of  the  electrons  doe  to  electric  and  magnetic  fields. 

;)fl 

at  col  represents  the  change  in  the  distribution  function  dee  to 
elastic  and  inelastic  electron  collisions. 

Considering  only  spatial  variations  in  the  direction  of  the  elec¬ 
tric  field,  then  78  reduces  to 

i»f\  (79) 


3f  ,  e  r  /  V  if 
v  —  +  —  E  ( x )  —  . 

X  lx  in  avx  Jt  /col 


Dividing  both  sides  by  v  ano  et(x),  Equation  79  can  be  transformed 

X 


in co  energy  space. 


if  +  Of  =  _  1 

3C  ev  E(x)  at  )  ol 


*) 


(80) 


where  e  J  E(x)dx  is  the  electron  potential  energy  and  t  -£nvx2  is 
;  h.t  electron  kinetic  energy  in  the  field  direction.  By  defining  the 
total  kinetic  energy  as  >  -  e  +  m,  then  f  <  ui  uc  ; cpi  esented  as  the  sum 
of  two  functions,  f  0  ,r.,7)  which  is  the  distribution  of  electrons 
with  vx<0  moving  against  the  field  and  f+  (<  , ;  , ■; )  which  is  the 
distribution  of  electrons  with  vy  >  0  moving  with  the  field.  The 
riable  represents  the  electron  rauooin  kinetic  energy.  Fractions 
7<3  and  8C  ran  be  used  to  mathematical  ly  define  equilibrium  ami 

nonequi  I  ibrioir  regions.  Equilibr  ium  ran  be  defined  as  occurring  when 

■  f  •  f  if  1  'f, 

,  <  ,  nr  when  ,  •  pv  r(x)  •  Thus  the  Mnetir  energy  gained 

X  C  0 1 

from  the  f  ield  is  b.ilarced  bv  the  loss  of  enere;  to  col  1  is.ion;.  In 
nonequi  1  ibrium  regions,  the  decrease  m.  tlie  potential  energy  of  the 
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electrons  is  converted  predominantly  into  an  'increase  in  their  kinetic 

ererqy.  In  Equation  79,  equilibrium  car.  similarly  be  defined  as 

f  o'  ( x )  if  j f  m  ;i f  \ 

occurnnc  when  v  M  nr  wren  '  trfTyV  .  In  other 

X  .X  r,  -vx  ,vx  Gf:ix/  col 

words,  equilibrium  exists  when  the  distribution  function  changes  faster 
in  velocity  space  thar  it  does  as  a  tunc  Lion  of  position.  Tins  was  the 
definition  used  in  Sectior  I. 

In  the  absence  of  collisions,  f  is  constant.  along  the  energy 
characteristics,  r  -  t  -  constant..  These  characteristics  correspond  to 
electron  trajectories  in  the  two-dimensional  energy  space  (£,<*).  The 
right  hand  term  ir.  Equation  80  represents  the  <  hanqe  ir  f  along  the 
energy  characteristics  when  collisions  are  included.  This  collision 


term  is  a  complicated  function  involving  i,  and  the  eross-sections  for 
all  the  various  collision  processes.  The  theory  and  derivation  of  this 
tern  is  described  in  Long's  report  (keie.ence  60:12-21).  He  showed 
that  80  can  be  put  in  operator  form 


jf  +  _a_? 

:>.j>  .it. 


(81) 


where  K  is  an  inteonl  operator  representing  a  sum  over  all  col  1  isicrial 
processes  and  Q(t  )  is  the  total  cross  section  tor  electrons  with  total 
kinetic  energy  .  Using  this  technique,  the  distribution  function  f 
does  not  need  to  be  stored  between  iterations,  only  Kf  must  be  re¬ 
tained.  The  Lollision  tern  ran  be  further  simplified  by  assuming  Kf  is 
nearly  isotropic  ir  velocity  space  and  expanding  Kf  in  legendre  poly¬ 
nomials,  P^(cos  0)  for  example. 


Kf  ^  E  A?(v)P  (coso) 

f  =0 


(82) 


where  A(v)  = 


(  --)  P  (cost)  d(coso) 
\  L/  ol  ' 


This  expansion  is  quite  different  troin  the  usual  two  tern:  expansion  cf 


hoi tzmann ' s  equation  in  spherical  harmonics.  This  expansion  is  an 


expansion  of  the  product  of  the  integral  operator  times  f,  anti  not  an 
expansion  of  Assuming  that  the  secondary  electrons  are  emitted 
i sotropical ly  and  neglecting  suferelastic  electron  collisions  and  the 


recoil  term  for  electrons.  Long  derived  the  loll  owing  form  of  the 


3f  Jf  1  1 

icf,  +  dr,  eE/NU)  coso 

2  e?.+i  {  ^  ) 

2.-0 

h 

—  f„('-+>.  \+)  Qh(>  +uh)  +  f  t  ftf  (c  '  ,<f  IQ,1  (e',G )de  1  1  P  (cose)  (83) 

h 

en-  . 

1 

-f(.  .£,*;{ q®(.  )  *  £  QS(i )  *  /  2 1  o’  (>  •  ')d,  ■  +  qJ(«)  +  ~  o(o 

1  h 

0 

where 

COSO  =  1  (--) 

f,(.  .0  -  -V1/ 

*  9 

f(>  cos  n ,f' )P  (coso)  d(cose) 

Q,  (■  )  =  4n  ~-J (cosy  )d(  COS.'  ) 

1  h  -  inelastic  collision  energy  loss 
i  1  -  ionization  energy  loss 

Q  -  collision  cross  secticr  integrated  over  all  scattering 
arigles 

superscript  h  refers  to  a  sum  ever  inelastic  processes 

superscript  i  refers  to  a  sum  over  ic  rnzatior  processes 

superscript  c  refers  to  a  sum  over  elastic  .  llision  pr;  cesses 

superscript  a  refers  to  an  attachment 

superscript  r  refers  to  recombination 

n4  -  pc  ’live  ion  density 


•  . 1 

>  a 


Thus,  ir  Equation  bi,  she  electrons  scattered-in  is  represented  by 

T— >  1  T  £  t  \  rfi  1  \  ,  \  +l  f  I  X  ^  ilfiS,  i. 


f|(t,',*)Qe(-  )  +  £  -y-  f  (>  +.  h: 

n 


+  / 


i  t  (,  ,.j.)Q(«  ,,  1  )d.  P  (i.os.>) 

r  it  ■ 


and  the  electrons  scattered  out  is  represented  hy 


Q(«.  )f=f( 


,£,:■){ q®<  )+  E  qJm  *  / 


2  Q’C.r ')*'«£(,.)♦£  Qq(-  )  } 


Legation  81  u-  solved  numerically  using  an  iterative  technique 
lor  evaluating  Kf.  It.  is  assumed  initially  that  there  is  no  contribu¬ 


tion  from  scattered  in  electrons,  ie  hi  -  0  and  Equation  83  is  then 


integrated  both  for  forward  moving  elector rs  and  backwura  moving 
electrons  for  each  value  of  r;=>-8 .  This  yielcs  a  first  approximatior. 
to  the  distribution  function  irrm  which  a  new  Kf  can  be  calculated. 


The  procedure  is  repeated  until  convergence  is  reached,  typically  when 
there  is  lets  than  .?"!■  choiige  in  the  distribut  icr.  function.  By  per¬ 
forming  the  integration  ir  a  specific  order  the  same  array  can  be  used 
to  store  both  the  old  and  the  now  values  of  Kf,  thus  minimizing  com¬ 
puter  storage  requi rements.  The  arrays  tor  Kf  anc  Kf  are  trape¬ 
zoidal  in  shape  since  f  ±  (>•,',* )  ••  0  for  6><J>+8  ,  where  t  is  the 

initial  kinetic  energy  in  the  iield  direction  of  the  electrons  leaving 
the  cathode,  figure  17  shows  how  the  arrays  for  Kf  are  arranged  in  core 
to  lurther  minimize  storage  requi  rements .  The  total  ,i  d  of  the  array 
in  Figure  17  is  48,240^  (136, 180^)  words.  For  a  2 CbV  cathode  fall, 
these  d  H.ioro  iens  gi-e  a.  ie\i  energy  ; evolution.  Ihe  storage  cf  only  Kf 
prever ■ s  the  distribution  for  c  1 1 on  from  being  know*  as  a  funct'on  of 
disicice,  except  at  the  cathode  and  anode  where  it  is  stored  in  sepa- 
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rate  arrays. 

boundary  conditions  are  usually  specified  at  the  cathode  and 
anode  At  the  cat  rode,  the  energy  distribution  of  electrons,  f  (.  ,0) 

is  a  variable  input.  The  functior  i  (  ,  ,  )  can  represent  the  dis¬ 
tribution  of  electrons  at  the  head  o*  the  positive  column  wher 

:  =  V  ,  the  vcltaqe  drop  across  the  cathode  fall,  negative  glow 

nia  x  c 

arid  Faraday  dark  space,  or  it  can  represent  the  distribution  of  elec¬ 
trons  at  the  anode  then  t  -  0,  or  it  con  represent  the  distribution 
along  an  axis  ut  symmetry  such  that,  t  -  f  +  .  Along  the  locus  of 
turning  points  where  the  electron  kinetic  energy  in  the  field  direction 
-  C,  the  forward  and  backwara  fluxes  are  equal,  f  (r ,0,<p)  =  f  (t,0,d). 

The  algorithm  used  in  Long’s  code  for  the  numerical  integration  of 
Fquation  L3  is  second  order  and  srabie.  The  distribution  function  is 


calculated  via 


+  C  (  ?  -  3  \ 

i  (  / _ _E/Ni  J  i 


Ac  E/N. 


where  C  =■  eE/N  {\)hl  and  Ac  is  the  energy  step.  The  relative  error  at 


each  step  is 


Error 


rel  2 


b.) 


which  implies  that  the  error  is  small  as  long  as 


Neat  the  cathode  this  con  :  i!.i(  r  is  easily  satisfice  since 
which  is  typically  LOeV  or  less. 

In  the  negative  9 1  cw  and  Ear-iday  hark  space  however  the  field  be¬ 
comes  <'ery  small.  In  this  region  m  the  discharge,  the  distribution 
function  is  changing  more  ^apidly  with  respect  to  position  (x)  than  it 
;S  with  respect  to  potertial  (  ;•  }.  Thus  ;  i  i.e  l"nget  an  appropriate 
variable,  are  /  becomes  the  variable  of  choice.  Charging  varieties. 
Equation  81  becomes 


--  +  eE(x)-'f-  =  N(x) 

.j  X  tw 


1/d 

( K-n(..  ))f 


(89) 


To  correctly  determine  r  through  these  regions,  the  code  first  starts 
with  coordinates  (■  0',;)  and  iterates  on  Equation  82.  When 
■  >L'L/tin,  then  the  coordinates  0 re  changed  to  (  ,  ,x)  anc  Equation 
89  is  i.ficr  iterated  on  until  convergence. 

The  cede  is  presently  configured  to  calculi 'e  the  electron  energy 
distribution  function  as  a  function  of  ti 4  stance  1  ><  n  ■  tie  cathode  tc  the 
ericde.  The  electrons  are  constrained  to  move  either  with  or  against 
the  field,  i.e.  all  integrations  are  in  the  plane  r>  =  0.  This  simpli¬ 
fication  significantly  reduces  computer  run  tin:e. 

The  electron  number  density  is  normal izea  to  the  electron  number 
drrsity  at  tht  cathode  such  that 


n  ( : )=n 
e  eo 


,:  )  +  f  (  ,:))> 


)  +  f  (>  ,1  )). 


0  0 


0  0 


(90) 


whr^e  subscript  0  indicates  the  boundary  condition  at  che  cathode. 
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Similarly  the  electron  current  density  is  normalized  tu  the  electron 
current  density  at  the  cathorii  end  is  giver  hy  the  following 


Townsend  rates  for  the  various  inelastic  processes  including  ionization 
are  normalized  to  the  electron  torrent  density  at  the  cathode  as  well 


In  order  to  obtain  a  completely  se it-consister 1  solution  ir  which 
Poisson's  equation  is  obeyed  in  addition  to  the  boundary  conditions  at 
the  cathode  and  anode,  successive  calculations  for  the  electron  dis¬ 
tribution  function,  electron  number  density  arid  the  resulting  electric 
^ i e 1 d  were  accomplished. 

long  attempted  to  do  this  for  the  case  where  the  anode  is  just 
i  side  the  cathode  fall-negative  glow  boundary,  by'  integrating  the 
(  ieclrun  d’Stribut  'or  lunctioi  if  obtained  the  election  number  dens  ty 
.  rd  electron  current  density  and  obtained  the  positive  irn  number 
density  usin<  current,  cei.strvat.ioir  . 
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(93) 


n  +  (x) 


_  J-Je(x) 

'  irr(x) 


v+(x)  =  vi+E(x)  (94) 

In  this  calculation,  however  he  had  to  use  hi;  initial  guess  for  the 
electric  field,  so  the  ion  drift  velocity  ( m  L )  was  net  consistent  with 
the  new  calculation  of  the  field. 

He  normalized  his  number  densities,  not  through  the  previous  defi¬ 
nitions,  but  by  the  ratio  of  E^/E^  where 


E  =  electric  field  at  the  cathode  prior  to  the  calculation 
o  r 

of  the  distribution  function 


E,  =  J(n,  -  n  )dx  which  is  the  electric  field  at  the  cathode 
i  X1  +  e 

after  the  unnornalized  densities  are  known,  is  the 
spatial  point  at  which  the  field  is  zero. 

This  approach  does  not  permit  the  secondary  emission  coefficient  to  be 
used  as  an  input  parameter,  however  the  coefficient  was  calculated  via 
the  following  equation: 


This  normalization  gc\e  good  resows  and  converged  as  lone:  as  the  anooo 
was  at  the  boundary  or  within  the  cathode  fall  region. 

Ji  was  found  by  this  study  that  Long's  technioi e  divergec  when  the 
anode  was  more  realistically  placed  outside  the  cathode  fa1!,  so  a  new 


technique  to  reach  a  self-consistent  solution  was  derived.  This 
technique  is  based  on  putting  Poisson's  eouation  in  cifferencc  form  so 


the  value  of  the  field  at  each  poim  is  based  on  the  value  at  the 
previous  increment  with  the  field  „t  i  ht  cathode  given  as  a  boundary 
condition.  Writing  Poisson's  equation  in  oiiierence  form  yields. 


E.=  E  .  .  -  —  (n.  +  n  .  ,-n  ,-n  .  -n  ,-n  .  , 
J  J-l  e  +J  +J-1  e/j  ej-1  -j  -j-1 


(x  ,-x  .  ) 

J  J-1 


(96) 


where  n  n  . ,  n  .  are  the  positive  icn,  electron,  and  negative  ion 

“J  "‘v/ 

number  densities  at  position  j,  and  x.  is  the  spatial  point  correspond- 

J 

ing  to  increment  j.  The  positive  icn  number  oei.;.’ t.y  can  be  written  in 
terms  of  the  current  conservation  such  that 


n 


+j 


ek  (E ,/p),/2 

•J 


(97) 


The  negative  ion  number  density  car:  be  calculated  from  the  negative  ior, 
current  density. 


A 

/ 


n  .  j  .dx  . 
J  ej  j 


(98) 


n~j  ek  (E7p)1/2 


(99) 


where  n  is  the  attachment  rate.  Equations  99  and  97  can  be 
substituted  into  96.  This  results  in  a  fifth  ordc*  polynomial  *  hose 
■elution  rar.  he  found  by  numerical  means. 
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X  + 


l?7? "  Ej-rc; 


-j-l 


) x"  +  a'^x  ■(  172  "  J-l  "  c  ‘y  j 
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where  a  = 


(J-Je.r2J-.i)p 
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c  =  —  ( n  .  +  n  .  .  ) 
■'  ej  ej-1 


E’. 


/  ■ 


Steinmarrs  nethod  (Reference  66 )  was  used  it  this  study  tc  solve  for  the 
roots.  Thus  E  is  calculated  based  on  n  ,  j  ,  j  aria  the  previous  value 
of  the  elec+ric  field.  This  technique  converged  after  successive 
iterations  cf  first  calculating  f(  -  ),  n  ,  i^  and  then  E(x).  Figure  18 
gives  an  example  of  the  convergence  as  the  eiectric  field  and  electron 
distribution  are  successively  iterated. 

lhe  initial  discharge  parameters,  Eq,  F^,  and  d^,  were  varied  to 
achieve  the  best  convergence.  These  parameters  describe  the  initial 
field  distribution  according  to  the  formulas 


E( x )  -  E  (1-  £)  x<d 
0  d  c 


(101) 


E(x)  =  E 


mi  n 


x  -  d 


(102) 


It  v.as  observed  that  the  optimum  convergence  was  achieved  when  the 
initial  L  . was  slightly  greater  than  (what  was  later  found  to  be)  the 
converged  ana  the  initial  slope  tQ/dL  was  equal  to  or  slightly 

shallower  than  the  shallowest  slope  of  the  converged  curve.  "T,is 
generally  was  the  slope  occurring  between  E^  and  3E  /4  or  7 Lq/8 . 

Figure  19  shows  the  effect  of  changing  the  value  o*  the  electric 
field  at  the  cathode.  Note  that  if  is  assumed  ‘or.  low,  in  this  case 
for  E  ^  757v/cn,  then  the  solution  to  Equation  ICC  becomes  imagi- 


CONVERGENCE  OF  THE  SHEATH  CODE  TO  THE  LOWEST  VOLTAGE 


ntry,  representing  an  unplysical  solution  for  the  field  between  the 
electrodes.  This  corresponds  qu'te  closely  to  the  phenomena  seen  in 
the  convergence  of  the  previous  GLOW  model  (Figure  4).  The  calcu¬ 
lation  for  Eo  =  758v/cm  was  deemed  to  Le  the  solution  most  representa¬ 
tive  of  a  discharge  since  it  represented  the  Irwest  voltage  in  which  a 
complete  real  solution  existed  from  cathode  to  .mode.  Although  E  was 
varied  manually,  a  technique  similar  to  the  modified  linear  interpo¬ 
lation  or  halving  the  interval  technique  could  be  programmed  to  con¬ 
verge  to  the  lowest  value  of  for  which  a  complete  solution  exists. 

Summarizing  the  changes  made  since  Long  publ  ished  his  report,  this 
study  has  added  attachment  as  an  electron  loss  process  in  the  numerical 
calculation  of  the  distribution  function,  modified  the  code  to  include 
a  variable  number  of  gases,  and  modified  the  code  to  include  a  variable 
number  of  collision  processes  per  gas.  The  original  program  was 
1  •’'rrited  to  analysing  a  pure  gas  with  one  momentum  transfer,  one  exci¬ 
tation,  arid  one  ionisation  cross  section.  The  'rownsend  ionization  and 
attachment  rates  are  now  calculated  from  the  dir ' ribution  function, 
instead  of  a  net  ionization  rate  from  the  electron  current  continuity 
equation.  In  addition,  in  the  calculation  of  the  positive  ion  drift 
velocity,  long  had  assumed  it  was  proportional  to  the  field.  This  has 
been  changed  so  that  the  positive  ion  drift  velocity  is  more  appropri¬ 
ately  represented  at  high  fields  by  Equation  44  in  which  it  is 
,, report i oral  to  the  square  root  of  the  field.  Also  a  new 
sel r -consistent  technique  was  found  for  calculating  the  electric  •f’eld 
from  the  electron  oer.r't.y,  the  electron  current  density,  and  the 
negative  lot  current  density. 

The  SHEATH  program  was  normal 1 )  run  on  a  Cyber  74  computer  system. 

A 


Each  iteration  to  converge  to  a  solution  for'  a  given  E^  took  or  average 
of  ISsec.  Normally  4  to  6  iterations  were  required  tc  reach  three 
significant  figure  convergence  in  E(d).  The  code  required  216K  of 
memory  to  load  arc1  run.  The  majority  u,:  t h ■* r  space  was  taken  up  by  the 
<c2G  x  TGI  array  described  ear1  ter. 

3.  RESULTS  OF  NUMERICAL  CALCULATIONS 

This  section  presents  the  results  of  calculations  using  the  SHEATH 
code  in  the  sane  rare  gas-HCl  mixtures  as  in  the  previous  chapter.  All 
the  data  presented  was  obtained  for  a  gas  pressure  of  1  ton .  Again 
the  results  for  he  will  be  be  presented  first,  followed  by  Ar  mixtures, 
and  then  Xe  mixtures. 

Only  the  cathode  fa1!  region  will  be  discussed  because  tie  anode 
was  placed  just,  outside  the  cathode  fall  region  where  the  electric 
field  leveled  off  to  a  constant  value  as  in  f.r.  obstructed  discharge. 

The  SHEATH  code,  while  improved,  still  had  or.  f  unity  converging  at 
very  low  electric  fields  where  the  coco  transformed  from  integrations 
over  energy  to  integrations  over  position,  ibis  transformation  was 
discussed  in  the  preceding  section.  All  the  SHEATH  data  preserved  in 
Ins  stuuj  was  obtained  before  this  transformation  of  variables 
eccurred,  ’e.  while  integrations  ate  conducted  in  energy  space. 

Further  comments  on  the  anode  fall  will  be  made  in  Section  V, 

Hel  lul  l 

This  section  presents  the  res i  Its  ot  the  nor, equilibrium  calcula¬ 
tions  for  He/HCi  mixtures.  A  summary  is  given  in  table  4  of  lie 
samt  sealing  parameters  listed  in  tables  in  the*  previous  chapter.  Ike 
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t  ne  Results  tor  fie/  HC 1  Mixtures 


Data  was  taken  f "on  Reference  of  (not  self-consister.l 2v  (not  self-consistent) 

Data  was  uken  from  Reference  79.  v  Data  was  taken  from  Reference 


cathode  fall  width  and  voltage  was  calculated  from  the  electric  field 
in  the  sane  manner  as  in  Section  III  for  the  GLCk<  model. 

There  is  a  wide  range  of  values  for  SF(y,  k+;  in  Table  4.  Some 
of  these  differences  are  due  to  the  ditierences  in  input  data,  such  as 
,  or  k+.  Also,  there  were  different,  approximations  used  for  the 
positive  ion  drift  velocity  dependence  on  the  electric  field.  However, 
there  were  also  important  model  differences  flat  influenced  the  results 
as  well.  The  data  from  calculations  by  Boeuf  and  Marode  (Reference  20) 
and  Tran  Ngoc  et  «1 .  (Reference  86)  are  not  self-consistent.  They 
picked  values  of  pd,  J/pz,  and  V  which  agreed  well  with  experiment  and 
assumed  a  linearly  decreasing  field  before  calculating  the  electron 
distribution  function  as  a  function  of  distance  through  the  cathode 
fall  using  a  Monte  Carlo  technique.  Segur  et  al .  (Reference  79)  calcu¬ 
lated  the  positive  ion  density  and  current  from  the  Boltzmann  equation 
for  ions,  so  an  analytical  fit  to  experimental  data  was  not  used  for 
the  positive  ion  drift  velocity.  Boeuf  et  al.  deference  21)  used  a 
Monte  Carlo  technique  for  both  electrons  and  positive  ions.  The 
positive  ion  "nonequi 1 ibr 1  urn  drift  velocity"  predicted  by  the  Monte 
Carlo  technique  was  within  20%  of  experimental  data  which  had  been 
measured  in  equilibrium.  In  the  SHEATH  coco,  Equation  44  was  useo  to 
calculate  the  positive  ion  drift  velocity.  This  analytic  approximation 
agrees  more  closely  with  the  drift  velocity  predicted  by  the  Monte 
Carlo  results  of  Boeuf  et  al.  than  the  experimental  data  they 
presented.  4s  a  result,  the  scaling  factor  SF(Y,k)  for  the  data  from 
the  Monte  Carlo  technique  and  the  SHEATH  coat  agree  within  1%.  This 
agreement  w'th  the  Monte  Carlo  technique  was  usee  .  .  a  benchmark  in 
comparing  the  accuracy  of  the  SHEATH  model  to  otho"  nonequi I ibrium 


techniques.  Note  that  'f  ,  .,  k4 ;  changes  very  little  for  the  different 

He/HCl  mixtures,  indicating  that.  tF(r,  k+/'  is  still  an  invariant  even 

when  the  ionizatirr,  efficiency  o4  the  gas  changes.  Ihese  changes  in 

the  ionization  coefficients  will  he  stir,  in  later  figures.  The 

observation  that  SF  (y,  k4 )  is  nut  a  luic’ion  of  the  ionization 

coefficient  is  consistent  with  Fouation  7b.  he  variation  in  SF  (y, 

k  )  in  this  author's  results  is  due  to  round  oil  and  extrapolation 

errors  ir.  pc?  and  V  . 

f  c  c 

As  in  Section  III,  the  theoretical  values  for  SF(y,k  ;  go  not 
compare  well  with  values  deriveu  from  experirrntal  results.  Again  this 
is  due  mainly  to  the  differences  in  the  theoretical  and  experimental 
definitions  of  V  and  dc  as  well  as  the  difficulty  ir'  making  precise 
tield  and  voltage  measurements  in  the  cathode  fall  region.  The  experi¬ 
mental  values  are  the  same  as  those  listed  in  Section  III.  The  SHEATH 
program  predicts  a  15i  decrease  in  the  electric  field  in  the  cathode 
fall  region  as  the  amount  of  EX  1  is  increased  i<,rrr  0  to  5%.  See  Figure 
?0.  This  corresponds  to  a  327  decrease  in  tin  rvthode  fail  voltage, 
this  trend  agrees  with  the  results  of  the  Glow  mouel  including  the 
linearity  of  the  field  near  the  cathode  in  both  cases.  The  shift  in 
the  electric  field,  however,  is  much  greater  in  the  ShEATH  model  arid 
can  be  observed  in  the  GLOW  model  in  Figure  6.  A  more  detailed 
comparison  of  the  two  techniques  will  be  described  in  the  next  chapter. 
Tins  decrease  ir  the  field  arc1  voltage  as  c  function  •  .f  Hr]  is  due  to 
the  increase  ir,  the  Townsend  ionization  coefficient  which  can  be  seer; 
in  Figure  21.  The  increase  in  the  peak  region  which  corresponds  t.o  the 
cathode  glow  region  is  due  to  the  fact  that  HC 1  ir  •  .,_h  easier  'o 
4c nize  than  he.  HC1  has  an  ionization  threshold  of  12.74eV  which  is 
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Figure  20  ELECTRIC  FIELD  IN  He/HCl  MIXTURES 


95* 


Figure  21  comparison  of  townsend  ionization  coefficient  in  Who  mixtures 


about  half  of  that  for  lie  whose  ionization  threshold  is  24.59eV.  Note 
that  the  spatial  threshold  for  the  Townsend  ionization  coefficient  is 
halved  when  HC1  is  added  corresponding  to  the  change  in  the  ionization 
cross  section  threshold.  This  indicates  that  very  close  to  the  cath¬ 
ode,  the  HC1  is  being  ionized  before  the  He.  These  results  are  consis¬ 
tent  with  the  results  ot  Segur,  et.  al.  with  Hg  in  He.  The  decrease  in 
the  linear  electric  field  with  the  addition  of  HC1  results  in  similar 
decreases  in  the  cathode  fall  length,  decreasing  6.4a  for  the  99/1 
mixture  and  21%  for  the  9b/ b  mixture.  These  results  are  summarized  in 
Table  4. 

The  curves  for  the  various  electric  fields  in  Figure  20  are  almost 
parallel  because  the  slope  of  the  electric  field  is  determined  through 
Poisson's  equation  by  the  positive  ion  density.  The  positive  ion 
density  at  the  cathode  is  determined  by  J,  and  the  positive  ion 
drift  velocity.  Since  these  variables  and  relations  remained  constant 
for  the  various  HC1  mixtures,  the  slopes  of  the  electric  fields  re¬ 
mained  constant. 

Figures  22  ana  23  compare  the  Townsend  ionization  and  attachment 
coefficients  as  a  function  ot  distance  for  59/1  and  95/5  concentrations 
of  lle/HCl.  The  magnitude  of  the  attachment  rate  is  much  smaller  than 
the  ionization  rate  throughout  the  cathode  fall  region  as  one  would 
expect,  except  immediately  in  front  of  the  cathode.  In  this  region 
which  correspond:,  to  the  Primary  dark  space,  the  elect  or  ns  have  just 
left  the  cathode  at  a  few  eV  of  energy  and  are  rot  yet  in  equilibrium 
with  the  field.  The  second  peak  seen  in  the  attachment  curve  at  l.c.r, 
is  probably  due  tc  the  electron  distribution  funct  cr  first  reaching 
the  ionization  energy  anti  producing  a  new  crop  of  slow  secondary  elec- 
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COMPARISON  OF  IONIZATION  AND  ATTACHMENT  COEFFICIENT  IN  99/1  He/HCl 


COMPARISON  OF  IONIZATION  AND 


trons . 


Figures  24  and  25  display  the  electron  and  negative  ion  current 
densities  as  a  function  of  distance  and  as  a  function  of  percent  of 
HC1.  The  increase  in  the  ionization  rate  also  manifests  itself  in  an 
increase  in  the  electron  current  density  throughout  the  cathode  fall 
region.  The  electron  current  density  stays  constant  for  the  first  few 
mean  ^ree  paths  as  the  electrons  are  being  accelerated  up  to  the  ioni¬ 
zation  potential.  It  then  increases  almost  exponentially  through  the 
remainder  of  the  cathode  fall.  As  expected,  the  negative  ion  current 
density  in  Figure  25  increases  as  the  proportion  of  HC1  increases.  It 
is  characterized  by  a  very  steep  onset  due  tc  the  large  attachment  rate 
near  the  cathode.  As  the  total  distribution  of  electrons  is  accelerat¬ 
ed  away  from  the  cathode,  they  reach  the  ionization  threshold  of  the 
gas.  Most  electrons  are  then  past  the  peak  region  of  the  attachment 
cross  section  and  the  negative  ion  current  density  grows  at  a  much 
slower  rate.  Note  the  logarithmic  scale,  indicating  the  growth  of  the 
negative  ion  current  density  is  still  slightly  faster  than  an  exponen¬ 
tial  function  of  distance. 

Figures  26,  27,  and  28  display  the  electron,  positive  ion,  and 
negative  ion  number  densities  as  a  function  of  distance  through  the 
cathode  fall  region.  The  result  of  nonequilibrium  electron  kinetics  is 
easily  seen  in  these  figures.  The  electron  density  initially  decreases 
because  the  electrons  are  rapidly  accelerated  away  from  the  cathode. 
Just  before  they  have  gained  sufficient  energy  for  ionization,  their 
density  reaches  a  minimum  and  begins  to  increase.  The  fact  that  by 
eld-cathode  fall  the  field  has  also  decreased  substantially  fcl'ows  the 
electron  density  to  increase.  Similarly  the  positive  ior  density  in- 


101 


DISTANCE  (cm) 

Figure  24  ELECTRON  CURRENT  DENSITY  IN  He/HCl  MIXTURES 
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Figure  25  NEGATIVE  ION  CURRENT  DENSITY  IN  He/HCl  MIXTURES 
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iyure  27  CHARGED  PARTICLE  DENSITIES  IN  99/1  He/HCl 
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CHARGED  PARTICLE  DENSITIES  IN  95/5  He/HCl 


creases  away  from  the  cathode  since  there  is  very  little  ionization 
close  to  the  cathode  and  the  positive  ions  are  being  accelerated  towards 
it.  Although  negative  ions  are  present  throughout  the  cathode  fall 
region  in  figures  27  and  28,  their  cerisity  is  generally  much  too  smal: 
to  affect  the  net  space  charge  to  any  significant  degree.  Note  that  the 
negative  ion  density  is  more  than  four  oraers  of  magnitude  smaller  than 
the  positive  ion  density.  The  positive  ier  density  actually  determines 
the  slope  of  the  field  throughout  the  cathode  fall.  The  magnitude  of 
the  electron  density  approaches  that  of  the  positive  ion  density  only 
after  the  negative  glow  is  reached  or,  in  these  examples  just  prior  to 
the  anode.  Equilibrium  is  never  reached  in  these  figures  because  a 
slightly  bimooal  electron  distribution  still  exists  at  the  anode. 

Since  the  cathode  fall  voltage  decreases  as  HC1  is  a deed,  the  ions 
arrive  at  the  cathode  at  a  lower  velocity.  Similarly,  the  electrons  are 
rot  accelerated  away  as  quickly,  sc  both  electron  and  positive  ion  dens¬ 
ities  increase  near  the  cathode  (20?  and  8?  respectively  for  up  to  5? 

HC1 ) .  The  opposite  trend  car  be  seen  in  Figures  20,  27,  and  28 
near  the  cathode  glow-negative  glow  boundary  (1.5-1.75  cm).  The  elec¬ 
tron  Gensity  increases  as  a  function  of  HC 1  in  this  region  due  to  t.he 
increase  in  the  ionization  coefficient  closer  to  the  cathode.  Since 
the  field  is  lower  in  th’s  region  and  the  ionization  coefficient  is 
lower,  the  positive  ion  density  decreases,  as  the  percentage  of  HC T  is 
increased,  which  in  a  self-corsistent  (or  interrelated)  manner  causes 
the  cathc.oe  tall  region  to  contract.. 


The  results  for  Ar/HCl  mixtures  are  summarized  in  Ta^e  5 
as  was  done  for  previcis  mixtures.  As  in  He,  there  is  a  wide  range  of 
values  for  SF(y,  k+)  in  Table  5  for  Ar.  The  only  other  nonequili¬ 
brium  data  available  for  Ar  is  that  cf  long  (Reference  60)  wirch  was 
not  entirely  self-consistent.  The  lack  of  sel t -consistency  of  his 
method  was  discussed  earlier  in  this  Section.  Also,  he  had  assumed  the 
positive  ion  drift,  velocity  was  proportional  to  the  field.  This 
assumption  is  only  appropriate  in  the  low  field  region,  i.e.  E/p  < 
60v/cm-torr.  Note  that  as  in  He  mixtures  the  scaling  factor  SF(v,  k+) 
remains  almost  constant,  as  the  concentration  of  I  Cl  increases. 

In  contrast  to  He  mixtures,  there  is  a  negligible  effect,  on  the 
electric  field  in  the  cathode  fal1  region  when  up  to  5%  HC1  is  added. 
This  can  be  seen  in  Figure  i9  and  these  results  agree  with  the  Gl OW 
model .  This  electric  field  distribution  throuoh  the  cathoae  fall  cor¬ 
responds  to  approximately  a  184v  potential  drop  across  the  cathode  fall 
region.  With  1%  HC1 ,  the  difference  in  total  ior’r.itiun  is  too  small 
to  have  any  observable  effect  on  the  electric  field.  The  ionization 
thresholds  for  Ar  (15. 7eV )  and  HC1  (12.7eV)  are  much  closer  than  for  He 
arid  FC1.  This  difference  is  too  small  to  be  observed  in  Figure  3C  in 
the  Townsend  icr.izat  on  coefficient.  With  5%  HC1 ,  there  is  sufficient 
F.C1  to  increase  the  ionization  rate  about  1.6%  throughout  most  of  the 
cathode  fall  region.  Figures  31  and  32  compare  the  T  wr^end  iuniza- 
tion  ana  attachment  coefficients  as  a  function  of  distance  for  99/1  and 
95/5  concentrations  in  Ar/HC1 .  As  in  He/hCl  mixtures,  the  attachment 
rate  is  larger  than  the  ionization  rate  only  for  a  ‘ ew  eiectror  colli¬ 
sion'  rinse  to  the  cathode  where  the  electrons  have  insufficient  energy 
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for  iorizat  ion. 


Tirures  33  and  z4  display  1  he  electron  end  negative  ion  current 
3ens.it.ier  as  a  fi.net  ior  of  distance  and  concu  tration  cf  MCI .  Again 
very  liuh  charge  is  observed  in  the  ilectron  current  density  since 
the  change  in  HU  concentration  hau  only  <.  minor  effect  on  the  ioniza¬ 
tion  coefficients  in  Figure  31.  lhe  electron  current  density  stays 
constant  as  the  electrons  are  accelerated  up  to  the  ionization  poten¬ 
tial.  Again  the  negative  ion  current  density  in  Figure  34  is  charac¬ 
terized  by  a  very  steep  onset  lol lowed  by  a  slightly  faster  than  expo¬ 
nential  growth  through  most  of  the  cathode  tali  region.  Even  though 
tht  negative  ion  current  density  increases  with  HC1  concentration,  its 
contribution  to  the  total  discharge  current  through  the  cathode  fall  is 
•-till  negligible. 

lhe  same  nonequilibrium  phenomena  that  occurred  in  He  mixtures  in 
the  number  densities  can  be  observed  in  Ar  mixtures  in  Figures  36,  36 
and  37.  The  positive  ion  density  increases  in  the  direction  moving 
away  from  the  cathode,  while  the  electron  number  density  initially 
decreases.  The  number  of  negative  ions  is  more  than  two  orders  of 
magnitude  smaller  than  the  number  of  positive  ions,  and  actually  their 
tensity  does  not  begin  tc  grow  exponentially  until  the  cathode  fall 
begins  to  me  rye  into  the  negative  glow. 

Fince  there  is  very  little  difference  in  the  fields  adjacent  to 
liie  cathode,  there  is  very  little  difference  in  the  electron  and  posi¬ 
tive  inn  number  densities  in  this  region.  The  positive  icr  density  is 
■..till  several  orders  of  magnitude  larger  than  the  electron  art  negative 
ion  densities  and  thus  determines  the  slope  of  ihe  electric  field. 

Only  in  the  last  third  of  the  cathode  fall  length  (.5-. 76  cm)  is  there 
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Figure  33  ELECTRON  CURRENT  DENSITY  IN  Ar/HCl  MIXTURES 
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a  charge  when  Figures  35,  3b,  and  37  are  compared.  As  in  t.e/HCl 
mixtures,  the  electron  number  tensity  increases  and  the  positive  ion 
number  uensity  decreases  but  to  a  lesser  decree  since  there  is  ies' 
change  in  the  Townsend  icrization  coefficient. 

Xenon 

The  results  Tor  Xe/HCl  mixtures  are  summarized  in  Table  6.  There 
H  s  very  Iffiie  data  for  comparison  in  Xe.  The  cathode  fall  voltage  V  . 
rd  length  are  sensitive  function.  ot  the  secondary  emission 
coeff  icient  v.  The  value  of  i-  .0U4  was  chosen  not  to  compare:  with 
an  iron  cathode  but  to  be  able  to  compare  with  the  tcui Ibrium  results 
rf  the  GLOW  code  which  were  compared  with  Ward's  oata.  Ever  though 
this  value  of  v  is  low,  it  can  not  totally  account,  for  the  targe 
discrepancy  between  theory  arid  experiment  in  SF(Y,k+).  lhe  main  source 
or  error  is  probably  in  the  value  used  for  k  .  The  value  of  was 
obtained  trori  Ward  (Reference  93:L791)  who  subjectively  fit  Equation  74 
to  experimental  data.  He  made  r:o  systematic  effort  to  obtain  the  best 
parameters.  In  addition,  only  the  first  term  in  Equation  74  was  used 
by  the  'HEATH  mooel.  In  effect.,  a  smaller  positive  ion  drift  velocity 
would  remit,  in  a  steeper  slope  for  the  electric:  field.  This  also 
weulc  tend  to  make  pd  and  V  smaller.  The  minor  variation  in  SI  Iy.L,! 
in  this  author's  results  is  due  to  round  off  and  extrapolation  errers 

in  pc  and  V  . 

c  c 

In  contrast  to  He  ar.o  Ar  mixture'.  there  is  a  slight  increase  in 
‘he  electric  field  in  fhr  cathode  fall  region  when  ip  to  5X  Hll  is 
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<  rider  1 1  s  ireno  o'  Me  field  «  l<  tiling  vertically  car  It  .ten  m  i  iy 
ua  oh.  It  result 5.  n  t.  slight  increase  ir  th  voltage  whuh  is  sun 
i  >a  ••  i ;  f  d  in  Tah  1 1  c.  ibis  trend  i*  also  the  .-dire  as  in  Ft  giro  13  f  rei 
the  Cl  TV  r  ceel .  Slightly  hiaher  fic'os  arc  nan  iced  tc  ft  tain  the  sane 
yet  t  on  i /at  ion  eea>  i  he  c.athode.  ihis  result  t  ;t.  o  un  nun  vertical 
shit'4  *  jit  fc. }  ec  r  r  c  field  throughout  toe  cathode  fa1!  region.  hole 
that  even  though  'he'  field  n  shifted  io  slightly  h  icier  values,  there 
i c  r-  observable  change  in  the  lonizaiien  coefficient  in  the  cathode 
fall  rector,  as  illustrated  in  Figure  it).  Again  cnly  near  the  minimum 
tielc  region  in  ftr  rtgative  glov.  is  a  difference  observed.  Xe/HC! 
r.iixti  res  are  unique  in  that  both  constituents  have  almost  the  seme 
.or  ration  potential  ( it.  7  4  for  MCI  anc  121. 13  tor  he)  and  siir.’lar  cress 
se  :  c<  nacnitudes .  f  oescr iptior  cf  these  cross  sections  occurs  in 
Appendix  f.  Thus  tht  s:ight  expansion  seer;  litre  is  probably  due  to  the- 
■formation  v  ■*  negative  i o r ;  close  to  the  cathode  where  there  is  a  peak 
•  i  tiie  attachment  rate.  f  loss  of  electrons  close  i,ci  the  cathode  is 
rimiiar  tc  the  effect  ed  reducing  the  secondary  emission  coefficient  cf 
the  cathode.  A  decrease  in  .  sindlarly  resu’ls  in  a  higher  cathode 
tai'  w ‘’facie.  However,  the  slope  of  the  electric  field  stays  the  same 
Ht  cause  the  positive  iot.s.  still  doiriraie  the  otitr  number  densities. 
Thus  the  'rend  of  the  field  shift  rr  vertical l>  with  an  increase  in 
■  (  T  r '  <  '.(i.  c  oe  expt  i  fed  to  err+inue  tor  ioiyer  concert  rations  ct  i-ll 
i  /t .  ir  Li.  •  *  1 >  dy  the  pt  .  itive  lor,  velocity  of  or  i>  the  ran'  g<=  s 
.1  :r.c  iv.ci  d.  For  Xr  'K  mixture  the  add  it  -.or.  ot  Ho!  wou.d  ’end  to 
rrrcasr-  J !.»  ,  verug-  ior  drift  ve  loci'.-  .  This  woiMc;  result  ir; 

a  uhclicv.fr  'lope  f  fir  c.ectru  field  arc'  ,  ’  i, .  uer  expansion  if'  pc  . 
r  ;<  u>  «"•.  40  ar  r  t  \  compare  lie  Towitserr*  i<  ri /.atioi  c.r.d  attachmord 


COMPARISON  OF  IONIZATION  AND  ATTACHMENT  COEFFICIENT  IN  99/1  Xe/HCl 


(ttO/I) 


31HH 


DISTRNCE  (cm) 

Figure  COMPARISON  OF  IONIZATION  AND  ATTACHMENT  COEFFICIENT  IN  95/5  Xe/HCl 


coefficients  as  a  function  of  distance  for  59/1  and  95/5  concentrations 
in  Xe/HCl .  Again  the  attachment  coefficient  is  much  smaller  than  the 
ionization  rate  throughout  the  cathode  fall  region  except  immediately 
in  front  of  the  cathode.  Note  that  the  region  where  the  attachment 
coefficient  is  greater  than  the  ionization  rate  is  smaller  in  Xe  mix¬ 
tures  than  in  He  or  Ar  mixtures.  This  is  due  to  Xe  having  a  lower 
ionization  threshold  than  the  other  buffer  gases.  The  first  peak  in 
the  attachment  coefficient  is  a  result  of  the  slow  electrons  leaving 
the  cathode  and  the  second  peak  is  a  result  of  the  formation  of  slow 
electrons  that  have  lost  their  energy  to  excitation  or  ionization.  The 
first  peak  has  the  effect  of  reducing  the  secondary  emission  coeffic¬ 
ient  of  the  cathode. 

Figures  42  and  43  display  the  electron  and  negative  ion  current 
densities  as  a  function  of  distance  and  concentration  of  HC1 .  Very 
little  change  is  observed  in  the  electron  current  density  until  near 
the  cathode  glow-negative  glow  boundary.  Here  the  electron  current 
density  is  observed  to  decrease  instead  of  increasing  as  a  function  of 
HC1  concentration  as  in  previous  rare  gas  mixtures.  This  is  another 
example  of  the  effect  of  attachment  in  the  cathode  fall  region.  Even 
though  the  negative  ion  current  density  in  Figure  43  increases  signifi¬ 
cantly  with  HC1  concentration,  its  contribution  to  the  total  discharge 
current  is  still  negligible. 

The  same  rionequi  1  ibrium  phenomena  that  occurred  in  previous  rare 
gas  mixtures  for  the  number  densities  can  be  observed  in  Xe  mixtures. 
The  results  are  ilUistrated  in  Figure  44  through  46.  Again  the 
trends  of  each  o*  the  number  densities  is  similar  to  the  trends  seen 
before.  The  positive  ion  density  is  still  several  orders  of  magnitude 
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larger  than  the  electron  and  negative  ion  densities  and  thus  determines 
the  slope  of  the  electric  field  in  Xe  as  well.  The  number  of  negative 
ions  is  more  than  three  orders  of  magnitude  smaller  than  the  number  ot 
positive  ions,  and  actually  their  density  does  not  begin  to  grow  expo¬ 
nentially  until  the  cathode  fall  begins  to  merge  into  the  negative 
glow.  As  in  Ar,  there  is  almost  no  difference  in  the  number  densities 
of  electrons  and  positive  ions  adjacent  to  the  cathode,  since  the  field 
changes  only  slightly  there  when  HC 1  is  added.  As  in  Ar,  only  in  the 
last  third  of  the  cathode  fall  length  (.45-. 65cm)  do  the  number  densi¬ 
ties  differ  as  a  function  of  HC1.  This  is  a  result  of  a  change  in  the 
net  ionization  prior  to  this  region. 

In  summary,  we  have  examined  three  cases  represented  by  He/HCl , 
Ar/HCl  and  Xe/HCl  mixtures.  In  all  three  gas  mixtures,  there  was  a 
peak  in  the  Townsend  attachment  coefficient  as  a  result  of  electron 
nonequilibrium.  However,  throughout  the  cathode  fall  region  the  posi¬ 
tive  ion  density  dominates  the  other  charge  densities  and  so  determines 
the  slope  of  the  electric  field  even  in  electronegative  gas  mixtures. 
He/HCl  mixtures  were  characterized  by  a  significant  contraction  of  the 
cathode  fall  length  and  a  decrease  in  the  cathode  fall  voltage.  This 
was  a  result  of  the  HC1  being  easier  to  ionize  than  the  lie.  As  a 
result  of  HC1  having  a  lower  ionization  cross  section  threshold  than 
He,  the  threshold  for  the  total  Townsend  ionization  coefficient  shifted 
closer  to  the  cathode.  There  was  a  negligible  effect  in  the  electric 
field,  cathode  fall  lergth,  or  cathode  fall  voltage  when  HC 1  was  added 
to  Ar  mixtures.  The  Townsend  ionization  coefficient  did  increase 
slightly  with  the  addition  of  5?  HC1 .  This  is  because  the  ionization 
rate  is  very  sensitive  to  the  gas  species  and  the  shape  of  the  electron 
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energy  distribution  function.  In  contrast,  in  Xe/HCl  mixtures,  the 
electric  field,  cathoae  fall  lergth,  and  cathode  fall  voltage  increased 
slightly,  probably  as  a  result  of  the  nor.equi  1  ibrium  electrons  being 
attached  to  form  negative  ions  in  the  primary  dark  space. 

In  conclusion,  the  formation  of  negative  ions  in  the  cathode  fall 
region  does  not  lead  to  a  contraction  of  the  sheath.  The  contraction 
that  is  often  observed  and  that  was  seen  by  Emeleus  and  Sayers,  was 
actually  due  to  an  increase  in  the  ionization  efficiency  of  the  gas 
mixture  when  an  electronegative  gas  is  added. 


SECTION  V.  COMPARISON  OF  RESULTS  OF  THE  EQUILIBRIUM 
AND  NONEQUILIBRIUM  ANALYSES 


This  section  compares  and  contrasts  the  GLOW  and  SHEATH  models. 
First,  their  equations  ana  assumptions  are  briefly  reviewed.  Second, 
the  results  from  these  models  are  compared.  In  order  to  determine  the 
influence  of  norequi 1 ibrium  without  further  complications  due  to  nega¬ 
tive  ions,  results  for  the  pure  rare  gases  will  be  emphasized.  These 
results  include  the  electric  field,  Townsend  ionization  and  attachment 
coefficients,  the  electron  current  densities,  and  the  charged  particle 
densities.  The  trends  predicted  by  both  models  are  summarized  along 
with  a  criticism  of  each  technique  on  its  applicability  in  the  cathode 
fall  region.  Finally,  the  applicability  of  each  model  to  the  anode 
region  is  discussed. 

1.  REVIEW  OF  BASIC  ASSUMPTIONS 

The  GLGW  and  SHEATH  models  are  founded  on  two  entirely  different 
methods  for  describing  the  kinetics  arid  fields  in  a  glow  discharge.  As 
described  in  Section  III,  the  GLOW  model  is  based  on  a  self-consistent 
solution  to  a  set  of  equilibrium  equations:  Poisson's  equation,  the 
hydrodynamic  current  continuity  equations  for  electrons  and  negative 
ions,  and  the  current  conservation  equation.  The  SHEATH  model  des¬ 
cribed  in  Section  IV  is  based  on  a  solution  of  the  Boltzmann  transport 
equation  for  electrons  as  a  function  of  distance  coupled  with  Poisson's 
equation  and  the  current  conservation  equation.  The  equations  are 
listed  in  Table  7  for  comparison.  The  fundamental  difference  behind 
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Comparison  of  Basic  Equations 


these  equations  is  that  the  glow  model  assumes  the  electrons  are  in 
equilibrium  with  the  field,  whereas  the  SI  1  HATH  model  does  not.  Another 
difference  exists  in  the  negative  ion  current  continuity  equations. 
Negative  ion  losses  due  to  ion-ion  recombination  were  not  included  in 
the  SHEATH  model  because  recombination  was  Tound  to  be  very  small  in 
the  GLOW  model.  Since  both  methods  are  describing  the  same  discharge, 
the  boundary  conditions  in  Table  7  are  almost  identical. 

li  addition,  both  methods  predict  results  that  appear  to  scale 
according  to  Equation  74.  In  the  equilibrium  mooel ,  limited  data 
was  taken  at  10  and  100  Torr  to  check  the  pressure  scaling.  There  was 
also  good  agreement  with  other  theoretical  data  when  the  same  value  cf 
the  secondary  emission  coefficient  and  the  same  relationship  of  the 
positive  ion  drift  velocity  to  the  electric  field  was  used.  Ir,  the 
nonecui  I  ibrium  model,  the  scaling  factor-  in  Table  4  for'  1 00%  He 
agreed  very  well  with  the  Monte  Carlo  results  of  Boeuf,  et  al.  (Refer¬ 
ence  21).  Scaling  agreement  of  both  models  can  be  expected  because 
both  models  support  the  assumptions  used  in  deriving  Equation  74,  which 
are : 

a.  Positive  ions  dominate  the  other  charged  particle  species 
through  most  of  the  cathode  tall  and  therefore  determine  the  slope  of 
the  electric  field. 

b.  The  positive  ion  drift  velocity  is  proportional  to  the  sauare 
root  of  the  electric  field  ut  high  fields. 

r.  The  electric  field  can  be  approximated  as  a  linearly  decreas¬ 
ing  function  through  most  uf  the  cathode  fall  region. 

Respite  the  differences  in  assuming  equilibrium  or  not  assuming 
equilibrium,  botti  models  scale  apprepri ately  and  predict  the  same 


136 


trends  iri  all  three  gas  mixtures. 


2.  COMPARISON  OF  RESULTS 

Comparing  the  results  of  the  (.LOW  and  SHEATH  models  from  Sections 
111  and  IV  respectively  indicates  that  both  methods  do  predict  the  same 
trends  for  each  gas  mixture.  In  He  mixtures  both  models  predict  a 
contraction  of  the  cathode  fall  viith  the  addition  of  HC1  as  seen  in 
Figures  5  and  20.  Data  from  both  models  indicate  tins  contraction  is 
due  to  the  more  efficient  ionization  of  the  HC1  in  the  mixture,  rather 
than  the  formation  cf  negative  ions  in  the  cathode  fall  as  was  specula¬ 
ted  by  Emeleus,  et  al.  (Reference  42).  Similarly  in  Ar/HCl  mixtures, 
both  methods  agree  that  the  addition  of  up  to  5°;  MCI  has  an  insignifi¬ 
cant  effect  on  the  cathode  fail  field  or  voltage.  This  can  be  observed 
in  figures  9  and  29.  In  Xe/hCl  mixtures  a  slight  increase  in  electric 
field  as  well  as  voltage  drop  across  the  cathode  fall  was  predicted  by 
both  models.  This  increase  can  be  seen  in  Figures  13  and  38.  Data 
from  both  methods  agree  that  this  field  enhancement  is  due  to  the  loss 
of  electrons  to  attachment  near  the  cathode,  although  this  is  perhaps 
most  easily  seen  from  the  theory  described  in  section,  "New  Equilibrium 
Analysis  oi  the  Cathode  Fall  Region  including  Negative  Ions"  in  Section 
III.  ltns  loss  of  electrons  to  attachment  close  to  the  cathode  is 
analogous  to  reducing  the  electron  secondary  emission  coefficient. 

Figures  47  through  4^  compare  the  distribution  of  the  electric 
field  through  a  discharge  lor  He,  Ar,  and  Xe.  In  each  figure,  the 
slope  of  the  electric  fields  are  approximately  the  same  as  a  result  of 
the  same  fielc  relationship  being  useo  for  the  positive  i or.  critt 
velocity.  Also,  the  results  of  the  SHEATH  mode!  are  consistently 
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higher  than  those  of  the  GLOW  model  until  the  field  begins  to  merge  in¬ 
to  the  negative  glow.  This  is  a  result  of  the  nonequi 1 ibrium  behavior 
of  the  electrons  near  the  cathode.  The  best  agreement  is  obtained  for 
Xe  ano  the  worst  +or  He,  since  the  electrons  have  shorter  mean  free 
paths  in  Xe  than  in  He  and  thus  spatially  reach  equilibrium  earlier. 

The  Townsend  ionization  coefficient  for  each  model  is  plotted  in 
figure  50  along  with  the  electric  field  calculated  by  the  SHEATH  model. 
These  coefficients  were  calculated  for  the  given  electric  field  using 
the  equations  in  Table  7.  As  can  be  seen,  the  plots  for  these 
coefficients  are  not  the  same  shape  due  to  the  r.onequi  1  ibrium  electron 
behavior  in  the  cathode  fall  region.  In  the  SHEATH  model,  the  elec¬ 
trons  require  several  mean  free  paths  before  they  reach  the  ionization 
energy  and  can  begin  to  multiply.  This  spatial  delay  in  ionization 
results  in  a  higher  cathode  fall  voltage  and  higher  field  at  the 
cathode.  The  ionization  coefficient  predicted  by  the  SHEATH  model 
actually  lags  the  field,  ie.  the  rate  is  low  where  the  field  is  highest 
and  then  overshoots  and  is  high  where  the  field  merges  into  the  nega¬ 
tive  glow.  The  GLOW  model  assumes  the  electrons  to  be  in  equilibrium 
with  the  field  and  thus  predicts  the  highest  ionization  at  the  highest 
fields  adjacent  to  the  cathode.  The  cathode  fall  voltage  required  to 
produce  sufficient  ionization  for  the  discharge  maintenance  condition 
(Equation  1)  is  correspondingly  lower  for  the  equilibrium  model  than 
for  the  nonequilibrium  model. 

The  Townsend  attachment  coefficient  for  each  model  is  similarly 
plotted  in  Figure  51  along  vjith  the  electric  field  calculated  by  the 
SHEATH  model.  These  coefficients  were  calculated  for  the  given  elec¬ 
tric  field  using  the  equations  in  Table  7.  Again  the  plots  for  these 
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Figure  51  COMPARISON  OF  ATTACHMENT  COEFFICIENTS 


coefficients  are  not  the  same  shape  Cue  to  the  nonequi 1 ibrium  electron 
behavior  ir  the  cathode  tall.  The  electrons  close  to  the  cathode  have 
a  low  energy  and  undergo  several  collisions  before  reaching  the  ioni¬ 
zation  energy  of  the  gas.  This  results  •in  a  large  peak  in  the  attach¬ 
ment  coefficient  close  to  the  cathode  in  the  resists  predicted  by  the 
SHEATH  model.  The  Townsend  attachment  coefficient  predicted  by  the 
SHEATh  model  actually  leads  the  field,  i.e.  the  rate  is  high  where  the 
field  is  highest  and  then  undershoots  ano  rises  again  where  the  minimum 
field  merges  into  the  negative  glow.  Note  that  the  attachment  coeffi¬ 
cient  predicted  by  the  nonequilibrium  model  is  more  than  an  order  of 
magnitude  smaller  than  that  predicted  by  the  equilibrium  model  except 
immediately  in  front  of  the  cathode.  Thus  the  ionization  and  attach¬ 
ment.  coefficients  predicted  by  the  ec;ui 1 ibrium  model  do  not  represent 
at  all  the  actual  spatial  distribution  of  the  coefficients  predicted  by 
the  SHEATH  model . 

The  growth  of  the  electron  current  density  is  portrayed  for  each 
of  the  rare  gases  in  Figures  52  through  54.  The  GLOW  model  consistent¬ 
ly  predicts  a  faster  rising  electron  current.  This  difference  is  again 
a  result  of  the  nonequil ibrium  behavior  of  the  electrons  close  to  the 
cathode.  The  electrons  first  have  to  be  accelerated  to  the  ionization 
potential  of  the  gas  before  they  can  begin  to  multiply.  This  causes 
the  electron  current  density  predicted  by  the  SHEATH  model  to  lag  the 
current  density  predicted  by  the  GLOW  model.  To  make  the  ecuilibriun, 
model  agree  more  closely  with  the  noriequi  1  ibrium  model ,  Tran  Ngoc,  et 
al.  (Reference  8 b)  have  suggested  that  the  mean  electron  energy  and  the 
Townsend  ionization  coefiicient  be  related  through  an  energy  balance 
eauation.  This  yielded  good  agreement  with  their  nonequilibrium 
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Figure  52  COMPARISON  OF  ELECTRON  CITRRENT  DENSITIES  IN  100?:  He  FROM  CLOW  AND  SHEATH  METHODS 


figure  53  COMPARISON  of  ELECTRON  CURRENT  DENSITIES  IN  1007,  Ar  FROM  CLOW  AND  SHEATH  METHODS 
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COMPARISON  OF  ELECTRON  CURRENT  DENSITIES  IN  100?:  Xe  FROM  CLOW  AND  SHEATH  METHODS 


ionization  coefficient  calculated  by  a  Monte  Carlo  technique.  Simi¬ 
larly,  relating  the  Townsend  ionization  coefficient  to  the  mean 
electron  energy  should  also  predict  a  slower  rising  electron  current 
density  more  closely  approximating  the  nonequilibrium  results. 

At  the  cathcce  tall-negative  glow  boundary  which  is  annotated  by 
x=dc  in  Figures  52  through  54,  the  electron  current  density  was  about 
90%  of  the  total  discharge  current  density  (88$  for  He,  89$  for  Ar,  and 
95%  for  Xe  in  the  SHEATH  model;  85%  for  He,  85%  for  Ar,  and  86%  for  Xe 
in  the  GLOW  model).  This  implies  that  a  larger  positive  ion  current 
density  exists  at  this  boundary  than  at  the  Faraday  dark  space-positive 
column  boundary.  As  a  result,  positive  ions  must  be  formed  in  the  neg¬ 
ative  glow  by  electron  collisions  from  the  few  remaining  fast  electrons 
that  have  managed  to  be  accelerated  across  most  of  the  cathode  fall 
region.  This  supports  Druyvesteyn  and  Penning 's  (Reference  38)  hypo¬ 
thesis  and  Ul'yanov's  (Reference  87)  calculations  indicating  a  signifi¬ 
cant  number  of  positive  ions  do  enter  the  cathode  fall  region  from  the 
negative  glow.  Thus,  Von  Engel  and  Steinbeck's  (Reference  90)  boundary 
condition  that  the  positive  ion  current  density  approaches  zero  at  the 
cathode  fall-negative  glow  boundary  is  not  a  good  assumption. 

Another  example  of  the  inability  of  the  equilibrium  model  to 
describe  the  nonequilibrium  behavior  in  the  cathode  fall  region  is 
shown  in  the  plots  of  the  charged  particle  densities  in  Figures  55 
through  57.  There  is  gooci  agreement  between  the  methods  on  the  posi¬ 
tive  ion  density  near  the  cathode.  The  GLOW  results  agree  within  13% 
for  He,  6%  tor  Ar,  and  5%  for  Xe  of  the  positive  ion  density  predicted 
by  the  SHEATH  model  at  the  cathode,  ihis  agreement  in  positive  ion 
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Figure  56  COMPARISON  OF  POSITIVE  ION  AND  ELECTRON  DENSITIES  IN  1 OOX  Ar  FROM  GLOW  AND  SHEATH  METHODS 
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Figure  57  COMPARISON  OF  POSITIVE  ION  AND  ELECTRON  DENSITIES  IN  100%  Xe  FROM  CI.OW  AND  SHEATH  METHODS 


density  results  in  the  electric  tields  being  almost  parallel  in  Figures 
47  through  4S.  however,  not  as  good  agreement  is  obtained  as  the  cath¬ 
ode  tail  transitions  into  the  negative  glow  region.  This  is  due  to  the 
different  rates  of  growth  of  the  electron  density.  The  GLOW  method 
consistently  predicts  an  electron  density  that  is  higher  than  the 
SHIATH  method  even  close  to  the  cathode.  At  the  cathoae  where  the 
electron  current  densities  agree  for  both  methods,  the  electron  density 
calculated  by  the  GLOW  method  is  larger  than  the  SHEATH  results  because 
the  electron  drift  velocity  calculated  by  the  GLOW  model  is  lower  than 
that  calculated  from  the  distribution  function  in  the  SHEATH  model. 

Ihe  difference  between  the  electron  number  densities  calculated  by  the 
GLOW  and  SHEATH  models  does  not  aftect  the  slope  of  the  field  until 
near  the  negative  glow  region  because  the  electron  density  is  5-4 
orders  of  magnitude  less  than  the  positive  ion  oensity.  As  a  result, 
the  equilibrium  theory  is  found  not  to  be  capable  of  describing  the 
particle  densities  as  a  function  of  distance  through  the  cathode  fall 
region.  Through  adjustment  oi  boundary  conditions  such  as  the 
secondary  emission  coefficient  and  the  electron  drift  velocity  at  the 
cathode  better  agreement  with  experimental  results  could  probably  be 
obtained,  but.  it  would  still  not  represent  the  actual  nonequi  1  ibrium 
kinetics  involved  there. 

The  macroscopic  parameters  such  as  discharge  current,  cathoce  tall 
voitage  and  length,  and  electric  field  each  contain  one  or  more  aver¬ 
ages  or  integrals  over  the  electron  distribution  function.  It  was 
thought  that  these  integrals  might  tend  to  dilute  the  effects  of 
nonequilibrium  on  these  macroscopic  parameters.  Thus  good  agreement  on 
the't  rncicroscopu  parameters  should  be  expected  between  the  GLOW  and 


SHEATH  models.  This  was  found  not  to  be  the  case.  Table  8  summarizes 


these  macroscopic  parameters  for  the  GLOW  and  SHEATH,  models.  The 
voltage  and  thus  the  electric  field  at  the  cathode  are  determined  by 
the  net  multiplication  of  the  electrons  as  they  cross  the  cathode  fall. 
For  an  electronegative  gas.  Equation  ?.  for  the  multiplication  of  the 
electrons  becomes 


M  =  exp  ^C(o.(x)-r,(x)  )dx  (103) 

o 

From  the  results  of  the  nonequilibrium  analysis  in  Section  IV,  the 
attachment  coefficient  is  always  much  smaller  than  the  ionization 
coefficient  except  immediately  in  front  of  the  cathode.  The  attachment 
rate  for  the  cases  examined  was  only  a  small  perturbation  on  M  in 
Equation  103.  This  also  agrees  with  the  results  presented  for  Xe/HCl 
mixtures  where  the  ionization  coefficient  effectively  did  not  change  as 
a  function  of  gas  mixture. 

Returning  to  Figure  50,  the  area  under  the  ionization  curves  are 
different  for  the  GLOW  and  SHEATH  models.  The  GLOW  model  assumed  the 
electrons  were  in  equilibrium  with  the  field  and  results  in  the  elec¬ 
trons  being  able  to  ionize  immediately  upon  leaving  the  cathode.  The 
SHEATH  model  allows  the  electrons  to  not  be  in  equilibrium  with  the 
field.  According  to  the  nonequilibrium  theory,  several  electron 
collisions  occur  before  the  majority  of  the  electrons  have  attained  the 
energy  necessary  tor  ionization  of  the  gas.  Thus  the  nonecui 1 ibrium 
model  predicts  a  higher  cathode  fall  voltage  than  the  equi 1 ibrium  model 
in  order  to  maintain  the  same  net  multiplication. 

Analyzing  this  through  the  mathematical  relationships,  the  slope 
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of  the  electric  neld  is  determined  by  the  mobility  (k  )  of  the  posit¬ 
ive  ions  through  Poisson's  equation.  This  is  because  the  positive  ion 
current  density  at  the  cathode  is  determined  by  the  boundary  conditions 
at  the  cathode  involving  the  secondary  emission  coeificient  and  the  to¬ 
tal  discharge  current  density.  The  length  of  the  cathode  fall  is  de¬ 
termined  when  the  multiplication  of  the  electrons  is  sufficient  to 
maintain  the  discharge  and  equals  1  +  l/f(Equation  1).  The  slope  of 
the  electric  field  and  its  extrapolated  x-axis  intercept  (d  )  then  de¬ 
termine  the  voltage  for  the  cathode  tall.  Thus  the  difference  in  the 
macroscopic  parameters  in  Table  £  is  a  result  of  the  nonequilibrium 
character  of  the  electrons  which  alters  the  ionization  coefficients 
from  those  predicted  by  the  equilibrium  (GLOW)  model.  Thus  any  agree¬ 
ment  of  the  equilibrium  method  with  experiment  is  entirely  fortuitous 
as  von  Engel  suggested  (Fveference  89:224). 

The  scaling  factor  SF(i>k+)  listed  in  Table  £  should  be  identi¬ 
cal  for  both  the  SHEATH  and  GLOW  models  since  they  both  used  the  same 
secondary  emission  coefficient  and  same  field  relationship  for  the  pos¬ 
itive  ion  drift  velocity.  The  slight  differences  in  the  data  between 
the  models  are  primarily  due  to  extrapolation  errors  in  determining  dc 
and  to  round  off  errors  in  d^  and  V  . 

Although  the  inability  of  the  CLOW  model  to  describe  nonequi libri¬ 
um  effects  accounts  for  most  of  the  disagreement  in  the  preceding 
figures  and  in  Table  8,  another  factor  affecting  the  agreement  of  the 
two  methods  is  the  use  of  two  different  data  bases.  The  GLOW  model 
bases  its  electron  kinetics  cn  curve  fits  of  empirical  data  for  the 
Townsend  ionization  ana  attachment  coefficients.  These  coerficients 
are  normally  derived  from  driff  tube  experiments  where  a  swarm  of  elec- 
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t ror> b  drift  between  two  electrodes  ird  are  in  equil  ibrium  with  the 
field.  The  SHEATH  model  ori  the  other  hand  bases  its  electron  kinetics 
on  scattering  cross  sections.  Under  ideal  conditions  both  determina¬ 
tions  of  t  ,  either  from  oritt  tube  experiments  or  cross  sections,  should 
be  consistent.  Some  attempt  was.  rude  in  Appendix  C  to  correlate  the 
sets  ot  cross  sections  used,  with  measured  transport  properties  using 
an  eoui 1 ibriuin  Boltzmann  code.  It  was  impossible,  however,  to  corre¬ 
late  them  over  the  whole  range,  since  transport  data  was  lacking  at 
some  high  E/p's.  In  addition,  the  equilibrium  Boltzmann  code  was  limit¬ 
ed  to  calculations  at  low  and  medium  fields  because  it  did  not  include 
the  newly  created  secondary  electrons  in  the  distribution  function. 

Thus  the  data  bases  for  E/p's  less  than  about  l(>0v/an-torr  dittered 
less  than  20?  but.  for  higher-  E/p  values  the  difference  could  be  more. 

In  summary,  both  the  SHEATH  and  GLOW  models  preoicted  the  same 
trends  in  He/HCT,  Ar/HCl ,  and  Xe/HCl  mixtures  and  both  models  scale 
according  to  Equation  74.  However,  the  discussion  in  this  chapter 
has  shown  that  the  GLOW  model  should  not  be  used  to  describe  any  of  the 
cathode  fall  parameters.  In  the  cathode  fall  region,  the  electrons 
have  been  shewn  not  to  be  ir  equilibrium  with  the  Held.  As  a  result, 
the  cathode  fall  voltage  and  length  are  larger  in  the  SHEATH  model  than 
the  GI0W  model  in  order  to  compensate  tor  the  nonionizing  region 
(Primary  bark  Space)  near  the  cathode  where  the  elect  i  ons  are  being 
accelerated  up  to  the  ionization  energy  of  the  gas. 

This  investigation  of  the  catlioae  fall  mien  sets  the  stage  tor 
further  investigations  or  the  negative  glow.  Due  tc  the  slightly 
bimodal  character  of  the  electron  distribution  function  at  the  cothoGe 
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fal  1  -nega  t  i  ve  glow  boundary  founc!  in  the  SH  EAT  I  i  results,  a  nonequi  1  ib- 


rium  technique  should  also  be  usea  to  investigate  the  negative  glow. 

In  this  region,  a  small  number  of  electrons,  that  have  managed  to  tra¬ 
verse  the  cathode  fall  with  one  or  no  collisions,  lose  their  energy  in 
the  negative  glow  forming  new  electrons  and  positive  ions  through 
ionization.  In  essence,  the  electron  distribution  function  is  still 
changing  faster  in  coordinate  space  than  it  is  in  velocity  space. 

3.  ANODE  FALL  REGION 

This  section  discusses  the  applicability  of  both  methods  to  the 
sheath  region  at  the  opposite-  end  of  the  discharge  iear  the  anode.  In 
the  GLOW  model  calculations  the  anode  was  placed  after  a  chosen  length 
uf  positive  column.  However,  due  to  a  limitation  of  memory  space  the 
anode  for  the  SHEATH  cases  was  placed  in  the  negative  glow.  This  cor¬ 
responds  to  an  "obstructed  discharge1’  (Reference  6).  The  properties  of 
the  anode  are  very  different  when  it  is  placed  in  a  slightly  negatively 
charged  region  such  as  the  Faraday  dark  space  as  compared  to  being  in 
contact  with  a  neutral  plasma  such  as  the  positive  column.  These  prop¬ 
erties  are  discussed  further  below. 

The  boundary  conditions  for  the  anode  fall  are  fairly  simple  in 
the  more  usual  case  when  the  anode  tall  and  positive  column  are  adja¬ 
cent.  Gutficient  electrons  must  be  accelerated  across  the  anode  fall 
and  create  sufficient,  ionization  such  that  the  number  of  positive  ions 
crossing  the  anode  fa  1 1 -pos i ti ve  column  boundary  is.  equal  to  the  number 
crossing  the  positive  column-! araday  dart,  space  boundary. 

It  is  commonly  believed  that  each  electron  entermc  the  anode  fei« 
must  be  accelerated  to  the  ionization  potential  of  the  gas.  Fach  elec- 


tron  has  to  create  an  ion  before  striking  the  anode  in  order  to  main¬ 
tain  charge  neutrality  lr.  the  positive  column.  It  was  found  experi¬ 
mentally  (data  is  summarized  in  References  I  and  6)  that  the  anode  lall 
voltage  was  usually  about  the  same  magnitude  as  the  ionization  poten¬ 
tial,  except  in  el ectroneaati ve  gases.  However,  it  is  also  known  that 
the  electrons  in  the  positive  column  have  an  equilibrium  distribution 
of  energy,  with  a  few  having  uti  energy  greater  than  the  energy  equiva¬ 
lent  to  the  ionization  potential.  Thus  the  potential  is  different  for 
each  electron  to  be  accelerated  to  the  ionization  potential  of  the  gas. 
In  addition,  charge  neutrality  ar.o  uniform  density  in  the  positive 
column  are  maintained  by  requiring  that  an.  ion  enter  the  Faraday  dark 
space  each  time  one  enters  the  positive  column  from  the  anode  region. 
Similarly  an  electron  is  required  to  enter  the  positive  column  from  the 
Faraday  dark  space  each  time  ore  leaves  the  positive  column  and  enters 
the  anode  fall  region.  Using  these  arguments,  it  is  not  necessary 
under  ideal  conditions  for  the  anode  fall  voltage  to  always  be  equal  to 
or  slightly  larger  than  the  ionization  potential. 

the  anode  acts  likes  a  probe  collecting  the  random  positive  end 
negative  currents  from  the  positive  column  and  so  its  area  is  also  very 
important.  If  tire  anode  area  is  so  small  that  the  negative  random  cur¬ 
rent  is  less  than  that  required  by  the  external  circuit,  then  a  region 
of  space  charge  lorms  to  pull  more  negative  charges  from  the  positive 
column  and  make  the  collecting  area  of  the  anode  larger.  This  reauires 
a  potential  on  the  order  ui  the  ionization  potential  or  less,  since  a 
potential  much  iaryet  than  the  icrizatior.  potential  would  tend  to 
create  numerous  electron-ion  pairs  reducing  the  effectiveness  of  the 
sheath,  /‘•node  tail  volt, roes  smaller  than  the  ionization  potential  hav1- 


been  observed  (Reference  6)  using  large  concave  anodes.  This  type  of 
anode  reduces  the  losses  to  diffusion  and  enables  more  negative  charges 
to  be  collected  from  the  positive  column. 

The  boundary  conditions  for  the  anode  fall  are  more  complicated 
when  the  anode  is  in  the  Faraday  dark  space.  Near  the  negative  glow 
boundary  there  may  be  sufficient  fast  electrons  from  the  cathode  fall 
to  carry  most  of  the  discharge  current.  These  electrons  are  already  of 
sufficient  energy  that  little  or  no  acceleration  is  required  for 
ionization  to  produce  positive  ions.  In  this  case  no  anode  fall  exists 
and  the  field  remains  almost  constant  right  up  to  the  anode.  As  the 
electrode  separation  is  increased,  the  anode  fall  region  must  of  course 
eventually  transition  into  the  normal  type  previously  described. 

Francis  (Reference  6)  claims  a  normal  anode  fall  develops  suddenly  at  a 
critical  distance  from  the  cathode  lying  near  where  the  Faraday  dark 
space-positive  column  boundary  would  normally  occur.  This  critical 
distance  should  be  approximately  equal  to  the  stopping  distance  in  a 
gas  for  an  electron  with  an  energy  approximately  equal  to  the  cathode 
tall  voltage. 

Other  factors  described  by  Francis  (Reference  6)  which  affect  the 
voltage  across  the  anoae  fall  are: 

a)  stepwise  ionization.  Stepwise  ionization  results  in  a  smaller 
voltage  since  the  electrons  do  not  have  to  be  accelerated  to  the  full 
ionization  potential  in  order  to  create  an  ion. 

b)  the  electron  current  density.  This  is  closely  related  to  the 
anoae  geometry.  Increasing  the  current,  density  beyond  what  the  anode 
can  normally  carry  is  tantamount  to  decreasing  the  anode  area  ard  hence 
increases  the  anoae  voltage.  At  high  currents,  regularly  spaced  spots 


form  on  the  anode  which  ate  regions  of  high  conductivity. 

c)  electronegative  gases.  Adding  an  electronegati ve  gas  to  an 
electropositive  gas  may  or  may  not  cause  the  anode  voltage  to  increase 
tor  the  same  reasons  Lite  cathode  tall  nay  or  may  not  contract  as  dis¬ 
cussed  in  previous  chapters.  The  potential  drop  across  the  anode  would 
depend  upon  the  partial  pressure  averaged  ionization  potential  of  the 
gas  mixture  in  addition  to  the  loss  of  electrons  to  negative  ion  forma¬ 
tion.  Generally,  a  high  coeh icient  of  electron  attachment  greatly 
reduces  the  ionization  efficiency.  For  example,  anode  fall  voltages  in 
the  halogens  have  been  observed  to  be  several  times  the  ionization 
potential . 

d)  artificially  produced  space  charge.  Specially  coated  anodes 
which  yield  positive  ions  upon  electron  impact  obviously  have  smaller 
anode  voltage.  If  the  number  of  positive  ions  is  more  than  that 
required  for  the  flux  entering  the  positive  column,  the  fall  can  be 
negative. 

t )  proximity  of  side  walls.  Changing  the  distance  of  the  side 
walls  can  change  the  rate  of  diffusion  and  raaial  fields  which  car 
result  in  increasing  or  decreasing  the  anoae  fall  voltage. 

The  anooe  fall  in  the  GlOW  code  calculations  was  adjacent  to  a 
positive  column  region.  The  anode  fall  voltage  in  Figures  5 ,  9, 
and  13  were  in  the  range  of  1  -  Tv .  These  values  were  determined  to 
be  the  difference  in  voltage  at  the  anode  arid  at  the  minimum  electric 
field,  ibis  is  much  smaller  thorn  the  normal  anode  tall  voltage  which 
is  or  the  order  ui  the  ionization  potential  of  the  gas.  The  use  of 
Ward's  formula  tor  tin  drift  velocity  (Keference  9^)  could  lead  to  a 
more  correct,  anode  tall  voltage,  bmoe  the  electric  field  near  the 


cathode  fall-positive  column  boundary  is  Uss,  when  Ward's  formula  is 
used,  this  would  also  result  in  a  smaller  electron  current  density 
entering  the  positive  column.  A  broader  anode  fall  region  would  then 
be  required  tor  the  electron  current  density  to  reach  the  total 
discharge  cm  rent  limit.  So  even  though  the  electric  field  g. adient 
would  also  he  less  in  the  anode  full  region,  the  lact  that  it  would  be 
over  a  lunger  distance  could  increase  the  anode  fall  voltage. 

The  anode  in  the  SHEATH  cd culations  was  placed  in  the  region  of 
the  negative  glow.  Any  correlation  in  the  data  close  to  the  anode  is 
i rideterminant  due  to  the  changing  location  of  the  anode  in  the  negative 
glow  region,  lhis  occurred  because  the  cathode  fall  width  changed 
while  the  electrode  separation  was  held  constant  for  each  rare  gas 
mixture . 

In  audition,  the  SHEATH  code  changes  from  constant  energy  bins  to 
constant  d'atial  bins  very  near  err  at  the  minimum  in  the  electric 
field.  At  present,  tins  change  of  variable  is  not  smooth  for  all 
parameters  calculated  by  the  cede,  so  any  predictions  by  the  SHEATH 
method  on  the  anode  fall  will  be  held  in  reservation.  Thus  in  their 
present  configuration,  neither  code  properly  describes  a  normal  anode 
fall  region. 

In  order  for  the  SliLAih  model  to  be  able  to  correctly  describe  a 
normal  anode  toil,  several  new  modifications  would  have  to  be  made. 

The  model  could  be  charged  to  reflect  the  boundary  conditions  at  the 
anode  instead  of  the  r.i'hode,  so  that  the  calculation  started  at  the 
anode  arid  preceded  towards  the  posit  ive  column. 
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This  also  requires  that  E  (ti‘ ,  j  (d),  and  the  electron  distribution 
striking  the  anode  be  known.  Potential  problems  with  the  point  at 
which  the  code  changes  from  equal  energy  increments  to  spatial  incre¬ 
ments  could  also  occur  again,  since  the  energy  gained  in  a  mean  free 
path  is  much  less  in  the  anode  fall.  Similarly  ttiere  would  he  poten¬ 
tial  problems  with  greater  than  or  less  than  statements  when  the  sign 
convention  was  changed  from  x  tc  -x. 

If  a  new  computer  system  is  installed  which  has  2-3  times  larger, 
accessible  memory  space,  it  would  also  be  feasible  to  expand  the  array 
so  that  discharges  with  larger  electrode  separations  could  be  modeled. 
This  is  perhaps  the  easiest  solution  to  modeling  the  anode  fall  using 
the  SHEATH  code. 


SEll JON  VI.  CONCLUSIONS  AND  SOME 
CONSIDERATIONS  FCR  FUTURE  STUDY 


1.  CONCLUSIONS 

In  the  present  study,  two  very  different  methods  of  analyzing  glow 
discharges  (GLOW  and  SHEATH)  were  developed  and  applied  to  three  dif¬ 
ferent  cases  represented  by  He/HCl  ,  Ar/HCl ,  anG  Xt/'HCl  mixtures.  Even 
though  the  two  methods  did  use  different  input  data  bases  founded  on 
the  prior  work  of  different  authors  using  different  experimental  ap¬ 
proaches,  the  results  of  the  two  methods  had  the  same  trends  and  both 
models  scaled  according  to  Equation  74.  The  following  conclusions 
are  based  on  the  results  presented  in  Sections  III,  IV,  and  the  compar¬ 
ison  of  the  models  in  Section  V: 

a)  Electrons  were  found  not  to  be  in  equilibrium  with  the  elec¬ 
tric  field  throughout  the  cathode  fall  region.  Close  to  the  cathode, 
the  whole  distribution  is  accelerated  up  to  an  energy  eoual  to  the 
ionization  potential.  This  predominantly  dark  region  is  called  the 
Primary  dark  space.  Next  the  electron  distribution  function  becomes 
birrodal  as  some  of  the  electrons  lose  their  energy  to  ionization  and 
excitation.  This  luminous  region  is  called  the  cathode  glow.  In  this 
region,  the  majority  of  the  electrons  have  a  mean  energy  greater  than 
the  ionization  eneryy.  As  more  and  more  of  the  electrons  lose  energy 
to  inelastic  processes,  the  cathode  qlow  merges  into  a  relatively 
darker  region  called  the  cathode  dark  space.  The  energy  difference 
between  the  decreasing  number  of  fast  and  the  increasing  number  of  slow 
electrons  gets  larger  as  the  fast  electrons  continue  to  bo  accelerated 
by  the  linearly  decreasing  electric  field.  As  the  electric  field  be¬ 
lts 


comes  very  snic.ll  and  almost  constant  in  the  negative  glow,  the  remain¬ 
ing  fast  electrons  resemble  an  electron  beam.  They  deposit  their 
energy  in  the  gas  creating  many  highly  excited  species  which  cascade 
down  and  produce  the  bright  luminescence  characteristic  of  the  negative 
glow  region,  ms  a  result  of  this  nonequilibrium,  the  GLOW  or  equilib¬ 
rium  model  was  found  in  Section  V  to  be  inadequate  in  describing  both 
the  microscopic  parameters  such  as  charged  particle  densities  and  ioni¬ 
zation  and  attachment  coefficients  as  well  as  macroscopic  parameters 
such  as  cathode  fall  voltage,  cathode  fall  length,  electron  current 
density  and  electric  field.  A  model  which  allows  the  electrons  to  not 
be  in  equilibrium  with  the  electric  field  such  as  the  SHEATH  model  or 
Monte  Carlo  technique  must  be  used  to  describe  the  parameters  in  the 
cathode  fall  region.  The  negative  qlow  region  should  also  be  modeled 
by  a  technique  which  allows  nonequilibrium.  The  model  for  the  negative 
glow  however  could  be  simplified  by  separating  the  electron  distribu¬ 
tion  function  into  a  low  energy  component  and  a  high  energy  component. 

b }  Better  agreement  between  theory  and  experiment  can  be  ob¬ 
tained  for  pure  gases  by  a.  judicious  choice  of  y  and  k^.  Equation 
74  can  be  used  to  relate  Y  end  k+  to  experimental  cathode  fall 
voltage,  currents,  ard  cathode  fall  widths.  These  values  ot  >  and  k+ 
should  then  be  used  as  inputs  in  the  SHEATH  model  since  the  SHFATH 
model  was  shown  in  Section  IV  to  obey  the  scaling  relationship  given  by 
F quo t ion  7b. 

r)  The  slope  ct  the  electric  field  in  the  cathode  fall  is  deter¬ 
mined  by  the  positive  ion  oersity.  If  Poisson's  equation  is  written  in 
terms  ot  the  ratios  ct  the  charged  particle  current  density  and  dritt 
velocity,  then  the  slope  of  the-  field  close  to  the  cathode  is  deter- 


mined  by  trie  mobil’ty  il  }  of  the  positive  lcrs.  The  positive  ion  cur¬ 
rent  density  is  constant  close  to  the  cathode  and  is  determined  by  t tie 
total  discharge  current  oensiLy  and  The  length  of  the  cathode  fall 
is  determined  L>  the  maintenance  condition  tor  the  discharge  involving 
the  multiplication  et  electrons  through  the  cathode  tail  (Equations  ] 
and  103).  The-  slope  or  the  electric  field  and  its  extrapolated  x-axis 
intercept  (d(  )  then  determines  t.he  voltage  drop  across  the  cathode  fall 
regioii . 

d)  The  formation  of  small  amounts  (n  <  ng)  of  negative  ions  in 
the  cathode  fall  region  was  found  not  to  he  the  reason  that  the  cathode 
sheath  length  contracts  in  the  axial  direction.  Actuary,  the  forma¬ 
tion  of  negative  ions  rauses  the;  cathode  tall  voltage  to  increase 
slightly  which  in  turn  results  in  a  slight  increase  in  the  cathode  fall 
length.  These  trends  could  he  seen  in  the  results  of  both  the  equilib¬ 
rium  and  nonequ i  1  ibrinirr  analyses  in  Sections  III  and  IV. 

e)  Whf  r  an  attaching  gas  which  has  a  low  ionization  potential  is 
added  to  he  or  he,  the  observed  i on traction  cf  the  cathode  fall  length 
was  found  to  be  due  to  an  increase  in  ionization  m  this  region.  This 
increase  in  ionization  is  a  result  of  the  lowering  of  the  ionization 
threshold  and  tie  difference  in  cross  section  shape  for  the  two  gases 
in  the  mixture. 

f)  The  bh LA T 1 1  method  which  allowed  the  electrons  to  not  be  in 
equilibrium  with  tht  tie  id  predicted  a  Townsend  attachment  coefficient 
which  contuinto  a  peat  at  >fe  oathocu-,  reached  a  minimum  at  the  region 
where  the  electron  i'i,m  energy  and  Townsend  ionization  peaked  and  then 
yie-v,  again  as  tht  tield  as  well  cis  the  electron  mean  energy  decreased. 

I  hut  t  he  cathrdr  tall  in  an  t  . ei  t r oneoa t i vn  yas  mixture  was  found  to  be 
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character ized  by  a  thin  region  cl  use  to  the  cathode  in  which  attachment 
dominated,  followed  by  a  rey.uri  in  which  ionization  of  HC1  dominated 
the  ionization  of  He  or  Ar  hut  not  Xe,  before  uniform  ionization  of  the 
mixture  Logan. 

g)  The  positive  ion  current  density  at  the  cathode  fa! 1 -negati ve 
glow  boundary  was  found  to  still  contribute  a  significant  amount 
(5-15/;)  to  the  total  current  density.  Thus  the  boundary  condition  that 
the  positive  ion  current  density  approaches  zero  at  this  boundary  is 
not  a  gooci  assumption. 

Although  this  analysis  was  mainly  concerned  with  modeling  phenome¬ 
na  in  the  cathode  tall  region,  the  conclusions  d  and  e  above  could  also 
be  applied  to  the  anode  tall  region.  In  contrast  to  the  cathode  fall 
region,  the  negative  chaige  density  is  much  larger  than  the  positive 
charge  density  in  the  anode  fall  region.  Thus  the  width  of  the  anode 
fall  would  depend  on  how  fast  the  electrons  are  accelerated  to  the 
ionization  potential.  Thus  the  relationship  of  the  electron  drift 
velocity  to  the  field  is  more  important  in  this  region. 

2.  THLORETIIAL  CONSIDERATIONS  FOR  FU11JRF  STUDY 

In  the  course  of  this  study  several  areas  can  be  identified  which 
warrant  further  theoretical  investigation: 

a)  Since  ionization  of  HU  caused  the  cathode  la'll  to  contract 
iri  He,  further  study  i  t  the  effects  resulting  from  changing  the  ioni¬ 
zation  cross-sect ioii  would  tie  of  value.  Both  changes  in  the  threshold 
anci  initial  slope  of  the  cross  sect-ic  r  should  be  examined  to  see  which 
affects  changes  ir  the  cathode  sheath  structure  ttie  most. 

b,-  /.  [arameter  napping  of  the  cpoondary  emission  coefficient  and 


positive  icr,  mobility  (1  '  is  recommended.  This  has  not  been  inves¬ 
tigated  using  noneqm 1 ihrium  techniques. 

t )  Dsing  the  present  results  from  the  cathode  fall  as  inputs, 
the  nonequi 1 ibrium  method  could  he  used  to  investigate  the  negative 
gl ow-Faraday  dork  space  region  provideo  the  field  does  not  go  negative. 

d)  The  linear  aria  lysis  presented  ir  this  stuoy  can  be  expanded 
to  describe  the  kinetics  in  different  electrode  geometries  represented 
by  thyratrens,  plasma  reactors  for  depositing  thin  films,  or-  hollow 
cathode  discharges. 

e)  The  results,  cf  this  study  can  also  be  applied  to  understand¬ 
ing  the  field  variations  and  charged  particle  kinetics  in  striations. 
Different  boundary  conditions  would  be  required  for  striations.  The 
charged  particle  densities  on  each  side  of  the  striation  would  replace 
the  boundary  conditions,  on  the  current  densities  at  the  electrodes. 

t)  Hot  cathodes  could  be  modeled  using  either  the  GLOW  or  SHEATH 
techniques.  Because  the  accelerating  voltage  is  less  in  the  cathode 
fall,  there  should  be  a  smaller  degree  or  nonequi 1 ibrium.  Thus  the 
GLOW  model  should  give  much  better  agreement  in  this  case.  In  addi¬ 
tion,  the  boundary  condition  at  the  cathode  relating  the  positive  ion 
to  the  electron  current,  density  would  have  to  be  replaced  by  a  single 
source  of  electron-.. 

y)  bue  tc  1  iii'”-‘-'d  compute!  storage  space,  the  only  electron  col¬ 
lision  fire. cesses  calculated  as  fence. ions  ol  distance  were  the  ionisa¬ 
tion  and  attachment  coef t  icients.  further  investigation  of  the  other 
excit.atice  coefficients  and  coefficients  lor  vibrations!  excitation 
would  bo  useful  to  comp  lei e  the  understanding  of  the  cathode  fall  re- 
g  i  on . 
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h)  Operation  of  t lit  SlttA'IH  cock  would  be  ninth  more  efficient  as 
well  as  convenient  if  a  numerical  technique  were  added  to  converge  to 
the  lowest  value  of  E  for  which  u  complete  solution  exists,  before  the 
field  is  iterated  tor  consistency  with  tht  electron  distribution  func¬ 
tion.  This  technique  could  he  similar  to  the  modified  interpolation 
technique  which  was  used  in  the  CLuW  code. 

l j  I  lore  detailed  iun  kinetics  options  would  make  the  SHEATH  model 
more  universal.  Electron-ion,  and  ion-ion  recombination,  detachment, 
and  two  step  irrigation  could  be  added  such  that  the  ion  or  metastable 
densities  could  be  updated  after  each  calculation  of  the  electron  dis¬ 
tribution  function  for  a  seif-consistent  solution.  Recombination  and 
detachment  will  be  more  important  in  the  negative  glow  and  Faraday  dark 
space.  These  processes  can  not  be  incorporated  in  the  present  GLOU 
mode  1  because  they  are  nonlinear  with  respect  to  j  . 

j)  Cathodes  can  be  coated  with  special  films  which  emit  negative 
ions.  Further  investigation  of  the  effect  of  a  cathode  as  a  source  of 
negative  ions  on  the  cathode  fall  voltage  and  width  would  be  interest¬ 
ing.  A  further  increase  in  voltage  arid  perhaps  width  of  the  cathode 
fall  should  be  anticipated  according  to  the  theory  in  Section  III  if 
negative  ions  are  released  at  the  expense  of  electrons. 

3.  EXPERIMENTAL  CONSIDERATIONS  I  UR  FUTURE  STUDY 

Several  experimental  topics  have  also  been  identified  in  the 
course  of  this  study  winch  warrant,  further  investigation.  At  present 
the  first  two  topics  are  beiro  pursued. 

a)  No  experimental  evidencr  could  be  found  for  the  non-contrac¬ 
tion  oT  the  cathode  ’all  in  mixtures  of  Ar  or  Xe  with  HC1  or  any  other 


attaching  gas.  A  smalt  discharge  lube  is  presently  under  eonstruction 
in  the  Plasma  Physics  Group  (Al  WAL/f  OCC-3)  which  will  allow  confirma¬ 
tion  of  this  phenomenon .  A  mass  • pet trometer  will  be  used  to  monitor 
any  change  in  the  neutral  particle  chemistry. 

b)  It  is  now  possible  to  consider  measuring  the  eiectric  field  as 
a  function  of  distance  through  the  cathode  tall  region.  Some  potential 
techniques  are  non  pc-r  tiubiiiy  boppler-f i  ee  laser  polarization,  Doppler- 
free  spectroscopy,  or  Stark  induced  mixing  of  different  parity  levels 
in  an  excited  molecular  electronic  stale. 

c)  The  energy  distribution  of  the  electrons  could  be  measured  in 
an  obstructed  discharge  by  sampling  the  electrons  through  a  snail  hole 
in  the  anoce.  These  resu'ts  could  then  be  compared  with  theory. 

d)  Similarly  the  energy  distribution  of  the  positive  ions  could 
be  sampled  through  a  small  hole  in  the  cathode.  Both  major  and  minor 
ions  should  he  sampled  because  of  the  differences  in  charge  exchange 
cross  sections,  by  knowing  their  energy  and  charge  transfer  rates, 
their'  urigin  in  the  cathode  fail  or  negative  glow  could  be  determined 
and  checked  with  theory. 

e)  Laser  induced  fluorescence  could  be  used  to  search  for  fast 
neutral  molecules  which  would  be  an  indication  that  charge  transfer  has 
occurred.  This  technique  would  compliment  those  described  in  c  and  d 
abovt  and  could  be  accomp I  i slice  simultaneously. 

it  is  hoped  tha.t  ihis  study  has  provided  c  better  understanding  ct 
the  Gathers  tall  ano  the  mtluenrr  oi  negative  ions  in  th 1 s  region. 

The  generic  if  Milts  of  this  study  nan  he  applied  to  txcimer  lasers, 
f hyratroris ,  or  plasma  reactors  tor  the  deposition  of  thin  films.  As  a 
final  guideline,  it  m  jest  as  important  to  understand  the  positive  ion 


kinetics  in  the  cathode  sheath  region  as  it  is  to  understand  the  e 
tron  kinetics,  because  it  is  the  ior,  mobility  which  determines  the 
slope  of  the1  electric  lie  Id. 
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APPENDIX  A:  REVIEW  OF  NEGATIVE 


ION  FORMATION 

In  order  to  pro  writ-  continuity  in  the  descriptions  of  the  models, 
t tie  mechanisms  and  important  feaiuros  or  negative  ion  formation  are  re¬ 
viewed.  Several  reactions  occur  which  form  negative  ions  and  these  are 
listed  and  described  below. 


X  -t  e  -*  X  '  h  v 

radiative  capture 

X  +  e"  +  M  -  X"  h  M 

3  body  electron  capture 

XY  +  e  -*  X  +  Y 

dissociative  attachment 

XY  +  e"  -* X+  +  Y"  -i 

ion  pair  formation 

(X  =  an  atom,  e~  =  an  electron,  h  =  a  photon,  M  =  an  electron  or  an 
atom,  XY  =  a  molecule  maae  of  atoms  X  and  Y) 

1.  DISSOCIATIVE  ATTACHMENT 

The  most  common  process  of  forming  negative  ions  is  dissociative 
attachment.  In  this  Case,  an  electron  collides  with  a  molecule  XY  and 
dissociates  the  molecule  into  a  neutral  atom  X  and  a  negative  ion  V  . 

To  aio  the  visualization  of  transitions  between  states,  dissociative 
attachment  can  be  divided  into  3  cases  usiny  molecular  potential  energy 
curves.  These  three  cases  v.ill  be  discussed  next.  The  discussion  will 
include  the  influence  of  vibrational  cxcitatior  on  dissociative 
attachment. 

In  all  three  cases  represented  in  Figure  F.8a ,  b,  c,  the  heavy 
solid  line  refers  to  the  ground  state  of  the  neutral  molecule  XY.  The 
nuclear  separation  oscillates  ttieri  between  P^  ard  R^.  The  thinner 


R,  R.  R.  R, 

DISTANCE  DISTANCE 

(b)  (c) 

Figure  58 

Potential  Energy  Curve  Allowing  Dissociative  Attachment 

curve  refers  to  the  molecular  ion  state  XY"  which  at  infinite 
separation  (X  +  Y~)  tends  to  a  limit  which  lies  below  the  potential 
energy  curve  of  the  neutral  molecule  by  an  amount  of  energy  equal  to 
the  electron  affinity  (E  )  of  atom  Y,  The  solid  region  represents  the 

O 

range  over  which  the  negative  ion  state  is  stable.  The  dashed  region 
represents  the  range  over  which  the  negative  ion  state  is  unstable.  In 
this  region  the  negative  ion  may  give  up  the  electron  and  return  to  a 
ground  state  neutral  molecule  XY,  via  autodetachment.  Ifr(R)/h  is  the 
rate  (sec  *)  for  autodetachment,  then  the  energy  in  this  unstable 
region  is  uncertain  by  an  amount  ^  (R)  and  the  dashed  line  actually 
should  be  a  band  with  a  width  v  (R).  The  molecular  ion  state  curve 
crosses  the  ground  state  curve  at  R  .  Applying  the  Franck-Condon 
principle  to  the  transition  from  the  ground  molecular  state  to  the 


negative  ion  state  leads  to  a  final  state  represented  by  some  point 
lying  between  a  and  b  on  the  pote.  tial  energy  curve.  1  he  threshold  foi 


dissociative  attachment,  E^,  is  the  difference  between  the  molecular 
dissociation  energy  D.  and  the  electron  a  finity,  F..  ,  of  the  negative 
atomic  ion.  Thresholds  have  been  tabulated  for  several  molecules  and 


they  are 


Er  A  =  D  -  E 
DA  o  a 


(106) 


compared  to  thresholds  for  ion  pair  formation  shown  in  Table  9. 


The  competition  between  autodetachment  and  dissociation  in 


molecular  negative  ions  leads  to  a  strong  dependence  of  the  cross 
section  on  the  mass  or  the  nuclei  and  the  rotational  and  vibrational 


states  of  the  neutral  molecule.  In  case  I  represented  in  Figure  58a, 
first  suppose  autodetachment  did  not  occur,  lf.cn  all  transitions  from 
the  ground  state  to  the  iuri  state  would  result  in  dissociation  into  an 


atom  X  and  ior;  Y~  with  a  relative  kinetic  energy  between  E ^  and  Ep.  If 
autodetachment  is  allowed  to  occur  at  some  internuclear  separation  R  at 
a  rate  ; (R)/h,  then  after  the  transition  tc  the  negative  ion  state, 
the  atomic  nuclei  begin  to  separate  with  a  relative  velocity  V(R),  and 
the  fraction  of  negative  ion  molecules  remaining  (Reference  b)  or 
"surviving"  is 


f  (R) 


(16/) 


If  the  sepai at  ;on  between  the  atomic  nuclei  rtuches  R  before  auto- 
detachmerit  takes  place,  then  the  neutral  molecule  dissociates  into 
X  +  Y  .  The  cross  section  for  this  process  can  then  be  approx ima tec  by 


Molecule 

HgT 

HgBr 

HqCl 
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Attachment  anci  Ion  Pair  Production 
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0 
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1 
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1 
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3.68 
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1 

2.08-S 

3.55 

1.47 
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10.43 

1 

.75-H 

2.80 

13.60-S 

1 

3.45-F 
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2.49 
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16.06 
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.75-H 
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17.42-F 

2 
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1 
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1 
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1 

1.27-C 
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2 
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2 

1.27-C 

9.82 
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’a(R)'"cexp-  „/'  S  ® 

where  11  is  the  cross  section  tor  the  initial  capture  (Reference  8). 
Because  heavy  atomic,  nuclei  typically  oscillate  more  slowly  than  light¬ 
er  nuclei,  the  heavier  nuclei  spend  more  time  near  the  turning  points 
and  have  steeper  potential  wells.  This  increases  the  probability  that 
the  internuclear  separation  R  ot  the  molecular  ion  will  exceed  be¬ 
fore  it  autodetaches.  Typical  peak  cross  section  values  from  the 

ground  vibrational  state  tor  dissociative  attachment  range  from 

-20  2  14 

10  cm  for  light  molecules  such  as  H„,  HD,  end  Dn  to  10  or 

C  L 

-15  2 

10  cm  for  heavy  molecules  such  as  HI,  and  I^.  (See  Table  10). 

In  case  II  represented  iri  Figure  58  (b),  only  electrons  with  ener¬ 
gies  between  Ej  and  E^  can  dissociate  the  molecule  ana  produce  atomic 
negative  ions  through  the  dissociative  attachment  process.  The  resul¬ 
ting  negative  ion  Y  and  atom  X  thus  will  have  relative  final  energy 
ranging  from  0  to  E,-E^.  Electrons  with  energies  between  and  E^  can 
be  captured,  but  they  will  form  excited  molecular  negative  ions  XY~. 
These  molecular  lers  will  eventually  autodetach  unless  the  excess  ener¬ 
gy  can  be  absorbed  through  a  collision.  Thus  the  total  attachment 
cross  section  for  this  case  will  be  independent  of  pressure  for  pres¬ 
sures  greater  than  some  critical  pressure,  but  will  be  proportional  to 
the  pressure  for  pressures  less  than  the  critical  pressure.  The  criti¬ 
cal  pressure  occurs  when  the  average  time  tor  the  vibrationally  excited 
molecule  to  transfer  its  excess  energy  to  a  gas  molecule  is  comparable 
to  the  time  for  au todetachment . 
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Table  10.  Thresholds  and  Peak  Cross  Section 
Values  for  Dissociative  Attachment 


Process 

Diatomic  Threshold  Peak  Energy  at 

Molecules  (eV)  (cm2) _ Peak  (eV)  Reference 


Table  10.  (cont.) 


Process 

Diatomic  Threshold  Lneray  At 

Mo  I  ecules _ (eV ) _ Peak  (cm2) _ Peak  (eV) _ Reference 


In  case  III  represented  in  Figure  58c,  negative  ion  molecules 
with  vibratioral  states  v-a'  are  stable,  but  those  with  >>.•  are 
unstable  and  can  result  in  autodetachment.  These  states  may  be 
stabilized  through  collisions  similar  to  the  process  described  for  case 
II.  Similarly,  the  cross  section  for  attachment  will  exhibit  resonance 
maxima  which  are  associated  with  the  vibrational  levels  with 
provided  the  lifetime  of  the  excited  vibrational  state  is  long  compared 
to  its  period.  This  situation  occurs  in  CL  and  NO,  resulting  in  the 
production  of  0  v  and  NO”  after  collisions  with  other  gas  molecules. 
Table  10  tabulates  experimental  data  for  thresholds  and  cross  section 
peaks  for  several  attaching  molecules. 

Increased  vibrational  excitation  enhances  the  Dissociative 
attachment  cross  section  for-  all  three  cases  since  electron  capture  can 
generally  occur  at  larger  and  larger  values  of  R.  Thus  the  molecular 
negative  ion  is  more  likely  to  "survive",  the-  closer  R  approaches  Rc 
through  vibrational  excitation.  Rotational  excitation  can  also  enhance 
the  cross  section  for  dissociative  attachment  due  to  the  centrifugal 
stretching  occurring  in  the  excited  states,  but  in  a  less  dramatic 
manner  than  vibrational  excitation.  Parasley  and  Wadehra  (References 
17,  91)  have  shown  that  the  magnitude  of  the  dissociative  attachment 
cross  section  from  the  =  5  level  in  and  increases  three  to 
four  orders  of  magnitude.  Allan  and  Wong  (Reference  la)  have  measured 
similar  increases  in  their  experiments  with  HC1 .  The  change  iri  the 
"survival  factor"  is  thought  (Reference  91)  tc  be  more  important  than 
the  change  in  the  Franck  Condon  factor  in  explaining  this  dramatic  ef¬ 
fect.  if  tfe  lifetime  of  the  vibrational!'/  excited  molecule  is  long 
compared  to  the  period  of  vibration,  the  attachment  probability  will 
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contain  resonance  peaks.  This  is  a  resell,  uf  the  initial  electron 
capture  cross  section  exhibitina  resonance  peaks  due  to  the  ratios 

of  populations  in  the  various  vibrational  levels.  Experimentally  this 
variation  with  respect  to  rotational  and  vibrational  population  appears 
as  a  temperature  dependence. 

The  remaining  processes,  including  radiative  capture,  three  body 
capture,  and  ion  pair  formation  are  generally  less  important  than 
dissociative  attachment  and  will  be  discussed  only  briefly. 

2.  RADIATIVE  CAPTURE 

Radiative  capture  is  the  simplest  process  +  or  creating  negative 
ions.  Electrons  are  captured  directly  by  neutral  atoms  with  the  excess 
energy  being  dissipated  in  radiation.  The  amount  of  energy  released  by 
this  process  is  equal  to  L  t  F.  where  E  if.  the  kinetic  energy  of  the 
incident,  electron  and  E  is  the  electron  affinity  of  the  atom.  This 
capture  process  produces  a  continuous  emission  spectrum  extending 
i ndt  finitely  tron  the  long  wavelength  limit  at 

In:  (109) 

E 

a 

-Zl 

A  lypiial  magni'ude  tor  this  process  is  on  the  order  of  10  cm  with 

a  Inrn  pt.dk  near  iurn  impact  energy  (Reterence  8).  This  is  two  orders 
of  i ’ctiji  i  tude  .i’<i  Mer  than  the  smallest  cross  section  lor  dissociative 
attachment.  !(n  relatively  large  low  energy  peak  recurs  because  the 
elect rui!  spends  a  longer  tine  within  the  range  of  tie  atomic  field  at 
lewer  incident  energies. 
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THREE  BODY  ELECTRON  CAPTURE 


In  three  body  electron  capture  the  energy  released  by  the  process 

is  converted,  not  into  a  photon,  but  into  kinetic  energy  of  the  third 

body.  The  magnitude  of  the  attachment  cross  section  for  this  process 

can  be  determined  by  arguments  based  on  two  body  collisions.  The  third 

body  will  be  effective  in  absorbing  energy  only  if  it  makes  an 

effective  collision  with  the  capturing  atom,  when  the  electron  is  also 

-8 

within  an  atomic  radius  (r  ~  1C  cm)  of  the  capturing  atom.  If  \  is 
the  mean  tree  path  for  effective  collisions  between  the  atom  and  a 
single  third  body,  then  the  probability  for  an  effective  collision  is 
to  a  good  approximation  r  />.  .  This  results  in  an  effective  three  body 
attachment  cross  section  that  can  be  approximated  by  the  following 
expression  in  terms  of  an  effective  cross  section  per  unit 
concentration  of  the  third  body. 


?  r 
1  o 
a  '  r  y~ 

0  A  (no) 

Thus  it  remains  to  determine  x  ,  which  varies  considerably  depending  on 
the  nature  of  the  third  body.  Tv/o  cases  will  be  considered: 

1)  when  electrons  are  third  bodies,  and 
Z)  when  atoms  and  molecules  are  third  bodies. 

When  electrons  are  the  third  body,  the  effective  cross  section  for 
energy  transfer  is  on  the  order  10'  cm,  and  the  mean  free  path  for  a 


unitary  third  bedy  density  simply  becomes 

=  1016cni 

and  the  th i rd  booy  attachment  cross  section  is  on  the  order  of 

O  -  10"4(W 

a 


(in) 


dll') 
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In  comparison  to  dissociative  attachment,  an  electron  number  density 

20  3 

greater  than  or  eoual  to  10  cm  would  be  required  to  produce  a 
comparable  number  of  negative  ions  per  sec.  For  most  purposes  this  is 
an  exceptionally  large  electron  density. 

khen  atoms  or  molecules  are  the  third  body,  tv/o  cases  exist: 
resonant  and  non-resor.ant  energy  transfer.  In  the  resonant  case,  the 
energy  transferred  is  absorbed  totally  as  internal  energy  iri  the  third 
body.  The  cross  section  for  this  process  is  known  to  be  quite  large  or 
the  order  10  cm‘  or  greater(keference  8).  This  yielos  a  mean  free, 
path  and  effective  attachment  cross  section  of 

A  -  1014cm  aa  =  10“38cm2  (113) 

18  -3 

Thus  third  body  densities  of  10  cm  or  larger  would  be  necessary  to 
make  three  body  attachment  comparable  to  that  of  dissociative 
attachment.  Thus  molecules  possessing  several  internal  degrees  of 
freedom  are  the  most  effective  as  third  bodies.  However,  energy 
transfer  into  vibrational  motion  often  is  not  as  efficient  as  transfer 
into  rotational  motion. 

4.  ION  PAIR  FORMATION 

The  last  process  to  be  discussed  which  forms  negative  ions  is  ion 
pair  formation  or'  polar  excitation.  During  this  type  of  collision,  an 
electron  excites  a  molecule  (XY)  to  an  unstable  electronic  state  which 
dissociates  spontaneously. 

XY  +  e"  -  X4  +  Y~  +  e"  (114) 

This  process  differs  from  the  previous  processes  in  that  the  electron  is 
not  captured  during  the  collision,  but  corresponds  to  an  inelastic 
collision  in  which  the  electron  acts  as  a  source  of  energy  for 
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electronic  excitation.  Thus  polar  excitation  occurs  at  much  higher 
electron  energies  than  the  capture  processes.  Minimum  threshold 
energies,  E^  usually  range  from  10  to  25  eV. 

A  typical  potential  energy  curve  for  this  process  is  depicted  in 
Figure  59 


X+Y* 

X+Y 


R.  R, 

DISTANCE 

Fig  59.  Potential  Energy  Curve  For  Ion  Pair  Formation 

Again  the  heavy  solid  line  refers  to  the  ground  state  of  the  neutral 
molecule.  The  thinner  line  refers  to  an  unstable  electronic  state  which 
dissociates  spontaneously  into  ions.  The  energy  at  infinite  nuclear 
separation  tends  to  a  value  of  E.  -  E  ,  where  E.  is  the  ionization 

la  1 

energy  of  X  and  Efl  is  the  electron  affinity  of  Y.  The  final  relative 
kinetic  energy  (T)  of  the  ions  is  given  by  the  following  equation: 

T  ■  Eth  *  Ea  -  Ei  '  Do  <  115  > 

where  Dq  is  the  dissociation  energy  of  the  ground  state.  Thresholds  for 

ion  pair  formation  for  several  diatomics  are  previously  listed  in  Table 

9. 

The  shape  of  the  cross  section  for  polar  dissociation  differs  from 
the  previous  capture  process.  Actually  it  resembles  the  shape  of  other 
electronic  excitation  processes  and  thus  is  typified  by  a  sharp  rise 
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near  threshold  to  a  maximum  on  the  order  of  10  at  2  to  4 

times  the  threshold  energy,  after  which  it  declines.  The  decline  being 
inversely  proportional  to  the  eneryy  or  more  slowly  if  optical 
transitions  are  allowed  between  the  ground  state  and  the  excited  state 
(Reference  68:497). 


APPENDIX  B.  PARAMETERS  FOR  GLOW 


CODE  CALCULATIONS 


This  appendix  describes  the  parameters  used  to  characterize  the 
ionization,  attachment,  diffusion,  and  ion-ion  recombination 
coefficients  as  well  as  the  ion  drift  velocity.  Most  of  the  data  for 
the  rare  gases  was  taken  from  the  literature.  Effort  was  taken  to 
obtain  the  best  parameters  for  HC1  from  drift  tube  measurements  of 
Davies  (Reference  33:31-55). 

The  parameters  for  the  rare  gases  are  summarized  in  Table  11.  A, 
B,  and  s  refer  to  the  parameters  in  Equation  41  and  were  taken  from 
Ward  (Reference  95:2790).  The  transverse  diffusion  coefficient  was  ob¬ 
tained  from  data  tabularized  by  Dutton  (Reference  39:6491  (He),  651 
(Ar),  654  (Xe)).  The  mean  free  path  a  was  estimated  using  the  formula: 

me) 

where  o  is  the  momentum  transfer-  cross-section.  Momentum  transfer 
cross-sections  have  been  tabulated  by  Kieffer  (Reference  54:2(he), 


Table  11.  Parameters  Used  to  Describe  Rare  Gas  Processes 

in  the  GLOW  Code 


cm^ton 


cnwtorr' 


cm  torr 


29.12 


cm  torr' 


26.64 


sec 

i — nTTfT 

i  3. tit 


volt*sec 


S.52E-3  8.25E3 


36.00 


1.64E6 


5.E-3 


4.0L3 


13  (Ar),  17  (Xe)).  The  value  of  the  momentum  transfer  cross  section  at  leV 
was  used.  This  energy  is  representati ve  ot  electrons  diffusing  racially  in 
the  positive  column  region  ot  a  discharge.  The  value  of  r  was  also  taken 
from  the  data  given  by  Hard  (Reverence  92:2791)  for  v+  and  v  . 

The  ion-ion  recombination  rates  were  estimated  using  the  approxi¬ 
mate  scaling  formula  developed  by  Hickman  (Reference  49:4875).  The  formula 


where  T  is  the  gas  temperature  in  degress  Kelvin,  m  is  the  reduced  mass 
of  the  diatomic  molecule  in  a.u.,  and  E.A.  is  the  electron  affinity  of 
the  negative  ion  in  eV.  The  values  of  these  parameters  and  the  result¬ 
ing  ion-ion  recombination  rate  for  a  3CC°K  gas  is  given  in  Table  12. 


Table  12.  Parameters  for  Calculating  Ion-Ion  Recombination  Rates 


m 

(a.u. ) 

E.A.  (CD 

(eV) 

Kr 

(cm/ sec) 

He 

4. 

3.613 

7.2E-6 

Ar 

39.90 

3.613 

3.15E-6 

Xe 

331.3 

3.613 

2.59E-6 

The  ionization  and  attachment  coefficients  for  HC'i  were  curve  fit 
from  the  experimental  data  of  Davies  (Reference  33:33).  Using  Equation 
41,  both  values  ol  s  -  1  and  s  =  1  were  tried.  The  value  s  =  1  gave 
the  best  fit,  which  also  yieldea  values  for  A  =  14.34  and  B  =  243.94. 

In  order  to  model  the  attachment  coefficient  over  the  complete  range, 


both  Equations  42  and  43  had  to  be  used.  A  fifth  order  (n=5)  polynomi¬ 
al  was  used  for  E/p  <  90v/cm-torr  v/ith  the  following  values  for  coeffi¬ 
cients 


a  =  -.748878814, 
a?  =  .1274582651, 
ah  =  -.0052731662, 
ah  =  9.48227228E-5, 
a]?  =  -7.86191615E-7 ,  and 
aZ  =  2.46552683F-9. 

For  E/p  >.  90v/cm-torr,  ihe  best  f't  was  given  by  C  =  .254193623029, 
D  =  9. 14964650482E-3,  and  t=l.  This  concludes  the  data  inputs  used  for 
the  GLOW  code  except  for  the  discharge  parameters,  such  as  current, 
pressure,  and  electrode  separation  which  are  discussed  in  the  Section 
III  with  the  results  of  the  calculations. 


APPENDIX  C.  PARAMETERS  FOR  SHEATH 


CODE  CALCULATIONS 

This  appendix  describes  the  parameters  used  in  the  SHEATH  code. 

The  initial  values  of  code  parameters  will  be  summarized  first  followed 
by  a  description  of  the  electron  iitipact  cross-sections  used  in  calcula¬ 
ting  the  electron  distribution  function. 

The  discharge  parameters  are  summarized  in  Table  13.  The  anode 
was  placed  just  outside  the  cathode  fall  region  where  the  electric  field 
leveled  off  to  a  constant  value.  This  was  done  to  keep  the  energy  step 
size  as  small  as  possible,  less  than  1.6eV,  in  order  to  maximize  the 
energy  resolution  in  calculating  the  distribution  function.  With  this 
step  size  limitation,  lengthening  the  electrode  distance  was  not  feasi¬ 
ble  since  most  of  the  storage  space  on  a  Cyber  74  was  already  being 
used.  For  100%  He,  the  electrode  distance  was  2.2  cm  and  for  the  re¬ 
maining  He/HCl  mixtures  the  electrode  distance  was  2.1  cm  due  to  the 
contraction  that  accurred.  J  is  the  discharge  current  density  and  Y  is 
the  secondary  emission  coefficient.  The  distribution  of  electrons 


Table  13.  Discharge  Parameters  Used  in  the  SHEATH  Lode 


Mixtures 

D 

(cm) 

J 

(  amp/cmr) 

Y 

He/HCl 

2. 1-2.2 

16 

.2 

Ar/HCl 

1.15 

10 

.0417 

Xe/HCl 

.87 

10 

.004 
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leaving  the  cathode  was  the  same  as  that  used  by  Lone  (Reference  60) 
and  is  displayed  in  Figure  60.  The  x  axis  is  plotted  in  terms  of  ener¬ 
gy  bins  which  were  typically  1.  -  1.6eV  in  width  depending  on  the  gar- 
mixture.  This  distribution  is  higher  in  energy  than  the  one  used  by 
Boeuf  et  al .  (Reference  14)  which  was  a  flat  isotropic  energy  distribu¬ 
tion  over  the  range  of  0  to  5eV.  Boeuf  and  Marode  (Reference  13)  later 
used  a  narrower  distribution  which  was  isotropic  and  flat  between  4  and 
5  eV. 

The  gas  parameters  are  summarized  in  Table  14.  The  gas  density 
was  constant  threughtout  the  discharge  and  was  determined  from  the  gas 
pressure,  pQ  and  temperature,  T  .  The  positive  ion  mobility,  k+  was 
obtained  from  Ward  (Reference  92:2791)  and  was  the  same  as  that  used  in 
the  equilibrium  analysis.  Qn]_x  determined  the  point  at  which  the  code 
transforms  from  energy  space  to  spatial  coordinate  as  described  earlier. 
It  was  varied  so  the  coordinate  transition  was  as  close  as  possible  to 
the  minimum  of  the  electric  field  within  the  bounds  of  Equation  88. 


Table  14.  Gas  Parameters  Used  in  the  SHEATH  Code 

po  To  k+ 

cm  -torr 

Mixtures _ (torr) _ (°K) _ volt  sec 

/HC1  1  300  4.1E+4 

/ HC 1  1  300  8.25E+3 

/HC1  1  300  4E+3 


The  cross  sections  for  He,  Ar,  Xe,  and  HC 1  which  were  used  in  the 
EATH  code  calculations  are  discussed  next.  The  most  important  cross 
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sections  in  the  calculations  are  the  ionization  cross  sections  of  the 
rare  gases.  This  cross  section  determines  how  quickly  the  electrons 
multiply,  and  thus  directly  effects  nfc  and  dE/dx  in  Poisson's  equation. 
For  the  rare  gases,  the  ionization  cross  section  was  obtained  from 
measurements  by  Rapp  et  al .  (Reference  7b;  over  the  complete  range. 

The  ionization  cross  section  for  HC1  had  to  be  extended  above  IGOeV. 
Since  HC1  was  used  in  small  concentrations  and  the  mean  electron  energy 
in  the  HC1  mixtures  was  always  less  than  87  eV,  any  error  in  the 
approximation  should  be  negligible,  hore  specifically  the  peak 
electron  mean  energy  ranged  in  He  mixtures  from  59eV  to  99eV  in  pure 
He,  for  Ar  mixtures  it  ranged  from  20.5eV  to  21. eV,  and  for  Xe  mixtures 
it  ranged  from  25.4eV  to  25.9eV. 

The  six  cross  sections  used  to  model  the  electron  kinetics  in  He 
are  plotted  in  Figure  61.  The  threshold  ana  references  for  each  range 
are  listed  in  Table  15  below.  The  electronic  cross  sections  were 
extended  using  the  formulas  listed  under  the  reference  column. 

The  cross  section  for  Ar  are  plotted  in  Figure  62  and  the  refer¬ 
ences  and  thresholds  ore  given  in  Table  16.  The  three  electronic  cross 
sections  were  derived  using  a  standard  technique  of  backing  out  cross 
section  from  transport  data  using  a  coae  which  solves  the  equilibrium 
Boltzmann  equation.  Ratios  between  these  cross  sections  were  obtained 
from  data  of  Eggarter  (Reference  41}),  Ganas,  etal  (Reference  44),  and 
Ghutjan  (Reference  28).  keeping  the  ratios  constant,  the  sum  of  the 
three  cross  sections  was  ther  varied  and  the  predicted  electron  crift 
velocity  and  Townsend  ionizction  rate  v.rre  compared  with  transport  data 
listed  in  Dutton  (Reference  29}  and  Specht,  et  al.  (Reference  83).  Up 
to  2UV,  the  agreement  with  transport  data  is  within  9V  and  tire  total 
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Momentum 

Transfer 

0. 

0.-6. 

6.-500. 

Kieffer  (Reference  54) 
Hayashi  (Reference  48) 

Electronic 

3  1 

rs  +  rs 

19.805 

0.-23.54 

24.-100 

100.-500. 

Kieffer  (Reference  54) 
Constant 
°  =  .624/ > 

3 

cP 

20.949 

0.-199.6 

200.-500. 

Kieffer  (Reference  54) 
o  =  .02755/e5 

?!p 

21.203 

0.-196.1 

200.-500. 

Kieffer  (Reference  54) 
o  =  1.179/e 

3XP 

23.071 

0.-207.8 

210.-500. 

Kieffer  (Reference  54) 
o  =.3605/ e 

Ionization 

24.586 

0.-500. 

Rapp,  et  al . (Reference 
75) 

Table  16.  Electron 

Impact  Cross 

Sections  in  Ar 

Process 

Threshold 

(eV) 

Range 

(eV) 

Reference 

Momentum 

Transfer 

0. 

0.-30. 

30.-500. 

Kieffer  (Reference  54) 
Long  (Reference  60) 

Electronic 

"P3/2<3P5 

)  11.6 

0.-500. 

Derived 

?P1/2^P6 

)  11.8 

0.-500. 

Derived 

Remaining 
Llec tonic 

13.2 

0.  -500. 

Derived 

Rapp,  et.  al .  (Reference 
75) 


or  izatiori 


15.759 


0.-500. 


electronic  cross  section  agrees  well  with  Jacob  and  Mangano  (Reference 
53)  up  to  17eV  where  their  cross  section  stopped.  The  shape  of  the 
total  electronic  cross  section  for  higher  cross  energies  was  obtained 
from  de  Heer,  et  al .  (Reference  35). 

The  cross  sections  for  Xe  are  plotted  in  Figure  63  and  the  refer¬ 
ences  and  thresholds  are  given  in  Table  17.  7':  for  Ar,  the  electronic 
and  a  section  of  momentum  cross  section  for  Xe  were  derived  by  backing 
them  out  of  transport  data. 

Table  17.  Electron  Impact  Cross  Sections  in  Xe 


Process  Threshold  Range  Reference 

(eV)  (ev) 

Momentum  CT!  0.-6. 5  k’ieffer  (Reference  54) 

Transfer  6.5-20.  Derived 

20.-500.  Hayashi  (Reference  48) 

Excitation  8.32  0.-500.  Derived 

Ionization  12.13  0.-500.  Rapp,  et  al . (Reference 

75) 

The  momentum  transfer  cross  section  was  lowered  in  the  region  from 
about  6.5  to  29eV  in  order  to  get  the  drift  velocity  to  agree  within 
20%  of  that  measured  by  Huang  (Reference  50)  and  the  data  listed  in 
Dutton  (ref  39).  Huang  found  that  the  electron  mobility  for  fields 
above  4Td  remained  constant,  given  by  uN  =  5.5E21 

molecules/(v-sec-cm)*.  The  shape  of  the  electronic  cross  section  below 
lleV  was  taken  from  Dowell's  unnormalized  data  (Reference  37).  The 
ionization  rate  predicted  from  the  Boltzmann  code  agreed  within  20% 
with  the  data  in  Dutton  (Reference  39)  up  to  60Td  using  this  electronic 
cross  section  up  to  20eV.  The  shape  at  higher  energies  was  obtained 
from  Ganas,  et  al.  (ref  44)  and  de  Heer  et  al.  (ref  35)  and  normalized 
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to  the  value  derived  for  20eV. 

The  cross  sections  for  HC1  are  plotted  in  Figure  64  and  Figure  65. 
The  thresholds  and  references  for  each  range  are  listed  in  Table  18. 

The  electronic  cross  section  for  the  b  and  C  state  was  assumed  to  have 
a  peak  at  lOOeV  which  is  slightly  less  than  the  energy  of  the  peak  of 
the  ionization  cross  section.  The  cross  section  was  assumed  to 
decrease  according  to  the  formula  0  =  18. /zK  The  shape  of  the 
ionization  cross  section  from  Compton,  et  al.  (Reference  79)  was  used 
over  the  range  100-280eV  and  normalized  to  Davies'  data.  At  higher 
energies  the  shape  of  the  ionization  cross  section  for  0^  (Reference 
75)  was  used  which  also  peaks  at  the  same  energy.  It  was  normalized  to 
the  HC1  ionizaton  at  280eV.  This  concludes  the  discussion  of  the  data 
inputs  used  for  the  SHEATH  code. 


Table  18.  Electron  Impact  Cross  Sections  in  HC1 


Process 

Threshold 

(eV) 

Range 

(eV) 

References 

Momentum 

Transfer 

0. 

G.  -100. 
100.-500. 

Davies  (Reference  33:45) 
Loqrithrimic  e.  u -polaticn 

Vibration 
v  =  1 
v  =  2 

.36 

.70 

0.-7. 

0.-7. 2 

Davies  (Reference  33:46) 
Davies  (Reference  33:47) 

Attachment 

C1~ 

H" 

.67 

5.6 

0.-2. 7 

0.-12. 

Davies  (Reference  33:47) 
Davies  (Reference  33:48) 

Excitation 

A 

B  +  C 

5.5 

9.3 

0. -11. 
0.-100. 
100.-50C. 

Davies  (Reference  33:47) 
Davies  (Reference  33:48) 
Derived 

Ionization 

12.74 

0.-100. 

100.-500. 

Davies  (Reference  33:48) 
Derived 
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